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DIGITAL COMPUTATION FOR STIFFNESS MATRIX ANALYSIS 


J. S. Archer! J.M. ASCE 
(Proc. Paper 1814) 


ABSTRACT 


A matrix method is described for the static stress and normal mode 
analysis of highly indeterminate structures on large scale digital computers. 
The application of the technique to the analysis of basic types of Civil Engi- 


neering structures is discussed. A simple building frame problem is solved 
for stresses and deflections. 


SUMMARY 


The stiffness matrix method of structural analysis has been coded for use 
on high speed digital computers by several investigators. This technique has 
proved to be very successful for the analysis of highly indeterminate aircraft 
structures. It is here proposed that the same technique be used for the an- 
alysis of Civil Engineering structures on high speed digital computers. This 
may be accomplished through the development of similar digital computer 


programs or through the midification of digital computer programs now 
available. 


INTRODUCTION 


Scope 


This paper contains an energy development of the basic theory for the 
Stiffness Matrix Method of structural analysis. It is shown that the method is 
‘ particularly adaptable to a systematic analysis of any structure which may be 


described as a composition of simple structural elements, such as beams and 
plates. 


Note: Discussion open until March 1, 1959. To extend the closing date one month, a 
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part of the copyrighted Journal of the Structural Engineering Division, Proceedings 
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The present application of the technique to the analysis of aircraft struc - 
tures at Convair, Fort Worth, Texas and at Convair, San Diego, California is 
described. The application of the technique to the analysis of Civil Engineer - 
ing structures is also discussed. An example is given for a Civil Engineering 
structure to illustrate how the procedure works in a typical case. Finally 
some suggestions are given for a practical method of developing a similar 
procedure for individual needs. 


History 


The stiffness matrix type of structural analysis was discussed in 1944 in 
a paper by G. Kron,(1) Consulting Engineer, General Electric Company, 
Schenectady, N. Y. No practical application of the technique using digital 
computers was reported until December 1951 by H. U. Schuerch, (2) formerly 
associated with Convair, San Diego, California. S. Levy, now associated with 
the General Electric Company, carried out an investigation, (3) of the stiffness 
matrix procedure which was published in July 1953. G. Best and M. P. 
Keating, of Convair, Fort Worth, Texas, described the technique based upon 
their experience with the IBM 701 in a report,(4) dated September 1953. 

D. Williams, of the Royal Aircraft Establishment, England, prepared a very 
comprehensive hracions ~ on the subject which was published in June 1954. 
The technique of using matrix algebra for the analysis of aircraft structures 
has been taught at MIT by Prof. R. L. Bisplinghoff for several years. His 
text was first published in 1955.(6) 

The author first made use of the stiffness matrix technique at MIT in the 
analysis of a load frame for dynamic tests on reinforced concrete and brick 
walls. The analysis was coded and solved on the WHIRLWIND Computer in 
1954 and 1955. After coming to Convair, Fort Worth, Texas, in 1955 the 
author worked with G. Best to improve the IBM 701 analysis. The analysis 
has since been recoded for use on the IBM 704, 


Theory 


Let U represent the strain energy stored within a structure which is loaded 
by the forces pj, p2, ----, Pj, ----, Py» The temperature of the material re- 
mains constant and the supports are rigid. Applying Castigliano’s Theorem 
of Structure Equilibrium (Theorem I) one obtains; 


where 6; is the deflection of the point of application of the load pj in the 
direction of pj. 


If the strain energy is evaluated in terms of the loads p; acting upon the 
structure we may expand Equation (1) as follows: 


= = 2 (2) 


If one assumes that the structure is elastic then Castigliano’s Theorem for 
Linear Structures (Theorem II) may be applied. 
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DIGITAL COMPUTATION 
5. = dp (3) 
= 


Substituting Equation (3) into Equation (2), 


(2p /26,) (4) 
+ 


The partial derivative Op / I$, represents the force developed at 
a 
point j due to a unit deflection of point i, all other points remaining fixed. 
This force is represented by the symbol Sji- The subscript j represents the 


point at which the force acts and the subscript i the point at which the unit 
deflection is imposed. With this substitution Equation (4) becomes: 


(5) 
444 


From Betti’s Law or the generalized Maxwells’ Law of Reciprocal Deflec- 
tions we obtain the relation 


(6) 


and hence 
= 


Writing Equation (7) in its expanded form 


+ + + %,,5,, (8) 


It is evident from the expanded form that Equation (7) is a superposition 
equation expressing the total load at point i as the sum of the loads developed 
by each deflection component 6; acting by itself. Each portion of Equation (7) 
describes an independent component of the structural behavior. The compon- 
ents may represent translation or rotation. The total number of components 
is the number of degrees of freedom which the idealized structure possesses. 
Using matrix algebra notation, Equation (7) may be rewritten as 


p= S68 | 


where p is a vector or column matrix made up of the load components pj, po, 
-++,»Pi,+++,Ppy- 6 iS a vector made up of the deflection components 61, 


Q,-++, 5j,+-+.+, 6). Sis a square matrix consisting of an ordered array 
of the stiffness influence coefficients sjj of Equation (7). Matrix S is called 
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the stiffness matrix of the structure. In the expanded form Equation (9) ap- 
pears as follows: 


5 


4 
4 


; (10) 
Ong 
An, he 
Equation (9) may be manipulated similar to an ordinary algebraic equa- 
tion. For instance it may be solved for the deflections by premultiplying 
both sides of the equation by the inverse of the stiffness matrix S. 


5 = (11) 


The inverse of the stiffness matrix is represented by the symbol F and is 
called the flexibility influence coefficient matrix. It is composed of an 
ordered array of the flexibility influence coefficients fjj. In the expanded 
form it appears as follows: 


n 


In the typical problem Equation (11) must be constructed and solved for 
the deflections in terms of the loads. The direct determination of the coeffi- 
cients fij of Equation (11) is difficult and impractical for a large indetermin- 
ate structure. However, the coefficients 8;j of Equation (9) are readily cal- 
culated. The usual procedure is thus, to compute the stiffness matrix 
coefficients directly from the known properties of the individual elements of 
the idealized structure. The coefficients of the flexibility matrix are obtained 
through the mathematical operation of inverting the stiffness matrix. 

After the deflections of the structure are found the stresses developed in 
the structure may be determined. This is accomplished by using the stiffness 
matrix of the individual elements of the structure. Using computed deflections 
and equations of the same type as Equation (9) the internal loads on individual 
portions of the structure are readily obtained. Having obtained the internal 
forces on the individual elements, the structures engineer can readily com- 
pute the stresses in the structure. 


| 


anbio| (Q) 


vA 
= 
oO 
< 
Q 


vee | 


quowBes wereg 


LNIWI13 IWNGIAIGNI 


ASCE 1814-5 
] 
DW 
== 
WWW 
=a 
wy WWW 
| 
UW) 
oD 
= 
Qa. 
3 
o™ 
5 
' 


October, 1958 


[eNpIAipuy uo 
$22404 (2) 


dx 4 aynduioy (p) 


ajajdwioy 404 


JYNGIIOYd XIMLWN SSANSSILS 


ST 6 
“ y 
wo 
a 
aa 
‘9 
- ~ & 
wow 0 ww 
| 
< O 
vt 
te 
A 
YW) 


ASCE DIGITAL COMPUTATION 1814-7 


Present Applications 


IBM 701 


Structural analyses have been performed on the IBM 701 at Convair, Fort 
Worth, Texas, for structures having up to 106 independent coordinates. The 
structures analyzed were idealized as a group of interconnected structural 
elements. These elements consisted of beams, torque tubes, torque boxes, 
and concentrated springs (See Figure 1). The procedure followed for a typi- 
cal stress analysis was as follows (See Figure 2): 


a) Compute the stiffness matrices for the individual structural elements. 
Store on magnetic tape. 

b) Assemble the stiffness matrix for the complete structure by bringing 
together the coefficients for all structural elements and adding corres- 
ponding coefficients. Individual coordinate identification numbers were 
indispensible for keeping the coefficients in order. 

c) Invert the stiffness matrix to obtain the flexibility matrix for the struc- 
ture. 

d) Compute the deflections by determining the matrix product of the ap- 
plied load vector and the flexibility matrix. 

e) Compute the internal forces on the individual structural elements by 
computing the product of the appropriate deflections with the stiffness 
matrices for the individual structural elements. 


The largest problems, of order 106, were found to require a total of four 
hours to assemble the stiffness matrix and eight hours to invert in order to 
obtain the flexibility matrix. Double precision arithmetic was used in this 
operation. Each stress calculation for a given load condition was obtained in 
an additional one and a half hours of computing time. Smaller problems were 
solved in considerably less time. 

The analysis of beam elements was performed considering bending deform - 
ation only. A linear variation in moment of inertia from one end of a beam to 
the other was assumed. 

In addition to the stress analyses, normal mode analyses were successfully 
performed. In such an analysis the natural frequencies and normal mode 
shapes were computed by an iteration technique. 


IBM 704 


The IBM 704 is much faster and has a larger high speed storage capacity 
than the IBM 701, therefore, a substantial savings in time can be realized by 
its use. Experience indicates that a matrix representing 100 independent 
coordinates can be assembled in one and inverted in three hours. 

At Convair, Fort Worth, the IBM 704 program incorporated the following 
additions and refinements to the aforementioned IBM 701 program: 


a) Shear deformation is included in the analysis of beams. 

b) Effect of tapering flanges is computed. 

c) Rectangular structural plates may be analyzed using a lattice analogy 
representation. 

d) Beam moments of inertia are automatically computed if desired. 

e) Length of beams and their directions are determined by the computer 
from basic geometric input. 
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Remington Rand 1103 


Convair, San Diego, California, has coded the stiffness matrix structural 
analysis on the Remington Rand 1103. This computer is comparable to the 
IBM 701. The sizes of the problems worked by San Diego were in general 
smaller than those worked by Fort Worth. However the speed with which a 
given problem could be worked was greater than for the IBM 701. 


Others 


English aircraft literature contains many articles regarding the application 
of the stiffness matrix structural analysis, one of which has been cited.(5) It 
can reasonably be concluded, therefore, that this technique has been coded on 
English high speed digital computers. 


Application to Civil Engineering Structures 


Indeterminate Trusses 


One of the most common Civil Engineering structures is the bridge truss. 
Such a structure could easily be analyzed by the stiffness matrix technique. 
However, the analysis of a determinate truss problem is quite simple and 
would hardly be worthwhile to analyze on a high speed digital computer. This 
is not always the case in the design of an indeterminate truss. 

Consider the truss pictured in Figure 3. The distortion and loads on such 
a structure can be described by two orthogonal displacements at each joint. 
Hence the analysis of this structure would require a matrix of order 17, for 
its solution. If a secondary stress analysis is required this may be obtained 
Simultaneously with the primary analysis. However, a rotation dispiacement 
at each joint must be represented requiring a matrix of order 27 for its 
solution as shown in Figure 4. The stress analysis of such a structure would 
give the forces in all members and the deflections of all joints of the struc- 
ture. For a structure of this type the order of the matrix is generally much 
greater than the degree to which the structure in indeterminate. 


Building Frames 


Consider the simple building bent illustrated in Figure 5. The distortion 
of this structure is completely defined by the five coordinates illustrated: 
one rotation at each joint and one displacement at the roof level. 

A more complicated broken story bent is illustrated in Figure 6. Even in 
this case, however, only twelve coordinates are necessary. 

It is evident from these examples that the analysis of building frames 
could be accomplished to advantage using the stiffness matrix technique of 
analysis. If three dimensional frames are to be analyzed the torsion of the A 
building could be included, as well as axial deformation of the columns. 

The computed answer in this type of problem gives the rotation and trans- 
lation of each end of the columns and beams, and the moments in the mem- 
bers. 


Miscellaneous Structures 


Consider the arch rib illustrated in Figure 7. The arch structure may be 
represented by a series of straight beam segments between the load points. 
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The distortion and loading of the structure may then be represented by the 
system of coordinates illustrated in Figure 8. Twenty one coordinates are 
required for a stiffness matrix analysis. Moments, shears, and axial thrusts 
may be computed at each end of each segment. 

The stiffness matrix technique may be used to advantage in the analysis 
of arch dams, suspension bridges and towers. The analysis of shell type 
structures by the stiffness matrix method should be particularly advantageous. 
A stiffness matrix analysis on high speed computers would make it practical 
to analyze shells subjected to complex asymmetrical air load distributions. 


Natural Mode Analysis 


Most Civil Engineering structures are not analyzed to determine their 
response to dynamic forces. There are possibly two exceptions. First, the 
structure which is analyzed to determine the effect of earthquakes upon it 
and second, the structure which is designed to resist the effect of a bomb 
blast. 

In either case, the calculation of the dynamic response of the structure to 
a transient load is easily performed if one first obtains the natural frequen- 
cies and mode shapes for the structure. The mode shape analysis is easily 
carried out using the flexibility matrix obtained for the deflection calculation. 


An iterative technique which converges on the lowest mode has given the best 
results. 


Example 


The building frame illustrated in Figure 9 has been analyzed on the IBM 
704 for the load condition shown. The input geometry, stiffness, and load data 
was placed on 5 punched data cards. The load and deflection coordinates 
used for the stiffness matrix are illustrated in Figure 10. 

The stiffness matrix illustrated in Figure ila and the flexibility matrix 
illustrated in Figure 11b were computed, assembled and printed in about four 
minutes on the IBM 704. The deflections in Figure 12a and the stresses in 
Figure 12b obtained for the given load condition required an additional two 
minutes of IBM computing time. A comparison of these results with the an- 


swers given by Wilbur and Norris(7) to the same problem show good agree- 
ment. 


RECOMMENDATIONS AND CONC LUSIONS 


A major drawback to the application of high speed digital computers is the 
initial difficulty of getting a given type of problem set up and coded. This 
operation is expensive because it requires trained coding personnel and 
generally an expert in the field in which the problem will be used. 

If the problem is large it will require a considerable effort for checkout to 
the point where the program is operational. Therefore, the use of a digital 
computer can only be justified if several problems are to be solved, or if the 
problem can only be solved using the computer. 

Several alternatives are available to surmount this difficulty. The first is 
to join with other interested parties to develop the required programs on a 
joint basis. On the other hand if one is located near a higher educational 
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CALCULATED RESULTS 


.0383 rad. 0125 rad. 


5.120 in. .0303 rad. 


.0304 rad. .0399 rad. 


.00448 rad 1.760 in. 


1.889 in. 


(a.) DEFLECTIONS 


Moment (l.e.) Shear Moment (r.e.) 
In.~ Lb. Lb. In.~ Lb. 


(b.) STRESSES 


FIG. 12 


a 
AB .4118 10’ 33,074 .1835 107 
BG -.1835 « 10’ 27,603 ~—,5148 « 10’ 
GD 5148 10’ 72,397 .3540 « 10’ 
DE .5224 « 10° 22,176 .4800 « 10’ 
EF .4800 « 10’ 100,000 0 
cD 0 16,126 .4062 « 10’ 
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center with adequate computing facilities, the required coding may be eco- 
nomically obtained through monetary support of a graduate thesis topic. 
Finally, the required program may be developed by modification of an exist- 
ing program obtained through the IBM SHARE program. 

In conclusion, structural analysis problems of all types may be solved on 
high speed computers using the stiffness matrix technique. This procedure 
appears particularly adaptable to the use of high speed digital computers. 
The application of the technique to the solution of Civil Engineering problems 


presents no greater difficulty than those already surmounted for the analysis 
of aircraft structures. 


BIBLIOGRAPHY 


G. Kron, Tensorial Analysis and Equivalent Circuits of Elastic Structures, 
Journ. Franklin Inst., Vol. 238, No. 6 (1944). 


H. U. Schuerch, Vibration Mode Analysis For Delta Wing Structure, Con- 
vair, San Diego, Memo DG-G-100, 12-18-51. 


S. Levy, Structural Analysis and Influence Coefficients For Delta Wings, 
Journal of Aeronautical Sciences, Vol. 20, No. 7, pp. 449, July 1953. 


G. Best and M. P. Keating, A Stiffness Matrix Method of Delta Wing Stress 
Analysis, Convair, Ft. Worth, SRG-16, 9-8-53. 


D. Williams, Recent Developments In The Structural Approach To Aero- 
elastic Problems, Journal of the Royal Aeronautical Society, June 1954. 


Bisplinghoff, Ashley and Halfmann, Aeroelasticity, Addison-Wesley Pub- 


lishing Company, Inc., Cambridge 42, Mass., 1955. 


- Wilbur, J. B. and Norris, C. H., Elementary Structural Analysis, McGraw - 
Hill Book Company, Inc., New York, 1948. 


= 


Paper 1815 October, 1958 


Journal of the 


STRUCTURAL DIVISION 


Proceedings of the American Society of Civil Engineers 


NUMERICAL SOLUTIONS FOR INTERCONNECTED BRIDGE GIRDERS 


Henry Malter,! A. M. ASCE 
(Proc. Paper 1815) 


SYNOPSIS 


Two numerical methods for the solution of the problem of the interconnect- 
ed bridge girder are considered for the particular case involving infinite 
torsional rigidity of the main longitudinal members. Formulas are developed 
for both deflection and moment for a two girder and four girder structure, 
and results from both cases are compared with those obtained by using rigor- 
ous formulas. Tables are included to compare solutions for various values of 
relative stiffness. The possibility of replacing the entire cross framing by a 
single center diaphragm for analysis purposes is clearly indicated. 


INTRODUCTION 


The interest shown by the engineering profession in the problem of the 
interconnected bridge girder is evident by the increasing frequency of articles 
devoted to the subject. This interest is engendered both by the theoretical 
aspects of the problem, as well as by the practicality of the results that may 
be employed in the expanding highway program both in the United States as 
well as throughout the world. Various theoretical solutions to the general 
problem have been proposed in recent years. These include moment distri- 
bution, (1) anisotropic plates , (2) and a method employing a Fourier Series 
loading. (3) In addition, many experimental and test results have been record- 
ed and appear in the literature. A comprehensive survey of the various as- 
pects of the problem has been included in an article by Professors Siess and 
Veletsos, (4) whereas a complete bibliography of recent articles appears ina 
discussion article by Professor Makowski. (5) 

The general problem consists of a series of main longitudinal girders 
(simply supported or continuous) interconnected by either a slab, or a series 
Note: Discussion open until March 1, 1959. To extend the closing date one month, a 

written request must be filed with the Executive Secretary, ASCE. Paper 1815 is 


part of the copyrighted Journal of the Structural Engineering Division, Proceedings 
of the American Society of Civil Engineers, Vol. 84, No. ST 6, October, 1958. 


1. Prof., Hd. Dept. of Civ. Eng., Robert Coll., Istanbul, Turkey. 
1815-1 


ee ST 6 
4 


1815-2 ST 6 October, 1958 


of diaphragms or both. A load placed upon any one of the girders will result 
in deflections of the adjacent beams of the structure as well as of the loaded 
one. Part of the load is thus transmitted through the diaphragms or slab to 
the other beams. 

The two limiting conditions for the general solution of this problem are 
a) negligible torsional rigidity of the longitudinal beams, and b) infinite 
torsional rigidity of the longitudinal beams. The first case would include the 
usual I-Beam type of bridge with concrete deck slab, and is by far the more 
important of the two from a practical point of view at the present time. The 
second type would be encountered in a bridge made up of heavy concrete 
girders that are able to resist rotation, 

This article will concern itself with but a limited aspect of the problem, 
namely, numerical methods for the solution of case (b), that of infinite torsion- 
al rigidity. The question of the applicability of any particular solution to a 
specific problem is more a matter of interpretation of the behavior of the 
structure under consideration, and becomes rather indefinite. Any practical 
conclusions must be governed to a great extent by experimental results and 
observations on full scale bridges, and no effort to consider the above will be 
undertaken in this article. 

A solution for the particular problem being considered has been obtained 
by Pippard(6) and by Hetenyi, (7) both employing similar methods, and basing 
their derivations on a rigorous solution of the differential equations pertaining 
to the particular problem. Pippard has tabulated the solutions for 2, 3 and 4 
girder structures, and the results to be obtained by the approximate methods 
will be compared with those obtained from a substitution in the Pippard 
formulas. 

Although the numerical methods to be employed will emphasize only one 
particular problem, that of a concentrated load at the center of one of the 
longitudinals, it should be realized that this is being done merely for ease of 
comparison with the results obtained by the more rigorous method of Pippard. 
Also, this is usually the most critical case, and is therefore of greater im- 
portance. The method to be employed is just as applicable to distributed load- 
ings or to concentrated loads placed not at the center. The ease of numerical 
solutions for particular cases will be obvious, although general formulas may 
also be derived, as will be done for several cases. 

The various aspects of the problem can be understood by considering a 
two-girder bridge in which the connecting members are either a series of dia- 
phragms or a concrete deck slab. The stiffness of the main members will be 
designated as E I, and the stiffness of the cross members per unit length of 
the structure as El,. This latter term may be obtained by assuming any 
cross member as being distributed uniformly over the cross member spacing, 
and adding to this the stiffness of the slab per unit length. (All terms will be 
defined where first encountered, and repeated in the Notation at the end of the 
article.) A load applied to any one girder will result in deflections along the 
length of the structure that will be different for each beam. 

To enable ease of comparison between the various methods employed, only 
symmetrical structures will be considered. However, after grasping the tech- 
nique employed, it will be obvious that no restrictions are placed on the type 
of problem that may be handled. Unsymmetrical structures, as well as beams 
of varying moment of inertia can be handled in the manner to be outlined for 
the symmetrical system. It should be remembered that the case of infinite 
torsional rigidity of the main girders is the case being considered. 


NUMERICAL SOLUTIONS 
Rigorous Method 


The general problem could be solved directly. However, for cases involv- 
ing symmetry, as is here assumed (E I, =E Ip), it is often preferable to 
utilize the familiar symmetrical and anti-symmetrical load systems, and to 
combine the results obtained. Thus if two anti-symmetrical loads, each of 
1/2 P, are applied at the centers of beams 1 and 2 (Fig. 1), the resulting total 
deflection of 2y between the two beams will be observed, where y is the de- 
flection of each beam from its initial position. This total deflection of 2y will 
produce in the cross beam a shear of 


where w is the shear per unit length of the structure, and 1 is the spacing of 
the main members. This shear now appears as a distributed load on the main 
girders, and is acting upwards on girder 1 and downwards on girder 2. Inas- 
much as this distributed load is a function of deflection, it may be recognized 
as an elastic foundation reaction, and the equation for its solution will be 


ont 
1°EI 


where k is substituted for the term in the parenthesis, and the spring constant 
r may be defined by r = 2kKEI. A solution to this problem may be easily ob- 
tained by recognizing it as one of the special cases considered by Hetenyi— 
that of a finite beam both simply supported and resting on an elastic foun- 
dation, and with a concentrated load at the center. 

A consideration of the symmetrical case (Fig. 2) will indicate that both 


beams will deflect an equal amount, producing no diaphragm action. The de- 
flections will be no different from that of a simply supported beam with a 
concentrated load at the center. If the assumed loads are in this case also 
1/2 P, combining the two solutions will result in cancellation of the load on 
beam 2, whereas a total load of P will be found to occupy beam 1. 
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If one were concerned only with center deflections, or center moments, 
then advantage could be taken of the special case solutions offered by Hentenyi. 


For the anti-symmetrical case with the loads of 1/2 P, the center deflections 
will be 


4 
PA fsinn AL - sin AL 
Co (= AL + cos AL A= LE I 


using the plus sign for the downward deflection of beam 1, and the minus sign 
for the upward deflection of beam 2, 

For the symmetrical case, with loads of 1/2 P, the deflections of both 
beams will be the same, each one being equal to PL / 96 EI. 

Combining the two cases will produce a total center deflection of 


Er |2 Ir | coshATL cosAL 


where the plus sign is to be used for beam 1, now loaded with a concentrated 
load of P, and the minus sign is to be used for the unloaded beam 2. Since 
k = 2A4, this can be simplified to a final form of 


sinh 
Er 2 cosh AL ¢ cos AL (1) 
It will be recognized that the expression derived above is similar to the one 
obtained by Pippard, however with the trigonometric terms omitted. The 
influence of these latter terms will at all times be small when AL is large, 
and their omission will result in small errors, as is evident from the fact that 
the value of the terms in the parenthesis is 1.08, 1.05, 1.01 for AL values of 
3, 4, and 5 respectively. 

To obtain center moments, advantage is again taken of the special problems 
solutions offered by Hetenyi, but this time using the formula for center 


moment. In a manner similar to that outlined above for deflections, one will 
obtain the following final form for center moments 


+ mAbs sin At) 

(=: L cos AL (2) 
Again the trigonometric terms may be ignored for large values of AL. The 
value of the terms in the parenthesis is 1.12, 1.00, 0.982, and 1.00 for AL 
values of 3, 4, 5, and 6 respectively. 

These derived equations have been employed to obtain both center de- 
flections as well as center moments for various values of relative stiffness k, 
and the results are recorded in Table I for comparative purposes with other 
methods, In all cases the results are shown as percentages of total deflection 
PL3 / 48 EI, or of total moment PL / 4. 

A rigorous solution for a four girder bridge will not be attempted. How- 
ever, using the formulas for this case as provided by Pippard, the results for 
both deflections and moments for the various girders have been tabulated for 
different values of k, and are included in Table II where they may be compared 
with the results obtained by each of the approximate methods to be outlined. 


= 
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Center Diaphragm Method 


The center diaphragm method (so called for want of a better name), pre- 
supposes that the entire cross diaphragm stiffness (whether continuous in the 
form of a slab, or consisting of distinct cross beams), may be replaced by 
three concentrated diaphragms, one at the center, and one each at the ends. 
The center diaphragm will have a stiffness of hEI,, and the end diaphragms a 
stiffness of 1/2 hEI,,, where as before, El, is the cross stiffness per unit 
length of the structure, and h is the subdivision spacing, in this case 1/2 L. 

It should be evident that the ends of the beams do not deflect, and consequently 


the end diaphragms are of no factor in the analysis, and will hereafter be 
omitted. 


Two-Girder Structure 


For the case of anti-symmetrical loading, with loads of 1/2 P, the con- 
ditions will be as seen in Fig. 3. 


The shear in the cross girder at the center will be 
12 Ele, 2th 
13 
where C represents the deflection of the structure from its original position 


at point C. This shear will now appear as concentrated loads as shown in 
Fig. 3c, and the center deflection of beam 1 will be 


i 


12 EI 1 
EI Ce 2Ch | 
13 


whereas the deflection of beam 2 will be only the second term and using the 
plus sign. Substituting k for the term in the bracket, and 1/2 L for h, and 
finally for simplicity using the term @ for kL4/ 96 there will be obtained 


) 
1 ] ii t 1 
‘ A D 
2h = L 1 L 
(a) (b) (c) 
Fig. 3 
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For the symmetrical case the deflection of each beam will be PL3 / 96 EI, 
and in combining the two cases the final deflections will become 


p13 1+ 


To obtain center moments by this approximate method, it was observed 
that for both concentrated center loads and for triangular loadings, the de- 
flections at the quarter points of simply supported beams is approximately 70 
per cent of the center deflection. For a diaphragm spacing of h = 1/4 L the 
shear in the center diaphragm will now be 

12 Ele L (0) Cy ) 
for a beam loaded with a load P on beam 1. Substituting the value for 
(C1 - Cg) found above, the shear will become 


P 1 
(x35) 


The concentrations due to the diaphragms at the quarter points will be 70 per 
cent of the above, and.the conditions prevailing will be as shown in Fig. 4. 
For this loading the center moments will be 


L 
( ‘) 


and the final center moments will be 


Results for both deflections and moments are tabulated in Table I. 


Four-Girder Structure: Load on Beam 1 


To show the effectiveness of this approximate method, its use will be 
demonstrated by application to a four-girder structure (Fig. 5). To make the 
problem more general, it will be assumed that the inner girders are each of 
stiffness EI, whereas the two outer girders will be assumed as having a 


ve SC“ 
1 
‘ 
Fig. 
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stiffness of E nl. The unit stiffness of the cross members will be as before 
EI, and the value of the relative stiffness will be defined as the ratio 

12 El,./ 13 E nl, thus identifying it with the outer girder. The symbol k will 
again. be used for relative stiffness. (This definition is the same as that em- 
ployed in the two-girder case for which n = 1.) For this example two solutions 
will be required, the first with a load P on girder 1, and a second general so- 
lution with a load P on girder 2. 

The method to be employed is similar to that used for the two-girder 
structure, however, that for each of the two cases, symmetrical and anti- 
symmetrical, two sets of simultaneous equations will be required to find the 
two unknown deflections C and E. For the symmetrical loads the situation is 
as shown in Fig. 5b, and the deflections at C and E will be 

12 El, L (C- 5) 3 
13 E nI 


Using the same substitutions as before for k and ¢, and solving, will give for 
the symmetrical case 


1 
C 56 nl | 


X 


where X has been substituted for the term (1 + @ + n @) appearing in the 
denominator. 

For the anti-symmetrical case the beam positions will be as seen in Fig. 
doc. The loading condition for beam 1 is the same as that for the symmetrical 
case, and its equation will be the same as that found above. For beam 2 how- 


ever, the conditions will be different. The resulting equations will be (after 
substitutions) 


pI? (C-£) 


6 nj 


= ng(C-£E) ng (2é) 


A solution of these will give 


PI? 14 3ng 
96 n. E af) 43ng)-n 


96 EnI Y 
where Y has been substituted for the term (1 + ¢)(1 + 3nd) = n@2. 


Combining both cases will result in the solution for deflections caused by a 
load P on beam 1. 


13 2 48 ET 
| 
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(c) 


Fie. 5 


These deflections will be 


Long 
= 96 Ent X Y 


(5) 
Ee PI? n g + n g 
(6) 
To obtain moments by this method, it will again be assumed that the 
quarter point deflections are 70 per cent of the center deflections. The shear 
in the diaphragms, now assumed as existing both at the quarter points as well 


as at the center, will be a direct result of the difference in deflection of ad- 
jacent beams. These differences are 


PL? én g 
53) | Y | 


[3 


(#3 - C),) 56 B ni Y 


The concentrated loads produced at any point will be 
Pc> : [ 


where the terms in the bracket are as seen immediately above. 


For this type of loading, the final form of the equations for center moment 
will be 


(a) 4 
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1 +14 2ng 
= - (2 (7) 
170g f1 + 1. 
q (x )] (8) 


Results from the use of these equations as well as those for deflections 


have been tabulated in Table II where they can be compared with the results of 
the other methods. 


Four-Girder Structure: Load on Beam 2 


For beam 2 loaded, a similar method could be employed resulting in the 
following equations: 


mf 3 
C PL + 
G, = [+ | 


(9) 
pi? 1 

(10) 
vw PL 1.70 1 1 
[ (: ++) (11) 


Finite Difference Method 


In all elementary textbooks on Resistance of Materials there appears the 
following equation for the solution of beam deflection problems; 
at 
re EI 
As in many problems involving differential equations, it has been found that 
difference operators may be substituted for the differential operators. Using 


the “molecule notation for central difference operators”,(8) the fourth differ- 
ence operator will be 


= 2, “lh, 96, 1) 


It is assumed that one is familiar with the general procedure of finite differ- 
ence operation, and only its application to the present problem will be empha- 
sized. 

A beam resting on an elastic foundation and subjected to vertical forces, 
will be deflected in some manner, and will result in foundation pressures that 
will vary along the length of the beam. This foundation pressure that is a 
continuous function, may be replaced by a set of concentrations at the various 
subdivision points of the beam. The concentrations may be represented in 
several ways. First, one may assume a “stepped” variation (Fig. 6a), result- 
ing in subgrade reactions of rhB, rhC, etc., where r is the subgrade modulus, 


| 
| 
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h is the length of the subdivision, and the letters B, C, D etc. are the de- 
flections of the beam at the points designated. Using a “straight line” vari- 
ation (Fig. 6b) would result in a more accurate form of concentration 

(B + 4C + D)rh/6, but would require more labor. The parabolic formula sug- 
gested by Professor Newmark for his numerical integration process, (9) could 
be employed to give concentrations of the form (A + 10B + C)rh/12. The 


simplest method to employ is the “stepped” method, and it is the one that will 
be used. 


Correction Factor 


Before considering any bridge structure, a set of correction factors should 
be obtained for later use. These can be found through the solution by the finite 
difference method of a simply supported beam divided into 4 subdivisions, and 
with a load P at the center (Fig. 7). To obtain the unknown deflections at B 
and C one fictitious point will be required in the first equation. The following 
two simultaneous equations will be set up 


-B- 4 6B = lO 4 B 


3 

EI 

where P/h has been substituted for w at point C. Since the deflection at point 
A is zero, a solution of these equations can be easily obtained, and will pro- 
duce the following deflections 


-8B + 60 


P13 py 
B= TT 


These values are about 11 per cent too large, an error that is inherent in the 


use of the four division spacing. The correct values, as can be verified from 
any handbook should be 


+ C) 


(a) Stepped (b) Straight Line 


WW 
rhp rhc 
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In the solutions to be derived below, similar errors will occur, and they can 
be corrected for by substituting as follows 


and i for 
Lo 


126 
Two-Girder Structure 


In the two-girder structure (Fig. 8), it is assumed that the entire cross 
framing is concentrated in the evenly spaced diaphragms at the quarter points 
and at the center. When the structure is loaded with a set of anti-symmetric 


loads, the shear in the diaphragms, and the concentrations at the division 
points will be 


12 El, 
1 


s 


Fo = a 2kh C EI 


The above structure has two unknown deflections, B and C, and consequently, 
two simultaneous equations will be required. These will be 


) 
-88 4 6C C 


replacing kh4 = kL4/256 by R, and recognizing that the deflection at A is zero, 
these equations become 


6B = 


+ 6C 


i. for | | 
(2h 3B) s 2kh BEI 
+ EI 
ORB 
Ph? 
|’ 
A 3 
4 ty 
e 
+ 
bh = L 
Fig. 7 
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A solution of these two equations will result in the deflections 


1 
64EI | 2 ((3 4 
Fa f aap 


where Z has been substituted for the term (3 + R)2 - 8 appearing in the de- 
nominator. It is obvious that the terms 1/64, and 3/128 are in error and 
should be replaced by the correct values as outlined above. 

For the symmetrical case, loads of 1/2 P placed at the center of each 
girder will result in deflections of 


16 96 EI ond 
The combined values will be 


~ 
~ 
w 


+) ] (13) 


3 
PL 14+1/3+R 
EI | 2 +2/ Zz )] 


Having found the deflections, these values may be translated into loads at 
the quarter points and at the center. For the anti-symmetrical case these 
loads will simplify to 


FR 


34R 
PRIS (225. 
(tet) 


" 


| 
| 
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and the combined value will be 


= PR E 


For the symmetrical case the center moments are PL/8, and combining 
both cases will result in the final form of the equation for moment 


025 + R (14) 
i = Z ) 


The results for both center deflections, and for center moments for various 
values of relative stiffness k are recorded in Table I, where the close 
agreement with the other methods is evident. 


Four-Girder Structure: Load on Beam 1 


A solution of the four girder problem by the finite difference method does 
not lend itself easily to the setting up of simple formulas. However, individual 
cases of particular relative stiffness can be easily solved through the solution 
of four simultaneous equations each for the symmetrical and the anti- 
symmetrical case. 

The structure to be considered (Fig. 9) will, as before, have EI as the stiff- 
ness of the inner girders, and E nl as the stiffness of the outer girders. Dia- 
phragms will be assumed at the quarter points as well as at the center, and 


the four unknown deflections are B,C, Dand E. For symmetrical loading the 
equations required will be 


6B -R (B 


-83 + 6C WR (C 


6D = 


-8D + 6E 


- >) 
(C E) | 
B Cc 
| | | Eni 
Fige 2 
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For the anti-symmetrical case the four equations will be 
68-hC -R D) 


pn? 


2E ni 


-8B + 6C 


“R(C- 5) + 


6D nR (B D) nt (21)) 


-5D + 65 e mk (C-£&) nk (22) 


Solutions of the above equations will give the deflections, and once the 
latter have been found, the moments are easily arrived at. Solutions of the 

above for both deflection and moment have been tabulated for various values 
of k in Table II, and may be compared with the results of the other methods. 


Tables 


To indicate the accuracy that can be obtained by the approximate methods, 
Tables I and II have been prepared to show comparative results for both 
moments and deflections. Table I shows results for two-girder structures of 
length 600 inches and 1000 inches. For the two-girder structure both longi- 
tudinals have the same stiffness. 

Table II has been prepared to show the distribution of moment and de- 
flection on a four-girder structure with a load on beam 1 only. The ratio of 
outer to inner girder stiffness is 1/2. The moments appearing on each beam 
are shown as percentages of the maximum value PL/4. For deflections the 
values shown are percentages of PL3/48 E nl for the outer girders, and of 
PL3/48 EI for the inner girders. It should be noted that in Table II the use of 
the Pippard formulas for k values as low as 0.10 x 10-8 is not recommended. 
For this low value AL = 2.64 and a correction factor is required. 

To indicate the general trend in results obtainable, Tables III and IV have 
been prepared, using the Pippard formulas, although it should be evident that 
either of the approximate methods might have been used. Both Tables III and 
IV apply to four-girder structures, both 600 and 1000 inches in length, and 
having n values of both 1/2 and 1. For each structure two sets of results have 
been obtained, first with the load placed on beam 1, and then with the load 
placed on beam 2. 

Table III has been prepared to show the moment percentage appearing on 
each beam for the various cases considered. 

If each girder of the structure is loaded with a maximum load at the center, 
then the value of the moment in each girder will be 100 per cent of PL/4. 
However, normal spacing often precludes the possibility of loading all four 
girders at the same time. Assuming that only three of the girders are fully 
loaded, the moments were combined for each case, and the total moments are 
shown in Table IV. In the preparation of this chart for the case of n = 1/2 it 
was assumed that the load appearing on each girder was the maximum design 
load, and the inner girders were given twice the load of the outer girder. Two 
sets of combined maximum values were obtained, one for beam 1, and a second 
set for beam 2. It is of interest that the maximum obtainable moment is be- 
tween 80 and 90 per cent of that on a fully loaded structure, or of the 
assumption that each beam carries its own load. 
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CENTER MOMENTS 


L 2000 inches 


Fipp Diaph, Diaph, 
ard Meth, Meth, 


(center moments as percent of PI/)) 


CENTER DEFLECTIONS 
51.8 
52.2 


$2.3 


53.9 


6h. 


(center deflections as percent of PL? /48ET) 


TABLE I ;: Comparative Results for Two-Girder Structure 


ASCE 
1.40 59.1 58.6 59.4 55.5 576 56.4 
1.10 59.7 8.38 59.9 55.8 57.6 6.6 
280 60.5 59.) 60,8 56.3 57.6 56.7 
61.8 60.4 62,2 54.0 5762 
220 65.5 64.2 66.5 58.9 58.5 59.1 
10 70.2 68,9 70.9 60.5 59.7 61.2 
1.10 50.4 50.2 50.3 
200 | 52.2 52.8 50.6 50.3 50.5 
250 53.5 Shel 50.8 50.5 50.8 
220 57.8 58.2 51.7 51.2 51.7 
10 63.5 63.6 52.8 5263 52.8 


1815-16 October, 1958 


CENTER MOMENTS 
Finite Diff, 
29.3 29.4 25.4 10.4 
2961 29.5 23.0 
26.7 29.7 23.5 21.6 9.7 
30.1 21.1 19.6 68,8 
15.2 5 2 
20.3 967 9.4 2.2 


CENTER DEFLECTIONS 


3420 

34.1 28.9 
2767 
34.9 25.3 
36.3 17.3 
38.4 6.6 


PL?/48 E nI for beams 1,  ) 


P13/48 & I for beams 2, 3 
Load on Beam 1 lL © 600 inches ; J——I—J——-1 


TABLE II; Comparative Results for Four-Girder Structure 


( center moment given as percent of PI/l ) 

1.0 | 22.2 21.4 34.5 2969 22.4 34.2 29.2 14.0 
1.10 23.3 313.7 22.6 3h.7 28.9 13.6 23.6 34.3 28.5 13.4 
25.1 12.9 2.6 35.1 2767 12.6 25.4 34.5 2723 12.6 
50 | 28.6 11.2] 28.4 35.6 25.0 10.8] 29.3 35.2 25.0 11.0 
220 ho.5 5.9 40.0 35.6 17.8 6.5 40.2 34.9 1729 6.8 
010 | 56.5 - 152.0 33.2 3.215163 32.4 12.2 he2 


NUMERICAL SOLUTIONS 


600 inches 1000 inches 


34.8 
Belt 
36.3 
37.6 
4) 


MM PO 


ONO 


Moments shown as percent of PL/l ; load placed on either 
beam 1 or beam 2 as indicated 


TABLE III ¢ 


Four Beam Structure - bercentace Moment Distribution 


ASCE 1815-17 
k 8 P P ‘ 
My Mo M3 M), My) Mo M3 M), 
1.40 1.3 22.0 18.8 17.6 20.7 
1.10 42.3 21.8 18.5 17.3 21.3 20 6h: 
50 43.8 21.6 18.0 16.7 20,6 20.0 
46.1 21.1 17.1 15.8 20.2 19.4 
«20 51.5 20.1 15.1 13.2 19.0 17.9 
P 
Mo M 3 M), Mo M 3 
1.40 22.0 38.3 20.8 16,8 23.2 33.0 22.5 21.3 
1,10 21.8 39.0 20.6 18.5 23.1 33.4 
280 21.6 40.) 20.4 18.0 22.9 3L.0 22.2 20.8 
250 21.1 42,2 19.7 27.2 22.7 3543 21.8 20.2 
220 20.1 46,1; 18.4 15.1 22.0 3709 19.0 
a) is, = 
P 
Ml Mo K3 My, My No My, 
1.40 348 29.3 2he3 11.6 2765 30.9 28.0 13.6 
1.10 35.9 29.1 2367 11.3 28.2 30.8 27.6 13.4 
20 37.6 28.7 22.9 10.8 29.1 30.5 271 13.2 
50 hO.] 28.1 10,1 30.7 30.2 26.3 12,7 
F P 
My Mo M3 M), My Mo M3 
1.10 14.7 47.0 26.2 12.3 29,1 1.0 
1.10 14.4 7.8 25.6 11.6 15. 42,0 28.8 13.6 
280 143.3 12.5 15.2 42.8 26.4 13.6 
220 10.8 15.1 43.9 27.8 13.2 
220 13.4 55.6 9.3 11;.7 6.6 26.) 12.3 
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10" | 600 1000 600 1000 Goo 1000 
inches 


Moments shown as percent of PI/k for beam indicated with loads 
on beams 1, 2, and 3 of structure. 


TABLE IV: Four Beam Structure - Moments for Combined loads 


The value of the relative stiffness factor k in each table is in units of 10-8 
inches 4. 

The Pippard formulas employed for the various structures considered are 
rather lengthy and will not be reproduced here. They may be found on pages 
107 and 108 of reference 6. 
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Notation 


- span of structure 


- length of subdivision 


- distance between main girders 


M - center moment of beam designated 

P - concentrated load applied at center of span 

r - subgrade modulus 2kEI 

El, - stiffness of girder 1; Elg = stiffness of girder 2 

n - ratio of El, /Elg in four girder structure 

El, - stiffness of cross framing per unit length of structure 

B,C, D, E - deflection of structure from original position at point indicated 
k = 12 El,/ 13E nl 

= kL4/ 96 

R = kL4/256 

xX { 
Y = (1+ + 3 ng) - 

Zz =(3+R)2-8 
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SYNOPSIS 


The present study concerns the analysis of rigid frames, whose joints are 


not just points but small undeformed areas, by the method of successive ap- 
proximations. 


INTRODUCTION 


The stresses in the members of a framework are usually calculated with 
the assumption that the connections between members are points to which they 
meet and are connected together by rigid or pinned joints. 

In reality however joints are small areas to which the various members of 
a framework are connected, 

Fig. (la) shows a typical rigid joint k, where k is not in reality a point but 
the shaded area aa' - bb'. We can consider this area as suffering no difor- 
mation because the moments of inertia in the directions 1-1 and 2-2 are very 
large compared with those of the sections of members 1-1 and 2-2. In point 
of fact, in the case of reinforced concrete structures, after calculating the 
bending moments of the various members meeting at joint k, we work out the 
bending moments at the sections of the faces aa', bb', ab, a'b' and from these 
values we determine the required dimensions of the sections and reinforcement 
of the members. 

Thus we accept that the area of joint k suffers no deformation. 

In Fig. 1b the joint k is the wall of reservoir. In this case, due to the di- 


mensions of the wall, joint k can safely be considered as an area suffering no 
deformation. 


Note: Discussion open until March 1, 1959. To extend the closing date one month, a 

written request must be filed with the Executive Secretary, ASCE. Paper 1816 is 
part of the copyrighted Journal of the Structural Division, Proceedings of the 
American Society of Civil Engineers, Vol. 84, ST 6, October, 1958. 


1. Dr. Civ. Eng. Asst., National Technical Univ., Athens, Greece. 
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We can now proceed with the stress analysis of rigid frames with straight 
members by the successive approximation method making the assumption that 
joints are small areas suffering no deformation i.e. absolutely rigid bodies. 


Elastic Movement of Frame 


As it is known, the elastic deflection of a rigid frame under load is charac- 
terized by the rotation and displacement of its joints. By knowing these and 
the external forces acting on the frame we can proceed with its stress analy- 
sis. Also deflections due to direct and shear stresses are very small com- 
pared with those due to bending and can be neglected. Neglecting therefore 
deflections due to direct and shear stresses the displacements of the various 
joints are not independent of each other and by making use of the “virtual dis- 
placement” we can express the displacement of all joints as functions of the 
displacement of a certain number of them. 

If a rigid frame is converted into a pin-jointed frame by replacing rigid 
joints by pinned ones the resulting formation will as a rule be capable of 
suffering geometrical displacement. We can make this resulting formation a 
just-stiff frame by the addition of a certain number of pin jointed bars (roller 
restraints) Fig. 2. The number of bars required for this purpose gives us the 
degree of the elastic freedom of the original rigid frame. 


Sign Convention(2) 


Bending moments at sections by the joints are considered positive when 
clockwise. 

Shearing forces are considered positive when the shearing couple tends to 
turn an infinite portion of the beam clockwise. 


Normal forces are considered positive when they produce tension to the 
members. 


Stress Analysis of Rigid Frames by the Method of Successive Approximations 


Stress analysis of rigid frames of the classical case (joints-points) is 
carried out either in two (1) or one stages (2). Below are given the new 
factors for both cases. 


Stress Analysis in Two Stages. Ist Stage 


Here also we assume that joints are not deflected and that the bending 
moments shear and normal forces acting on the joints are those of a member 
built-in at both ends. We determine below the new distribution and carry-over 
factors when the joints are not allowed to deflect. 


Reactions on Joints When Fixed 


Assume all joints fixed (rotation not allowed). In this case all members are 
built-in at either end and exert to the joints bending moment, shearing and 
normal forces M,,, , Ne 

Due to these systems of forces the joint will rotate around the instantaneous 
center of rotation by an angle 0, (Fig. 3) _ 


By replacing these systems of forces /,, , 0. ; N,, by the equivalent ones 
we shall have 


ST 6 October, 1958 
Mai My Qu Oni Not (1) 


It is obvious that joint a will suffer the same elastic rotation 0, by either 
of the above two systems of forces. 


Therefore is this case the bending moment which acts on joint a is: 


where Ma represents any existing external couple on the joint. 

The fixing moment of joint a is equal and opposite to M, acting on joint a 
as calculated by Eq. (2). 

The forces on the joints by which we calculate the fixing forces f°. are: 


R=3,0,,+2; (3) 


where Pa: any existing external force on joint a. 

In Eq. (3) the summation is geometrical. 

The shearing forces 0; and 0,, are calculated by considering 8M, ,™,,,/%, 
and external load of span 1 and also 5, /,,, M,,and same external load of 
span (; . In other words the shearing forces of the built- in span /, are taken 
as being the same as those of an imaginary built-in span {; with the same ex- 
ternal load, and end fixing moments M/, and M,, (Fig. 3). 


Distribution Factors 


Joint a will, under the action of the B.M, Ma, suffer a rotation say ©, 
around the instantaneous centre of rotation a. The section «, of any member 
i will turn through and angle ©, due to the rotation of joint a and will therefore 
be displaced at right angles to its axis by a distance 4,xa@,; (Fig. 4). The 
section 6, will suffer no rotation or displacement because joint 6 remains 
fixed. 


Bending moments are given by the well known equations: 


a 


ds 
t 


where: «,; is the bending moment which should be applied on section o, to 
obtain rotation of it equal to unity i.e. the factor of stiffness to rotation of 
member i. ,,, is the carry-over factor of member i for the rotation of the 
section only and: u,,,4,; are the B.Ms at the sections a, and bj of member i 


for a ration Si = i.e. “factor of stiffness to displacement” of member i. (2) 
Shearing forces are given by the equations: 
=O = [Oxi (14444) + Up: ) Soi) a (5) 


Joint a is in equilibrium under B.Ms, M,,,M, and siaeaiian forces Q,, 
B.M. ™,, is given by the equation: 


May = Mi Qui (6) 


If in Eq. (6) we replace the values for Mai and Qai from Eqs. (4) and (5) we 
shall obtain 
‘ 


_ 
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a al d 
=i I+ (1444 ) Uai (1+ ) (8) 
The equilibrium of joint a is expressed by the equation: 
M,+2M,, =o (9) 
and by replacing the values for M,, from (7) we obtain: 


at 


From (7) and (10) we have: 


Mai=(-M) Kai (11) 


K 
bai (12) 
The factor K,; is the new “factor of distribution” for the case of rigid 
frames examined. 
The factor 6; is the new index of stiffness of the member and /6,, the 
new index of stiffness of joint a. 
If in Eq. (9) we replace the value for ™,, from (11) we have: 


=t! (12") 


Special Cases—Straight Members with Constant Moments of Inertia 


The factors bai for the case of straight members of constant moment of 
inertia have the following values. 


(a) Member built-in at both ends. 


! 
E d. q 
b,, +3 ( : (13) 
(b) Member built-in at one end—Pin joint at section 6, (Fig. 5) 
BE, d,,\? 
bi= [ 142 + ( (14) 


(c) Member built-in at one end—Pin joint at section ox, (Fig. 5) 
SEI, 
b (-Aai)’ (15) 


at” 


Carry-Over Factors 


We express bending moment ,, in terms of the calculated bending 
moment M,, 


+ 


1816-6 1958 
where: 
| 
where: 
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and by replacing the values for M,, and Q,, from (4) and (5) we obtain 
] 
Ms, = (Hi + Up; ) Uy + [> ai )+ + Up, ) a } On (17) 


and replacing the value for the angle 0, and working out we obtain: 


My Mab (18) 


a 


Axi 
Hap (9) 


Ap, + Up, ) 


ae 


which gives us the required carry-over factor of member i from joint a to 
joint b. 


Special Cases—Straight Members with Constant Moments of Inertia 
(a) Member built-in at both ends 


3 Max + +6 Ay 
t 


! 


al (da 
143 43 
l ai) 


(b) Member built-in at one end. Pin joint at section 6, 


(c) Member built-in at one end. Pin joints at section a, 
dy, 


Mab= 


Calculating of Bending Moments 


at 


The bending moments M., of the rigid frame are calculated by successive 
approximation as done in the classical case, 


For the calculation of bending moments M,, we work out at first the shear - 
ing forces Q,,, 


These are obtained from equation: 


Mar + Mp, 


(23) 


where (Q,; is the shearing force for the freely supported span [, 
Having obtained the value of shearing force (,, bending moments M,, are 
calculated from the equation: 


Qui (24) 


In the above equation B.M. ™,, andS.F. Q,, are entered with their ap- 


propriate signs. 


where 
(20) 
Ay, 
(21) 
l, 
(22) 
— 
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2nd Stage 


Having calculated the bending moments M,, of the Ist stage, we work out 
the shearing and normal forces Q,.,N,; . With these values and any other 
existing external forces on the joints we calculate the fixing forces F of the 
1st stage. The stress condition of the 2nd stage is the stress condition of the 
rigid frame acted upon by the forces on the joints Q,,, N,, and P or by an 
equivalent loading by forces Z equal and opposite to the fixing forces F of the . 
1st stage. 

We call “group of condisplaced joints” and “group of condisplaced 
members” all the joints or members of the frame which will be displaced to- 
gether should one only roller restraint be removed. (2) 

We number each roller restraint by I, Il, Il..... 

The solution of the 2nd stage in the classical case of frame is obtained, as 
it is known, by calculating the “partial stress condition” of the frame of a 
simple elastic freedom and obtaining a system of linear equations of the same 
number as the degree of elastic freedom of the frame by the help of which we 
obtain what are called the “convertion factors”. The required “stress con- 
dition” of the 2nd stage is obtained by the addition of the partial stress con- 
ditions of simple elastic, freedom having multiplied them by the respective 
“convertion factors”. (1) 

On calculating the above “partial stress conditions” by the method of suc- 
cessive approximation, the fixing moments are worked out with the supposition 
that the joints of the frame can be displaced without rotating. The relative 


displacements of the joints are obtained by the help of the virtual displacement 
diagrams. (2 


Virtual Displacement Diagrams 


As it is known from solid mechanics, any solid moving in a plane can have 
generally the possibility of two displacements and one rotation in the plane of 
the frame. 

The movements of the joints depend therefore on the number of members 
which connect the joint with the ground or other rigidly fixed neighbouring 
joints (Fig. 5a, b). 

Thus, joints a and b of Fig. 5a, by the addition of the member I have the 
possibility of rotation only. 

We call joints of double elastic freedom those which have the possibility of 
displacement and rotation. Joints a and b of Fig. 5a and joint a of Fig. 5b are 
of “double elastic freedom”. 

We call joints of “single elastic freedom” those joints which can only rotate 
around an instantaneous centre 0. 

In the rigid frame of Fig. 6, joints b and c are of a “single elastic freedom 
and joint a of “double elastic freedom”. By the addition of member I pinned 
at both ends, a can suffer only rotation whereas and c can suffer neither ro- 
tation nor displacement. 

For any movement 6., we can, by making use of the diagrams of virtual 
displacement, calculate the movement of any point of the joints of the frame. 


Fixing Moments of 2nd Stage 


In the classical case of rigid frames, i.e. when joints are taken as points, 
in order to calculate the fixing moments of the 2nd stage, we suppose that 


| 


ASCE SUCCESSIVE APPROXIMATIONS 1816-9 


joints do not rotate and by using the diagrams of virtual displacement we work 
out the relative displacement between two adjacent joints in the direction 
perpendicular to the member joining the joints under consideration, 

By making use of these relative displacements of the ends of the members, 
we calculate the fixing moments at all sections by the joints of the rigid 
frame. 

In the case considered and for the 2nd stage, we suppose that joints of 
“double elastic freedom” suffer no rotation during displacement, and joints of 
“single elastic freedom” rotate around the instantaneous center or rotation of 
the joint. 

Thus, when member I is removed (Fig. 6), joint a is displaced without ro- 
tation, and joints b and c rotate around their instantaneous centres of rotation 
Op ando. . 

Below are calculated the fixing moments of a member connecting joints of 
“single elastic freedom”. The equations given are general. Special cases of 
joints of double elastic freedom and also of one joint of single elastic freedom 
and one of double elastic freedom are also considered. 


(a) Joints of Single elastic freedom. 

In Fig. 7 we have joints a and b of single elastic freedom connected to- 
gether by member i and moments ,, and M,, at sections o, and 6, are re- 
quired as well as the moments at «° and ¢ distances a,, and a, from the 
sections ~, and 6, and for a displacement 6 of the point where the “imaginary” 
pin jointed member of joint b. 

Using the diagram of virtual displacement and in connection with the dis- 
placement placement 6, we calculate the rotations 0, , é, and 2, of joints b, a 
and member i respectively, by using the equations: 


4-019, (25) 
where: 
! ! e. 
Po €b0 Pa €b0 Can (26) 


In Eqs. (26) constants are either positive or negative depending on whether 
the rotation 6 around the instantaneous centres o,,o,0, are clockwise or 
counter-clockwise. 


Displacements 46,,,¢, ,4,, and d,=¢,, (Fig. 9) are found from the 
equations: 


Os, = Ca, 6, 6, Oni = Co; 6, Car 4, (27) 


where 


Special Cases 
(a) Members 0, 1. . . which are connected to joint b are parallel. 
In this case the constants 0 and c take the values: 


= 


Ch. San (26") 


Cy=tt Car=t Pala: Cr=*Pa ar (28") 
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b) The members connected to the joints a and b when taken as a total are 
parallel 


(26") 


Calculations of Moments Ma;, Mg, 


Sections oc, and 6, (Fig. 7), due to the rotation of the joints a and b by 0), 
and ©, will rotate themselves, in accordance with the rules of solid me- 
chanics, by the same angles ©, and 0, and in the same direction. 


The equations with the help of which we calculate the moments M,, and 
M,, are: 


5 
Mas = +H, ) + Uni i, 
re) 
Me, (Hi Xai 2, + Up, (29) 


In the above Eq. (29), the moments due to the rotations by ©, and 3, are 
entered with negative signs because a clockwise rotation of the sections by 
d,and 4, causes negative moments to the sections «, and 6, . 


Substituting the values for &,,4, and 3, from (25) and (27) in (29) we ob- 
tain: 


where 
‘ Xf 


x a 
Uni — Us, [ - Mi 


(31) 
The shearing forces are given by the equation: 


d, 


For straight members with constant moment of inertia, Eqs. (31) take the 
form: 


1) Member built-in at both ends. 


EJ, 
t t 
6EJ, 
=| (2 Pa ) (33) 


Us, = + 2Op) ci] 


2) Member built-in at one end—Pinned at section 


3EJ 


(34) 


= 
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3) Member built-in at one end—Pinned at section 
=0 


3EJ, 


Calculation of Moments Maiand M4, 
The moments M., and Ms, are given by the equations: 

Mat =Mai + Qari 

= Moi + (36) 


If in the above Eqs. (36) we substitute the values for M,,,/,, from Eqs. 
(30) and Q,,, Q,; from Eqs. (32), we obtain: 


(37) 


where: 


= US, (1+ + up, dae 
t 


Up, = Ua, + ug, (1+ tee) 
Comparing Eqs. (30) and (37) we draw the conclusion that Eqs. (37) are the 


general form of those of Eq. (30) because when a, = d;, = 0, by substituting 


these values in Eq. (37) we obtain Eq. (30), i.e. the moments M,,,M,, at 
sections and 6, . 


In accordance with the above, the moments M,,and™,, are calculated using 
Eqs. (37) where, for d,, = o, = 0, we have 


Ugi = “ai = Us, (39) 


Joints of “Double Elastic Freedom” 


Joints a and b of Fig. 8 are of “double elastic freedom” and are connected 
together by the member i which has an instantaneous centre of rotation 0, . 


The diagram of virtual displacement has been drawn for a given dis- 
placement 6,. 


Angle 0, can be expressed in terms of On (for infinitely small dis- 
placements) by the equation: 


(40) 


(41) 
where: 


(42) 


Joints a and b when displaced suffer no rotation. By making this suppo- 
sition, sections «,and 6, will take the positions «! , 6, (Fig. 8). The relative 


displacement of sections a, , 6, normal to the direction of the member is that 
calculated 6 . 


= tai 
Ms, = Up; 
Sb 
t 
From Fig. 8 we obtain the equations: 
Com gt |_| 
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The difference of the constant c, for joints of a single and double elastic 
freedom results from the fact that in the first case we take the span |, and in 
the second the span {. . 

The moments M,,,%,,M,,; and M,, are calculated by making use of the 
general Eqs. (37) and (30) in which the factorsu,, G,,, u%, , ui, take the 
values: 


For = = 0 we have and 38)) 


=C; Us Us; = G Ubi 

(43) 
Uy, = [Hae Up, (44) 


In the case of straight members of eeteihank moment of inertia, Eq. (43) and 
(44) take the forms: 


1) Member built-in at both ends 


(45") 
_ / 


= FA ad 


C; 
Joint A of Double Elastic Suintee~debs B of Single Elastic Freedom 
In this case (Fig. 9) we have: 


where 


Pp Cr: es, (48) 


The displacements are given by the equations: 
Oa = Ca 6, — Oar =C (49) 


ar 


where: 


(50) 

The Moments ™,,,™,, and M,, are calculated by using the general Eqs. (30) 
and (37) where the factors u®, _u%, and «, take the values: for @, = 0 we 


have: (Eqs. (31) and (38)) 


51 
Us = Up (- (51) 


ra=+ far 


(1+ )+ Up, Sai 
t 


J 
° 6E J, 
Uni =U 
(45) 
Up, = (142 ) 
2) Member built-in at on i 
¢ 
Up = Ubi | 
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For straight members with constant moment of inertia we have: 


1) Member built-in at both ends: 


6E/, 
(— 2p 


/ 
E 


2) Member built-in at one end—Pinned at », 


Up, 


3) Member built-in at one end—Pinned at o 


3E J, 
EX 1, +c.) 


Fixing Forces of 2nd Stage 


The fixing forces of the second stage are calculated in a more simple way 
by the principle of virtual work. 

We denote the diagrams of virtual displacement obtained by the removal of 
each of the imaginary pin jointed members (roller restraints) I, II, 
the general sy mbols (/7') (N”) respectively. 

Also by 6, , 0? c} ...., we denote neem 6 and 
constants c of the respective es. for the removal of I, 0. 

The displacements 5; | 6” ... andthe constants c’ , ci ... are calculat- 
ed with the help of the virtual displacement diagrams (7’), (1”), in conjunction 
with the geometry of the figure and the respective Os Gg ees 

Also we denote by ,... .... the rotations 0 and the 
constants o of the joints of single elastic freedom. 

If we now detach the members of the frame by taking sections close to the 
joints and apply to the joints forces Q), N,, Ff’. ... andthe moments 
M., for joints of single elastic freedom or ™! for joints of double elastic 
freedom, each joint will be in equilibrium under these system of forces. 

The equilibrium of the joints (taken as a whole) can be expressed by the 
principle of virtual work. 

If we apply the principle of virtual work, e.g. for the stress condition which 
has resulted by the removal of the imaginary member I and (virtual) defor- 
mation the geometrical one which is taken from the diagram of virtual dis- 
placement (/7'), we shall have:(2) 

E'S, +3 On Min — =0 (55) 


a*a 


where: 0; Qh, 


>,,: summation extending to all joints of single elastic freedom. 

2, } Summation extending to all members of the group of condisplaced 
members I. 

The sign of the third additive in (55) is minus because positive shearing 
force and positive displacement 6) result in negative work. 


| J (53") 
us, (54) 
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The work done by the forces £' , £’, is equal to zero because the dis- 
placements of the points of application of II, III. . on the joints are normal to 
the directions of the imaginary pinned membersI, I... 

The forces N’ do not appear in Eq. (55) because the work done by them, 
being equal and opposite, cancels each other. 

Also the moments M., at the joints of double elastic freedom taken at each 
joint, give work equal to zero. 

If in Eq. (54) we substitute the values for On and 6; as a function of 6, , 
we obtain 


Using Eq. (55) we calculate the fixing forces of the 2nd stage. 


Stress Analysis of Rigid Frames by the Method of Successive Approximations 
in One Stage(2) 


In the case of stress analysis of rigid frames by the method of successive 
approximations in one stage only, after the calculation of the fixed end 
moments and reactions, for joints suffering neither displacement nor rotation, 
we free to rotation successively the joints of double elastic freedom. Using 
the distribution and carry-over factors, the moments are calculated. By the 
help of the restraining factors(2) which will be calculated below, we work out 
the fixing forces f° ~E*. ..., Which are developed at the im- 
aginary pin-jointed members due to allowing the joints to rotate. 

After the successive freeing to rotation of the joints of the frame and the 
calculation of the fixing forces, we remove successively the imaginary pinned 
members applying at their points of connection equal and opposite forces to 
the calculated fixing forces. In this case we consider that joints of double 
elastic freedom are displaced without rotation and joints of single elastic 
freedom rotate around the instantaneous centre of rotation of the joint. With 
the help of the distribution factors due to displacement and shear carry-over 
which we shall work out below, we calculate the moments at the joints of the 
rigid frame and the fixing forces at the imaginary pinned members. The 
above work is repeated successively and thus we obtain the stress analysis 
of the rigid frame by successive approximations in one stage. (2) 


Restraining Factors 


In the stress analysis of (Fig. 10) all joints of double elastic freedom suffer 
no rotation except joint a on which acts the moment Ma. 


To generalize the equation, we denote moments with the general symbol 
where: 


(M2, ) = for joints of single elastic freedom 
Ms, )'= for joints of double elastic freedom 

Lengths of members (; are also denoted by the general symbol (ly where: 
(4 ) = (, for member joining joints of single elastic freedom 


for joints of double elastic freedom 


(4) =t+a,, for joint a of double elastic freedom and joint b of single 
elastic freedom. 


= 
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Moments (4, )'and (/{,)'are calculated from the equation (Fig. 16) 

(MS )'=(-Ma) 

(Mo Ma) Xai Hab (57) 


The shearing force is given by equation: 


\ 
The calculation of the fixing force F}* is given below. The remaining Ff 
, are calculated in a similar way as f° 
” If we express the equilibrium of the forces on ‘the joints M2 QQ? =N? and 
Fa;, Fa; ... (taken as a whole) by the principle of virtual work, taking as 
virtual deformation the ( 7’ ), we shall have the equation (Fig. 10) 


F° 6, cosd 2, Q d+ 2, Mg. 9, =0 (59) 


The moment Mm’ ,asa function of the factor of distribution due to rotation 
K,, and the carry-over factor ,,,., is given by the equation: 


=m ( M, ) Mar Kar (60) 


If now in Eq. (59) we replace 3; and 0) as functions of 6, and substitute the 
values for Q* and ™° from Eqs. (58) and (60) we shall have: 


(-M,) (61) 
where 
6a: Kail!* Mat) 2, Max (62) 


cos 0, i 


The summation extends to the joints r of single elastic freedom next to the 
joint under consideration (Fig. 10) 


For 0, = 0 we have: 


By the Eqs. (62) or (62') we obtain the restraining factors of the joints of 
double elastic freedom of the rigid frame. 


The dimensions of the “restraining factor” are 1/length. 


Distribution Factors Due to Displacement 


After the removal of the imaginary pinned member, e.g. I, we apply on the 
joint at the point of connections, force equal to -F , i.e. equal and opposite to 
the calculated fixing force at position I. 

where: 
F° is the fixing force at position I for the case of perfectly rigid joints 
(neither rotation nor displacement allowed) of the frame. 

Due to the force -Fy, the joints of the rigid frame will be displaced. We 

suppose that during this displacement, joints of double elastic freedom will 


not rotate and these of single elastic freedom will rotate around the 
instantaneous centre of rotation, 


| 
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Due to the displacement of the joints of the rigid frame, forces Q) , N 
and moments ( rt, )' will develop at the sections next to the joints. 
The moments (/™,, )' and ( ™,, )', are given by the general Eqs. (37) 
(Mai ) 
5 (64) 
(Mo) = Gb 
The shear forces are given by the general equation: 
Q, Qa as ( Ly (65) 
If we express the equilibrium of the forces acting on the joint of the rigid 
frame by the “principle of virtual work” taking as virtual deformation the 


elastic one given by the diagram of virtual displacement (7'), we shall have 
the equation: 


CF) 6, cos 8-2, Q; dj +2, 0' 5, (Mj, )=0 (66) 


r 


where 2, summation extending to all member i of group I. 


>,. summation extending to all joints r of single elastic freedom of 
group I. 
If in Eq. (66) we replace the values for QI,., and (MI,,)' from Eq. (64) and 
(65) and express 3; and 9! as functions of the displacement 5’ we shall ob- 
tain: 


(67) 
Combining (64) and (67) we have: 
(Max) = (-F) (68) 


Eqs. (69) and (69") give us the distribution factors due to displacement. 


Check 


The shearing forces as function of the moments (M, )' are given by the 
equations: 
Q' + Bou 
t ( ) (4) 
If we apply the principle of virtual works for the calculated, shearing 
forces, normal forces and moments of group I, due to -Fy and virtual defor- 


mation” the geometrical one given by the diagram of virtual displacement 
(n'), we obtain: 


cos 
where: ‘ 
In the case of angle 0, = 0 Eq. (69) becomes: 
| 
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re) 
‘ 4 (Ly t r r Bre 
t 


If now in the above equation we place the values for 3) and 3! in terms of 
the displacement 5 and the constants c! and «9! , we have the equation: 


1 
t Or 


For angle © = 0 the above relation becomes: 


1 


(4) ‘ x rk 
‘t 


Remarks 


The factors u), take the following forms depending on whether the joint is 
of single or double elastic freedom. 


1) For joints of single elastic freedom: 


(Mi, )'=M; (equation 31) 


ut at “cut 


2) For joints of double elastic freedom: (/! 


) ° 
t Up, 


(equation 38) 


Shear Carry~-Over Factors 


In the rigid frame, besides the pinned member I, exist also II, IIl,... 
Due to the removal of I and application to the joint where I is connected of 
force -Fy, fixing forces Fly, Fly, . . are developed on the remaining 
members Il, III, . . . of which we require the value in terms of -Fy. 

We give below the calculation of the fixing force Fly. 

If we express the equilibrium of the forces acting on the joint, by the 
principle of virtual works, taking as virtual deformation that given by the dia- 
gram of virtual displacement ( (7° ) we shall have the equation: 


If in Eq. (70) we express the displacements 5’ , & in terms of the constants 
> a e” and displacement d, » and solve it with regards Fy, we obtain: 


-2 


a 


r 


(71) 


cos rol 


The shearing forces ©’ and moments (™!,)' are given by the equations: 


1 


(73) 
(™,, )=( 


If in Eq. (71) we replace the values for Ql and (M,, )', from (72) and (73), 
we obtain: 


( F 


(74) 
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"= | (Feu * ~2r e; | (75) 
For the case of angle 0, = 0 Eq. (75) becomes: 


where: 


2, summation extending to members x connected to joint r of single elastic 
freedom. 


Layout of the Work of Successive Approximation 


The layout of the work of successive approximations is given again here in 
a table form (example 4th) as in the case where joints are taken as points. (2) 

The final total moments are obtained by adding up all values, in the 
squares, of the bending moments at the sections by the joints of the members 
of the rigid frame. 


CONCLUSION 


The present study constitutes the general form of the method of successive 
approximation for rigid frames whose joints are not points or infinitely small 
areas, an assumption which is made in all up to day classical methods of 
stress analysis of rigid frames, but areas of limited dimensions, which at the 
analysis of the rigid frame, are taken as solid areas suffering no deformation. 

We can say that the equations given, both for the analysis in one stage or 
two stages, constitute the generalized form of the method of successive ap- 
proximation. The equations known up to day for the analysis of rigid frames 
can be considered as a special case of the generalized ones. 

Besides the theoretical interest of the study, in certain cases either due 
to the interposition of partition or due to large dimensions of sections of 
certain members of the rigid frame, it is possible to have types of rigid 
frames for which the joints cannot be considered as being either points or 
infinitely small areas. The stress analysis of such types of frames by the 
method suggestion is in our opinion, simple and rapid. 

In order to make more clear the application of this method, we give below 
four examples. 


1st Example (Fig. 11) 


Moments 


2000x65* 
—7080 Kg™M 


2000x5.5* 


12 =SOSOKg™ -M. a 


~ 2 
2000 xS 
= 4160 KgM 
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Fig. 1 


Fig. 12 
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2000x5.5 ° 2000x658 
= 5500 KQ maa -~m6500 KG 


° 2000 x5 
50009 


Fixed-End Moments 


— 
505045500 x0.25=6425KgM 


7080+ 0.30x 6500 =+ 9030 kgm —7080-025x6S00—-8700 kg 


+4160+0.25X SOOO = +5410 


Distribution Factors 


d,-025 A,=035 
fell as 
d O2s = 0.045 =0087 dp, 0.30 =0054 Gee x 

55 le 400 t, 55 6 
d 0.35 d 0.25 ad 0.25 035 
1. 400 65 500 400 


0046 


Working out for b,, we neglect the quantity ( $a )2 as being very small: 


Joint a: 
4Ex275xI 2.27 0.643 
= (14 3x0045)—2 = 0-94: 
4E x] 1261 0.357 
Goo” (143 x0 087) £1 3531 1000 


4Ex275] 
(1#3 X0.054) —233 El 


4Ex2751 192 

(1+3x 0.046) =I92 El ‘ 6.27 0.306 
4Ex16!1 1+3x0.087) = 

( 627E! 


Joint c: 
4EX278! (143%0038) —169E! 


65 


4Ex275!1 
(14 3x0.05) m2 53 El 
126 El 


4eEx! 
(14+3x0.067) 


Carry-Over Factors: 


0.25x0.30 
"5.52 


Shearing Forces 
Joint b: 
233 
| 
189 
0,332 
568 
253 
568 0.446 
126 _ 0.222 
568 1000 
Hap 2 
Had 1261 
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030x025 25x030 
030x025 , ,_0.25XC 


oo 


035 
$.. 
126) 
03040 2S 0 25x00 30 
+6 — 
— 568 
1114 
O35 
1+ 
5S ) j om 4 
1261 
The successive approximations and the values of _, are given in Table I 
Shearing forces 
Q. = 5500+ 
al 
9468 — 7674 
Q,,—- 6500 705 
+ 3408 
»000 § 25 
Q 191-95 
bs 435 
Total moments: 


Mo 2254 —4320K0 25+ 1174 
Mp3 9468 — 6754 x 0.25= 7783 
M. +6921 —-5670X0 25 —550! kQI!) 


2nd Example (Fig. 12) 


Distribution Factors: 
Joint a: 


4Ex4] 


550 


=291E!1 
OBOE! 
500 


Joint b: 


7BE]I 


4Ex3lf 1 \? 
33 +3 -SBIEl 


0.30.) 
4.70 / 
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, 2.25+0 30 |, 025x030 
55 
Mba 2 1162 
Mp, 2 1138 
My me --9277 +0080 XK 30—=-6277 KOM 
M 7674+ 6246 —-6144 
le 
! 291 795 
arn 39) 785 
0.80 0.215 
4Ex4! ! 478 
0.501 
4E] hol 0087 
470° «(('+3 1166 
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October, 
Joint c: 


300 
= X 620 


4Ex3/ 240 386 
—=—=240E! K3= 620 oO 
OBOE! 080 _ 0.129 


=500  620EI 620 ~ 1000 


=—0485 


Moments M,; 


000x4? 
x 


=2660 Kg 


2000x5* _ 
a3 12 


Shearing forces: 


7 ° 
Q, = Qs, 2 6860 XQ 


2000x40 


— +4000 KG 
2 


Fixed-end moments: 


+6300 xg -13180 


‘ 
M,,— + 2660 +4000 x l= 6660 — 2660 


4160 
Carry-over factors: 


Vab=> L =0.773 


4 (+355) 


470 
> 


(1+ 3700) > 

1 3 


For moments M,, see Table II: 
Shearing forces: 


1925-12895 
—6880 + 8570 


1680-176! cago 
0,.— +4000 +4 11689=!76! sogo x 
ji 500 9 


Total moments: 
= 1925 = Kg 
al 


M, +11690~5980x10=-S700 KQ == —1761 x9 


g= + 2366 Kg Mys= 


43230) 
S00 
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0.372 |0322 


be 


16425 6700 |} +5410 


-2295 |-4130 |- 2400 


26 | 


+620 +1460 


200 | - 190 


500| 0 SOO 


+4160 


67 18 - 
+3) 
15 ~4 t2 


ASCE 
Table I 
= 
Kai | — |0.357 |0643 | 0 306 |0 332 |0 222 |0.446 | — | — | 
q ad ab ba | & cb | cf 
Hab | 
0.500 |O 580 0565 |0500/0558 
-4160 
148 | 
+15 + +12 
+364 | +730 
| | } | 
|-130 |-230 |- | - 106 | - 100 
| | 
+20| +41 | +74| +43 | 
-14] -13 ? | -7 
Table 
Hab 
|_| 
+6300]-13180 +6660 |- 2660 4160 
-675 |-1350 |-4950|-3830 
+2000 #4260 |+ 900 + S190 #2350 +450 
215 |-430 |-1570]-1210 ~1640 |-1870 |- 500 |-1480 | -740 - 250 
+ 550 #1170 |+ 250|+ 1430 |+ 645 +125 
- SO - 432 |- 334 - 274 |- 313 83 |- 249 12S -42 
+127 |+ 250 | + 53 |+30S +27 
- 28 27 |-100 77 - 59 - 26 -9 
— 
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ist Stage 


Joint b: 


Joint c: 


2nd stage: 


a) For 


Distribution factors 


Carry-over factors: + 


The moments of the Ist stage are given in Table III 
Restraining forces of 1st stage: (Fig. 16) 


Moments relative to b: 2H+!0375 -11.199-16.104x8+16x4 


Moments relative to 0a: 


Calculation of «2, and uj, 


ST 6 


Example No. 3 (Fig. 13) 


Va*+ 2? =447m 


Active moments 


_ 4Ex3I,_ 200EI, 


October, 1958 


22m 


l, 


10375-8.714 
16 


= 164 


=—16.104t 


H=3282¢ 


8.94 F =32.82x8+10375 F=3056¢ 
F=- 32.91 + 32.47¢ 


Virtual displacement diagrams. 


removal of imaginary roller restraint I: (Fig. 


c,-0895 @= $94 400 0.224 C,=—0.224x8 
1 1 

6, =-1.846, Cj=+184 6, = +1846, 


14) 


= 0.695 5, 


22=-1 84 


22-1065 
2 
Mes=—Mae= =6 tm 
a 
0.500 Kp3= 0500 
E I a 
_ 4ExSi,_ 
= 2.44 6.44 = 0.448 
4ExI, 100 
=100E], KXea™= Baa =O. 184. 
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General formulas: 


6EL=100im 


447.) 6El, 


447 


El, 
= (2x0. n2x49°+ 0695) tm 10tm 
(2x0.112 x 922-184) SEX Sb — 746 tm 
8.22 
(0.12 x - 184) -93.4tm 


6ESI, 
— 
u ue 1.64 112.0 tm 


9 
M,2=-74.6 s5>=-9.07 tm M,=-234 =-1135 tm 


The moments for removal of the imaginary roller restraint I are given in 
Table IV. 


Restraining forces—Shear forces 


I i 
+6660kgmM +5830KgM 7450 xgM 
Mp -8859 xgm Mp=-10928 xgm 


1 


Mig -5077 xg Mya-2536x9m Mie Mi Mn 

Q) = =2400 xg 7016 #102928 2.260 Kg 

Of = =2590- 1269 _ _9s50x9 


i 
= 2800+0.695 x 3490+2400x 184+ 2260x184 +0.112 (6660+ 7450 ~ 8859) 


F'= 


=-2400 x 2065 —2260 2055 -0.950x 1=410.53¢ 


29.8z15=666 tm 26 =583 im 
= 298-=745 tm tm 
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For removal of imaginary roller restraint II (Fig. 15) 


1816-29 


— d,=-20554 C=-2055 &=0 


« 
+2055 


6E[, = 1004 


6EX5lg 2055 
Me 822 E22 


50 xQgm 


=15200 


= 15200 xgm 


The moments for the removal of the imaginary pinned member II are given in 
Table V 


Moments: 


M 


ai 


M,,= 0 


zx Zz 
+14672 Mp +14144 


Shear forces 


14672+14144 14144-1149! 


= 3510 XQ 9 


“ 349441749 


=873 Kg 


250+6250 


Restraining forces: 


184X 3510 —-184x3.20+0.112x 14672=-10.54t 


3 510X2 05543120 X 2.05543 780433 125—2684Ft 


Convertion factors: 


+15.19 X, 


10.54 X, + 3056=0 


-10.53 X, + 2684 X —3247=0 


X,=-1.62 9570 


Total moments: 


eo 


-~6660 X1.62——-10780 KQ™M K9M 


M,;~ 


+8859 162 + 14672 X0.57+ 10375 = + 33085 


+10928 1624 14144XO.57-11199= +1456! KQM 


400M 
| . 
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-9070 5 5 + 13600 
-6100 

+200 + 200 
~ 90 


Table Il 
[owe [ows [ome [= [= | 
+10650 |-10650 |+10650 |-10 650 +6000 |~6000 4g 
+ 1040 |+2080 \. 860 #1710 + BSS + 430 
- 260 S20 | - S20 | - 260 | 
| + 58 #116 +48 + 96 +48 + 24 
Table V 
Ra [= [ease [ows [oom | 
~2500 |- 5000 2500 |-1250 
+22 +23 +12 
-6 -2 -4 
Table V 
| — | 2800 | | | | | | - | 
+15200 |+#15200 |-15200 |4+62S0 +6250 
+2000 |+ 4000 |+ 1650 |+3300 |+ 1650 |+ 825 
~500 | -1000 |- 1000 | — S00 
+112 | +224 | +92 | +104 +92 +46 
- 28 =-56 - 56 -28 
+13 | +S +10 +5 +3 
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14144 xO°* 411799 


7616 X162—11491 XO S7?-8714 = 27574 
+2539 x 1624+ 7997x0574 908 

> + 7124 +454 6564 xqm 
50O77X162+ 3494057 + 7806=+4 18006 Kor 


2536 16241747 xOS?7- 


K6 da? 


Example No, 4 


The rigid frame of example 3 is worked out again by the one stage method. 


Constants of diagrams of virtual displacement: 
Group I: 


Group II: 


Distribution factors due to rotation: 
Joint b: Joint c: 


Khe osoo 


! 


> 


Carry-over factors: 
Restraining factors: 


For the case of the rigid frame considered we have for all members: 


General formula: 


Group I: 


ASCE 
M 
c3™ 
Ma 
i 
Chm +l C0895 Chm 
e) Cj=-205S Ca = 
c= 
0448 
k 360 
I 
(/ i) ‘ai 
3 Gs 
= 
6,,—-2; > Xai i, 
164 o500+ [&.3..0500 (-o1n2) 229°.+0028 
8 « 
Group II: 
2055 
0448 x =-+ O rf 9 
8 ‘ 2 4 aa 
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2000 xgim 


COT 


10375 11199 


| 


+2004 400 —+ 


R% 
yi | 
4 4° 
Va V, 
800 + 8.00 
/ 
a/ | 
16t 
BO 
A® 
ty 
{0650 
800 + 800 
Fig 17 


Group I: 


ol 447 
Uso =(2x0112x— 


47 


47 6) Sh =—/790EI, 


+1 


4a? 


400 
(2x0.112x — + 


22 
us (2x0 


bem (0.112 


Up, 164 
or 
Ugo + 


4.48 


a o > 
l, 
un 
i > 
n2( 


382 
+O 3€ 
6 72 /822 


Check: 


a 


a4 


(- 0.112) (039040436 


Group II: 


on 


on, 
Us, = B22 


on a 
— 


Cc 1.790 41.565 


90+ 1 


822 


0.3904034) 0436+0 382 


-0.531)— +1004 


Uno Uy 


ar os a a 


0.895) = =—1790EI, 


400 


~184 448 El, 


822 


=—672El, = 


447 447 


565 O895 


~ 400 400 


+560 


184 
822 


6.72 184 


447 4 8 22 


Calculated «/, : In the case examined we have: i’ 


El, = + O275El, 
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Distribution factors due to displacement: 


Ug, Ug 


El,=+0.168 EI, 


EL= + 18 7 é I, 


+0367EI, 
ElL= 


, L79EL _ 4.48 El, 


O3IEI, 


—(+0.997+003) +1026El, 


1565 447 
+0341 = 
1028 ar 
448/822 
= 53) Q,.= 
ae 1028 tbe 


——~ x O895+ 


(—2.055 El, 


X10 El, 


8 22 


os 
Uy, = O ise = Use 


0 531+0 664 
822 


6E 
8 


x164+4 


x5iI, 
22 


os 


2x7 50€El, 2055 
— 2 . =0 458 El, 


< 


cs 


= ISOE 


_179 |4.00 


22. 


0.795+0 795 


2055==7 S50 EI, 


| 
bo 
ar ar 
wie 
822 
or or or 
ii. 
r 0-664 
i 
de~ Gne™ 
822 | 84 
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B22 B22 


=0 458 El, 


x 0.1875 ET, 
400 


(2x0458+4 x0 1875) I, 1666 El, 


2 
Qao™ Yeo™ 941: ™ © 


Qa2™ 1 666 0548 


4x xX 2055+8x 
Shear carry-over factors: 
General formula: 


P & 2x0795 
mj = (0531-0 664) = 


xX 2055 =-0697 
822 8 22 


548 OS548X184x2 


—~0112) x 0546= -0.430 
822 B22 ( 


Fixing forces: (Fig. 17) 


6000 Korn 


2H=—8x8 H=32t 32% 8+106S 


= Se =29802 


R=-3200¢ 


The successive approximations and the arrangement of the operations is 
tabulated in the Table VI. 
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GOLDEN GATE BRIDGE VIBRATION STUDIES# 


George S, Vincent,! M. ASCE 
(Proc. Paper 1817) 


SYNOPSIS 


The records of bridge movement obtained over a period of years by means 
of instruments at ten locations on the main and side spans of the Golden Gate 
Bridge are analyzed. They are plotted against wind velocity and correlated 
with theoretical analyses of wave form and against model test indications of 
the response to winds of different velocities. It is shown that the oscillations 
are almost always in a symmetric mode with the amplitude inhibited by the 
sharply downward deflection of the prevailing west wind as it passes over a 
high hill before striking the north side span. It is shown that during one forty 
minute period the oscillation was in an asymmetric torsional mode, uninhibit- 
ed, and reached amplitudes predicted from the model tests. A bottom lateral 
system has been installed. Tests indicate that it will prevent a reoccurrence 
of the objectionable torsional motion. 


INTRODUCTION 


The Inception of the Project 


The failure of the Tacoma Narrows Bridge in 1940 led to a close scrutiny 
by engineers and scientists concerned with the problem of all evidence of 
wind-forced oscillation of suspension bridges. Oscillation sufficient to be of 
interest had been observed on the Golden Gate Bridge although even in the 


Note: Discussion open until March 1, 1959. To extend the closing date one month, a 
written request must be filed with the Executive Secretary, ASCE. Paper 1817 is 
part of the copyrighted Journal of the Structural Division, Proceedings of the 
American Society of Civil Engineers, Vol. 84, No. ST 6, October, 1958. 

a. A report on a cooperative research project conducted by the Golden Gate 
Bridge and Highway District and the Bureau of Public Roads, U. S. Dept. of 
Commerce. A contribution to the work of the Advisory Board on the In- 
vestigation of Suspension Bridges. 

1. Bridge Engr., Division of Physical Research, Bureau of Public Roads, 
Washington, D. C, 
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severest storms in which a man could not stand erect the amplitude of oscil- 
lation had been insufficient to cause a computed stress of as much as 1,000 
pounds per square inch in the chords of the stiffening truss, which is less than 
the stress due to temperature variation alone. Little concern was felt for the 
safety of the bridge, but the Advisory Board on the Investigation of Suspension 
Bridges expressed interest in a cooperative research project between the 
Golden Gate Bridge and Highway District and the Bureau of Public Roads. A 
cooperative undertaking for this purpose became effective May 15, 1943. 
Under this agreement, instruments for measuring vertical movement of 
the bridge were developed and ten of them were installed upon the bridge. 
They were operated on a cooperative basis until June 30, 1950, at which time 
the cooperative agreement was terminated. All ten instruments have been 
operated continuously since that time at the expense of the Golden Gate Bridge 
and Highway District alone. The District has made available for the purpose 


of this report all data collected since the termination of the cooperative 
agreement. 


Early Vibration Observations 


The effect of earthquakes was taken into account in the design of the Golden 
Gate Bridge, and those responsible for its design and construction were alert 
during all phases of the work to the possible effects of vibration and reso- 
nance, Advantage was taken of the opportunity to measure any vibration which 
might occur at several stages of construction, using equipment and personnel 
provided by the United States Coast and Geodetic Survey.(1) Seismological 
instruments were set up on the piers, towers, cables and deck. 

The piers showed a weak vibration at about the frequency of the waves, and 
this was attributed to wave impact. They also responded to the vibration of 
the towers, 

Observations on the towers showed some variation in frequency from tower 
to tower, and at different times on the same tower, but with all of these there 
appeared rather broad bands of characteristic frequencies. In general, these 
frequencies responded to the changing mass and elastic properties as the 
towers were built up. 

As the cables were completed and other elements of the superstructure 
were erected, their contributions to the mass and overall stiffness were re- 
flected in altered frequencies. With the bridge completed, the instruments on 
the tower top showed a frequency of about 8.6 cpm (cycles per minute) for a 
movement parallel to the longitudinal axis of the bridge. This is near the 
frequency later observed for symmetric oscillation of the bridge which neces- 
sarily involves motion of the tower tops. Instruments on the cable and deck 
at mid-span showed a longitudinal oscillation of a few inches at frequencies of 
5.46 to 7.5 cpm which range includes the frequency computed later for the 
first asymmetric vertical mode. It should be noted that longitudinal motion at 
mid-span indicates the first asymmetric mode. 

The Coast and Geodetic Survey instruments showed also a lateral oscil- 
lation of small amplitude at mid-span at frequencies from 2.7 to 3.2 cpm. No 
such oscillation has subsequently been reported on the bridge. The instru- 


ments installed under the project herein reported respond only to the vertical 
component of motion. 
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Early Wind Excited Oscillation 


The Golden Gate Bridge, in common with other suspension bridges, exhibit- 
ed wind- or storm-excited oscillations during various stages of construction, 
from the cat-walk, used in cable construction, to the completed bridge. Only 
those occurring on the completed structure are considered in the present 
report. 

Early observations were approximate and incomplete. In those recorded 
before the special instruments were installed the period of vibration was esti- 
mated or measured approximately by pocket watch. Several wave forms were 
noted and remembered in accordance with the impression of the observer. In 
some instances amplitudes were estimated by sighting on bridge elements. At 
other times, sightings were taken with a transit on rods attached to the bridge. 

Mr. R. G. Cone, resident engineer during construction of the bridge, has 
recorded observations of two storms, one on February 9, 1938, some eight 
months after the bridge was opened to traffic, and the other on February 11, 
1941 (see Appendix I). 

The movements on February 9, 1938, evidently were a multi-noded vertical 
mode, with a phase difference which gave the effect of a running wave. No 
mention was made of evidence of a transverse or torsional oscillation. The 
high frequency observed, 20 to 30 cpm, corresponds to computations for the 
six-noded vertical oscillation (25.6 cpm)*, The computed loop length of the 
six-noded vertical motion averages 600 ft. (Cone estimated a motion of “as 
much as 2 ft. vertically in 300 ft. of bridge”), No other recorded observations 
describe movements having as high frequency or as short loops as these. 

The movement during the first recorded period of the storm of February 
11, 1941 suggests the second symmetric vertical mode having a computed 
period of 6.3 seconds (9.5 cpm) and a long end wave with the maximum ampli- 
tude near the quarter-point.**The change in frequency during the second re- 
corded period to about 7.5 cpm (8 second period) with dominant motion still at 
the quarter-point suggests a shift to coupled torsional and vertical asymmetric 
oscillation such as occurred on December 1, 1951, described hereinafter, ex- 
cept that Mr. Cone mentioned no evidence of torsional motion. 


Early Instrumental Observations 


A. Anemometer. In January, 1942 a recording anemometer was installed 
on the bridge. A wind velocity transmitter (pressure cup type) and a wind di- 
rection transmitter were installed on a steel tower at mid-span, approximate- 
ly 10 ft. west of the west cable and 5 ft. 3 in. above the top of the cable (14.5 
ft. above the roadway). These were wired to a recorder placed inside the east 
leg of the south tower. Both velocity and direction were continuously recorded 
against time in ink on a paper roll, each roll covering a period of two weeks, 
The time scale was 1 in. equals 20 minutes, 

B. Accelerometer. In March 1942, Professor Perry Byerly of the De- 
partment of Geology, University of California, was invited to visit the bridge 
and, on his recommendation, an instrument known as the Hall Recorder, 
loaned by the University to the Bridge District, was installed at the middle of 
the main span on the sidewalk. A photograph and description of this instru- 
ment have been previously published, (2) A brief description is in order since 
the instruments built later were patterned after it in some respects. 


*See computed frequencies and wave forms, Fig. 6. 
**See computed frequencies and wave forms, Fig. 6. 
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The Hall Recorder has two conical pendulums for measuring the two hori- 
zontal components of motion and a mass on a helical spring for measuring the 
velocity component. Each of the three has a weight of about twenty-five 
pounds. The spring suspension can be adjusted to give a considerable range 
of natural frequencies, permitting the selection of frequencies most favorable 
to the recording of the expected vibrations. The record is made by stylus on 
smoked paper on a drum which is revolved by clockwork in a spiral motion. 

A timing signal is also recorded on the drum at intervals. 

This instrument was installed on the west sidewalk at the middle of the 
main span on March 18, 1942, and remained until replaced by the new instru- 
ments designed to record only the vertical motion. 

Professor Byerly was employed by the Bridge District as a consultant to 
study and recommend the type of instrument to be used on the bridge. The de- 
sign was developed by the late L. H. Nishkian, Consulting Engineer for the 
District, Professor Byerly and Mr. Otto W. Hahn, Consulting Mechanical and 
Electrical Engineer, and the first instrument was completed in Mr. Hahn’s 


shop in December 1943, as the first major step under the cooperative 
agreement. 


Development and Verification of Accelerometer 


The Instrument 


A diagram of the instrument is shown in Fig. 1 and Fig. 2 shows a partial 
view of its final form. The weight hung on the helical spring vibrates verti- 
cally with a natural period, Tp. When the instrument is mounted on a body, 


such as a suspension bridge, which oscillates vertically with a period, Tj, 
and a steady amplitude, y, the sprung weight will oscillate at the forced peri- 
od, Ty, and will reach a steady state amplitude, a, measured with respect to 


the frame which supports the spring. The relation between a and y is ex- 
pressed by the equation:(3) 


a/y = 


4 92€ 
1.862 + 


In this equation € is the damping ratio which is the ratio of the first to the 
second of any two successive single amplitudes in the same cycle after the 
sprung mass has been displaced and released to oscillate under the influence 
of no external force or acceleration other than gravity. If there were no damp- 
ing € would be 1.0. 

When the spring weight moves a vertical distance, a, the horizontal arm 
attached to it will rotate about its pivoted end and the vertical arm, which 
carries the recording pen at its upper end and is rigidly attached at its lower 
end to the horizontal arm, will rotate through the same angle about the pivot. 


The movement of the pen, ap, is related to the movement of the sprung weight, 
a, by the expression 


= in which 


T> = distance from the pivot to the pen 
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Figure 2 - Accelerometer. 

Rocker arm and damper are pivoted to the 
spring-supported weight. Drive motor 
mounted on drum support is concealed by 
the drum. Relay for interrupting the 
recording is mounted on support partly 
under drum. 
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r = distance from the pivot to the center of sprung weight 


Then the ratio of the recorded amplitude, a,, to the amplitude of the oscil- 
lation of the supporting structure, y, is given by the expression: 


a/y = a/y (r,/r) (2) 


and the value of a,/y may be determined by multiplying the second term of 
Eq. (1) by To/r- The ratio, ao/y> is referred to as the total magnification. 

If T, exceeds T,, the displacement of the sprung mass is opposite in di- 
rection to the displacement of the spring support. If T,/T, = 2 or more the 
instrument is a “vibrometer” or displacement-measuring device. 4) Such 
instruments are used to record earth tremors. 

If, on the other hand, Ty exceeds To, the sprung mass moves always in the 
direction of the movement of the spring support; i.e., a and y have the same 
sign. If Tj/T, exceeds 2, the instrument acts satisfactorily as an “acceler- 
ometer” and the extension of the spring is nearly directly proportional to the 
acceleration acting on the instrument. Except for a phase lag caused by the 
damping the relation is nearly expressed by: 


acceleration = ®t g (3) 
ast 


in which: 


at is the instantaneous displacement of the sprung weight from its 
rest position 


Aagt is the static extension of the spring due to the sprung weight, and 
g is the acceleration of gravity. 


For the low frequency oscillations of the Golden Gate Bridge with computed 
periods as great as 10.5 seconds (5.7 cycles per minute; see Fig. 6) the instru- 
ment with its period, Ty, of 1.07 seconds acts as an accelerometer. Fig. 3 is 
a graph showing the ratio, ap /y plotted against the period, Tj, for different 
values of €. The ratio, ro/y of the instrument is 0.92. 

Since the acceleration acting on a vibrating body increases as the square 
of the frequency, higher frequency vibrations are recorded at a much greater 
scale than are low frequency vibrations. To prevent exaggerated recording of 
transient high frequency vibrations, such as are produced by traffic, damping 
is introduced in accelerometers. 

The drum upon which the record paper is wrapped is turned one revolution 
per hour by an electric clock motor through a train of gears. The axle is 
spirally grooved and supported on skewed knife edge rollers so that the drum 
advances longitudinally 0.4 in. per revolution. Twenty-six revolutions of the 
drum can be recorded on the 11 in. wide paper. The circumference of the 
drum is 37.43 in., giving a time scale on the record of 37.43/60 = 0.624 in. 
per minute. One of the gears of the driving train is mounted on a moveable 
bar which can be unpinned at one end to unmesh the gear and permit free ro- 
tation of the drum for adjustment or for damping measurements. 


Calibration 


The first model of the accelerometer was calibrated on a Scotch yoke at 
several frequencies in order to verify the applicability of Eq. (1) and Fig. 3. 
In the great majority of tests the actual and indicated amplitudes checked 
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within less than one per cent. Discrepancies of about 3% were recorded in 
three cases. Those were all for short periods of 2.9, 3.5 and 4.8 seconds for 
which the prediction is more sensitive to € as may be noted on Fig. 3. 

Of course, if the motion is not simple harmonic the determination of the 
amplitude of the motion of the bridge may be uncertain. However, since the 
instrument responds to the applied acceleration, a fairly correct interpre- 
tation of the record is possible for motions in which a more complex oscil- 
lation can be represented by a combination of two or three harmonic motions 
(See Appendix II). 


Damping Studies 


It is, of course, essential that the damping ratio, € , be known before Eqs. 
(1) and (2) or Fig. 3 can be used to determine the amplitude, y, of the oscil- 
lation of the test mount or the bridge from the measured amplitude, ap, of the 
pen record. A simple test to determine € was made frequently during the 
calibration runs and was made occasionally to check the damping of the instru- 
ments on the bridge. Fig. 4(a) illustrates the method. 

First, the drum drive was thrown out of gear and the drum was rotated by 
hand so that the pen described the centerline, A. Then, with the drum sta- 
tionary, the sprung weight was pressed down or lifted up till the pen was de- 
flected about an inch to one side as shown by B, and released. Point C repre- 
sents the amplitude reached by the pen one-half cycle after the release of the 
sprung weight. The damping ratio, € = a. 

Eq. (1) was derived on the assumption that the damping is viscous (the 
damping force being proportional! to the velocity of vibration) so that ¢€ is 
constant for all amplitudes, i.e., any two successive amplitudes of oppasite 
sign will have the same ratio as all other pairs of amplitudes. When the damp- 
ing was materially increased in order to cope with extraneous vibration it be- 
came desirable also to know whether or not the damping was viscous. The 
test for this is illustrated by sketches (b), (c), and (d) of Fig. 4. After the 
drum was rotated to draw the centerline, a straight pencil line, D, was drawn 
from one side intersecting the centerline at an angle of about 45°. Then the 
drum was positioned so that with the weight depressed an inch or so the pen 
would be on the pencil line. Upon release of the weight the longest record 
shown on any of the three sketches would be drawn by the pen. Then each of 
the shorter records just touching the pencil line was similarly drawn. With 
viscous damping and constant ¢€, the end points of the recording indicating the 
amplitude of the first damped half cycle in each test must fall on the straight 
line, E, as shown by sketch (b). 

If the damped amplitudes fell on a curved line as in sketch (c) it indicated 
that the damping effect was greater for large than for small amplitudes which 
would be true, for example, if the damping force varied as the square of the 
velocity of vibration. Records of this type were obtained when the damping 
was increased by means of discs which were nearly solid and nearly covered 
the cross sectional area of the damping-pot. 

It was thought that viscous damping could be obtained by exploiting fluid 
friction by forcing the damping oil to pass through a maze of small channels. 
A simple device for testing this theory was provided by sandwiching one or 
two copper “chore balls” of fine copper turnings between two perforated metal 
discs. In some tests with these the copper turnings bulged beyond the disc 
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FIGURE 4 —- SKETCHES ILLUSTRATING THE 
MEASUREMENT OF THE DAMPING RATIO. 


D 
C 
(a) (b) 
D 
¢ ¢ 
(c) (d) 
7 


ASCE VIBRATION STUDIES 1817-11 


and could be heard rubbing along the sides of the damping pot. In this con- 
dition the damping test record took somewhat the form of sketch (d) indicating 


the greater damping effect at lower amplitudes. This is recognized as charac- 
teristic of Coulomb friction damping. 


Operation of Instruments 


Installation 


Following the verification tests the first available instrument was operated 
for a few months in 1944 on the San Francisco, Oakland Bay suspension bridge, 
on the bottom chord of the north truss at the middle of the 2310 ft. main span. 
No strong winds occurred during this period and the only significant vibrations 
recorded were caused by the passage of trains. 

The instrument was then returned to the shop to facilitate preparations for 
making the additional accelerometers. Thirteen of these were made, nine for 
use on the Golden Gate Bridge, two for the San Francisco, Oakland Bay Bridge, 
and two for the Port of New York Authority. All were completed by the end of 
July 1945. Conduits and wiring circuits were installed on the Golden Gate 
Bridge and plywood houses to protect the instruments were placed at selected 
locations on the sidewalks. Nine instruments were installed in November and 
December, 1945. One instrument, No. 7, was kept at the shop and used for 
check tests which showed satisfactory performance. It was installed on the 
bridge March 1, 1946. The positions of the ten instruments are shown on Fig. 
5. Number One was the first instrument made and differs in non-essential 
details from the others. A synchronizing clock was placed inside the east leg 
of the south tower. This was set to operate relays to lift the pens of all instru- 
ments for five seconds at five-minute intervals. 


Operation and Development 


Routine servicing of the instruments occupied one man’s time for most of 
each forenoon. The accelerometer record sheets were changed daily and filed, 
and every two weeks the anemometer record roll was changed. The records 
were inspected regularly and memorandum notes were made of any sizeable 
amplitudes and the winds which caused them, but detailed analysis was de- 
ferred until many months of records were available. 

The maintenance personnel of the Golden Gate Bridge and Highway District 
constantly studied the behavior of the instruments to detect deficiencies and 
devise improvements in many details. The windows of the houses were 
covered with plywood to exclude the direct sunlight which caused considerable 
variation in temperature, thus altering the damping. Reflector lamps were in- 
stalled to maintain more uniform temperatures and prevent the absorption of 
atmospheric moisture which caused the paper on the drums to swell and 
wrinkle. Through steady experimentation the pens and their mountings were 
improved to correct defects which marred the earlier recordings. 

The variability of the damping between instruments and, over a period of 
time, of the same instrument, was sometimes enough to affect the interpre- 
tation of the record. After some experimentation with various weights of 
motor oil, a silicone fluid was substituted in 1947. This gave damping that was 
quite constant but rather high, nearly critical. Values of € of 100 to 1000 were 
used to interpret the record after this change in damping fluid. 
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In the fall of 1947 revised gear trains and stronger motors were installed 
to turn the drums and a spring action was provided in the meshing of gears to 
eliminate play and backlash and improve the uniformity of the record. 

When the time scale of the record was selected, it was thought that it would 
only be necessary to determine the frequency over several cycles with suf- 
ficient accuracy to obtain T; in order to compute the amplitude and also to 
correlate the observations with theoretical calculations of the frequency of the 
modes of vibration of the bridge. The scale was adequate for these purposes 
and also for identifying pure vertical or torsional motion by noting which 
stations moved in phase or in opposite phase. However, subsequent experi- 
mental work(5) and theoretical analysis(6) revealed that a truss-stiffened 
suspension bridge such as the Golden Gate Bridge might be expected to oscil- 
late in coupled vertical and torsional motion at an altered frequency and with 
a distinctive phase difference in the two motions. Evidence of such motion 
developed in the records, but the time scale was too small to permit the de- 
termination of phase differences. 

In 1953 at the direction of Clifford E. Paine, Consulting Engineer for the 
Bridge District, a change was made whereby the throwing of a switch, mounted 


in the southeast tower leg, would treble the speed of rotation of all of the 
drums, 


Observed Bridge Behavior 


General Observations 


Table I is a summary of major recorded bridge movements taken for the 
most part from the memorandum notes made from time to time by Mr. S. J. 
Dickson, Maintenance Superintendent, when he examined the records of wind 
and of bridge movement before filing them. General observations concerning 
the movements may be summarized as follows: 


a) The oscillation usually had a frequency of about 8 cycles per minute al- 
though small oscillations of other frequencies were occasionally noted. 

b) With the exception of the action on December 1, 1951, the predominant 
mode was symmetric, i.e., the wave form was essentially symmetrical 
about the center of the main span and, generally, the greatest amplitude 
was at the center. 

c) The motion was principally vertical, but on all occasions there was also 
a transverse torsional motion, sometimes very marked. It should be 
noted that torsional oscillation of very small magnitude could be readily 
detected by sighting across the deck and noting the angular movement 
with reference to a distant sighting point. The instrument records also 
revealed the more pronounced torsional oscillations. 

d) The amplitudes of vertical displacement were almost invariably greater 
along the leeward side than along the windward side of the bridge. This 
was noted also in laboratory tests of a section model of the bridge.(7) 

e) Oscillation was likely to occur in any wind above 20 to 25 miles per hour 
and the amplitude increased with an increase in wind velocity. 

f) On most occasions the amplitude increased sharply within a few cycles 


in response to a sudden increase in wind velocity and died down as 
quickly when the gust had passed. 
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TABLE I 
SUMMARY OF MAJOR RECORDED OSCILLATIONS 


Maximum 
Double* 


Amplitude Velocity Instrument Showing 
Inches mph Direction Prominent Movements** 


15 70 W. Midspan 
15 Th W. Midspan 


12 Midspan 


14 and 8 


10 


q 
a 
W 


/ 
/ 
/ 


Ss 
W 
Ss 


West 
amplitude usually occurred at Instrument No. 5 with a notable 


exception on December 1, 1951, when the maximum amplitude occurred at 
Instrument No. 3. 


** See Figure 5 for instrument locations. 


: 
1943 
12/9 
1945 
“10730 - 
1946 
1/5 
194 
1728 = 55 W/S.W. 5 
11/1 14 4o S.W. Sand 8 
1948 
2722 27 50 West 3, 4, and 6 
2/27 9 46 $/S.W. 5 
3/2 14 ho West 5 and 8 
3/19 16 30 West 3, 4, and 6 
5/13 16 34 West 3, 4, and 6 
5/17 27 46 W/S.W. 3, 4, 5, 6, and 8 
10/29 17 35 West 3, 4, and 6 
199 
3/25-28 13-18 32 West 5, 6, and 8 
4/12 14 30 West 5 and 8 
5/17 19 34 S.W./W. 3, 4, 5, 6, and 8 
10/5 21 38 W/N.W. 3, 4, 5, 6, ana 8 
1950 
te ks 45 to 56 West 3, 4, 5, 6, and 8 7 
10/27 a1 58 to 70 South 3, 4, 5, 6, and 8 
11/16 14 60 S.W. 5 and 8 
1952 
iio 39 65 WW. 3, 4, 5, 6, and 8 
4/29 25 .W. 5 and 8 
4/30 22 . 3, 4, 5, 6, and 8 
10/24 33 79 WwW. 3, 4, 5, 6, and 8 
12/1 132 69 emt all 
1952 
* The maximum 
i: 
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g) A wind of any velocity produced its maximum effect if it approached 
normal to the bridge axis. The amplitude for the same velocity was 
progressively less as the horizontal angle of approach deviated from 
90°. 

h) The maximum double amplitude was 45 inches except on December 1, 
1951. In that storm the amplitude was three times as great and several 
special features of the motion were observed, as described later in this 
report. The only significant structural effect of this motion was damage 
to the bronze facing plates of the wind shoes which necessitated re- 
placement of the plates. 


The detailed study of the records was directed toward three major ob- 
jectives: 


a) To correlate the observed motions of the bridge with its computed natur- 
al modes of oscillation; 

b) To determine the relation between the bridge motion and the velocity and 
direction of the wind which caused it; 

c) To correlate the observed behavior of the bridge with the behavior of a 
75-scale section model of it which was tested in the wind tunnel at the 


University of Washington on the initiative of the Washington Toll Bridge 
Authority. 7 


Some attention, also, was given to such evidence as could be obtained con- 
cerning the general character of the wind—its possible vertical angularity and 
degree of uniformity. 


Wave Form and Frequency 


In Fig. 6 are shown the wave forms and frequencies of several of the natur- 
al modes of motion of the Golden Gate Bridge computed by formulas developed 
in the cooperative research projects of the Washington Toll Bridge Authority, 
University of Washington and Bureau of Public Roads(8) and of the Oregon 
State Highway Department, American Institute of Steel Construction and 
Bureau of Public Roads.(6) The calculations took into account the vertical 
stiffness of the trusses, the extensibility of the cables and the torsional stiff- 
ness of the towers where these factors influenced the mode of motion. The 
torsional stiffness of the suspended structure without bottom laterals and the 
bending stiffness of the towers were found to be insufficient to influence the 
motion, In the study of the asymmetric oscillation which occurred on 
December 1, 1951, it was found necessary to develop an analysis taking into 
account the effect of the short hangers near the middle of the bridge. That 
analysis is described briefly in Appendix III. 

Table II shows the modes of oscillation most frequently identified on the 
bridge together with their observed frequencies and their computed vertical 
and torsional frequencies. The observed frequency induced by wind action was 
consistently between the computed vertical and torsional frequencies which 
agrees with model tests and the theory of flutter (the wind-forced coupling of 
vertical and torsional oscillation). 


Storm of December 1, 1951 


a) Instrument Record, In the storm of December 1, 1951, the motion of the 
bridge differed in several respects from that usually observed. The wind di- 
rection was quite steady, changing gradually from southwest at 3 p.m. to west 


October, 1958 


Ny = 9.5 Ny =12.1 


+ 


3RD SYM. | Ny =15.8 Nt =20.1 


+ 


| 4 TH SYM, Ny=18.6 23.7 


5TH SYM. Ny =25.6 Ny =32.6 


1ST ASYM Ny=5.7 
Ny =7.8* 


TOWER 
BRIDGE ¢ 


2 
* WITH HANGER PULL EFFECTHFREQUENGCIES IN CYCLES PER MINUTE 
Nv iS FOR VERTIGAL AND Nr FOR TORSIONAL OSCILLATION 
FIGURE 6- COMPUTED WAVE FORMS AND FREQUENCIES 
OF GOLDEN GATE BRIDGE WITHOUT BOTTOM LATERAL 


SYSTEM 


1817-16 ST 6 Po 

| | ST SYM Ny=7.5 Ny=96 
| | | | 

| 

| 

| } 


ASCE VIBRATION STUDIES 1817-17 


at 9 p.m. During mid-afternoon, with the wind averaging nearly 50 mph, the 
usual symmetric oscillation built up to about 2 ft. double amplitude at mid- 
span, being a little greater at the east or leeward side than at the west. The 
motion was somewhat random; i.e., the peak amplitudes rose quickly but sub- 
sided in two or three cycles and often were not matched by comparable motion 
at the opposite side of the roadway. Occasional bursts of the first asymetric 
mode showed at the quarter-points. From 4:45 to 5:00 p.m. the wind died 
down from a 60 mph gust to an average of about 44 mph and then rose steadily 
to a 64 mph gust at 5:15. The bridge became relatively calm during this inter- 
val but as the more violent part of the storm, 5:15 to 6:20 p.m., began a 
change in the mode of oscillation took place. Between 5:20 and 5:39 the asym- 
metric mode became dominant, beginning as a rather spasmodic and random 
movement but changing quickly to a steadily increasing oscillation of quite 
constant frequency (about 8.4 cpm). At the south quarter-point where there 
were instruments on both sides the amplitude was greater on the east 
(leeward) side, as indicated by instrument No. 3 than on the west (windward) 
side, instrument No. 4. Also the two sides were usually out of phase from 
one-fourth to one-half cycle, indicating a coupling of vertical and torsional 
modes. Such comparisons cannot be made for the north quarter-points be- 
cause there was no instrument on the east side. 

By 5:40 the double amplitudes at the south quarter-points of the main span 
had reached nearly 80 inches on the east and about 45 inches on the west and 
a new movement developed. While the asymmetric oscillation continued to in- 
crease, a symmetric oscillation at 16.8 cpm (double the frequency of the domi- 
nant asymmetric mode) appeared with well defined amplitudes at the centers 
of the main and side spans. In this symmetric mode the side spans generally 
moved at opposite phase to the main span movement but much of the time the 
two side spans themselves were not in phase* 

In contrast to this unsteadiness in phase along the bridge and to the strong 
element of torsion in the dominant asymmetric motion, the symmetric oscil- 
lation was almost pure vertical motion. The reason for this became 


*The small time scale of the records makes it impossible to be precise about 
phase observations. In a few instances during the more violent movements 
in this storm, the records for instruments Nos. 2 and 9 show obvious dis- 
tortion possibly due to slack in the drum controls and some of the variation 
in apparent phase may be due to small slip of this kind. The statements in 
the text concerning phase relations are based upon enough readings to show 
clear trends in spite of individual imperfections. 


TABLE II 
Co ison of C ited and Observed Frequencies (c 
Computed Frequency Observed 
Mode Vertical Torsional Frequency 
First Symmetric 7.5 9.6 8.0 
First Asymmetric 7.3 6.5 


Second Asymmetric , 14.3 12 to 12.4 
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evident when the cause of the symmetric oscillation was recognized as dis- 
cussed below. 

Two characteristics of the first asymmetric mode give rise to a lifting 
motion at mid-span twice during each cycle of the oscillation. First, the so- 
called node at mid-span is not stationary but must rise slightly at each ex- 
treme of the oscillation because of the distortion of the cable from a sym- 
metrical curve. It is this slight rise of the center of gravity of the bridge 
which stores the potential energy needed to maintain the oscillation. This 
slight distortion of the mode is shown in exaggerated form in sketch (a) of 
Fig. 7. For the Golden Gate Bridge a single amplitude of 5 ft. at the quarter 
point (equalled or exceeded from 5:50 to 6:10 p.m.) causes a computed mid- 
span rise of about 2 inches. 

The second and greater factor tending to lift the mid-point is illustrated by 
sketch (b) of Fig. 7. In the first asymmetric mode the mid-point of the cable 
moves back and forth longitudinally, being shifted toward the depressed 
quarter-point. For the Golden Gate Bridge this longitudinal displacement 
from the normal position is about 17 inches for a single amplitude of 5 feet at 
the quarter-point and is proportional for other amplitudes. If the suspended 
deck structure remained stationary longitudinally, the hangers would be pulled 
out of plumb. At the short middle hangers this would lift the suspended 
structure 1 inch for a 5 foot quarter-point amplitude. If, however, the sus- 
pended structure should move in the opposite direction to the cable motion 
and to the maximum limit of the slot of the wind shoe at the tower, the lift at 
mid-span would be 4.5 inches. (That such relative longitudinal motion did oc- 
cur will be demonstrated later.) 

Since the vertical movement of the midpoint is upward at both extremes of 
the dominant asymmetric motion, it follows that the two sides of the roadway 
must rise together whether the asymmetric mode is vertical or torsional and 
the 16.8 cpm symmetric motion must be a vertical mode as observed. 

These two factors cannot directly account for the single amplitudes of 8 to 
12 inches at mid-span which accompanied single amplitudes of 6 ft. at the 
quarter-points. However, this lifting of the suspended structure near mid- 
span exerts periodic forces on the cable which could build up a symmetric 
vertical oscillation whose natural frequency happens to be about 17 cpm. Fig. 
6 shows that the third and fourth symmetric vertical modes have frequencies 
near this value (15.8 and 18,6.) If the third symmetric mode had been excited 
(see Fig. 6) it should have shown its characteristic frequency and nearly maxi- 
mum amplitude on the records for instruments Nos. 4 and 6 at the southwest 
quarter-point and three-eights point, respectively, of the main span, whereas 
those records actually showed the dominant asymmetric mode frequency and 
their amplitudes were consistent with the wave form of that mode. On the 
other hand, the fourth symmetric mode as computed has its minimum ampli- 
tudes at the three-eights points and its relatively small amplitude at the 
quarter point would be obscured by the large asymmetric motion. It would, 
therefore, be revealed only at the middle of the main and side spans where the 
double-frequency oscillation was actually measured (on instruments Nos. 1, 2, 
5, 8, 7 and 10). This mode, with all parts of the bridge rising or falling simul- 
taneously, involves relatively high cable stretch.(9) Periodic vertical pulls 
confined to the middle portion of the span might be expected to excite just such 
a mode, The records (typified by No. 2 of Fig. 8) suggest strongly that the 
actual movement may have been a combination of two vertical modes at 
frequencies of 8.4 and 16.8 cpm since larger and smaller amplitudes alternate 
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Relative shift of deck and cable = 


(b) 


FIGURE 7 - SKETCHES ILLUSTRATING LIFTING AT 
MIDSPAN DUE TO FIRST ASYMMETRIC MODE. 
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and this supposition is well supported by the study shown in Appendix II and 
illustrated in Fig. A-1. The two modes involved might well have had the 
forms of the first and fourth symmetric. The phase of the side span oscil- 
lation differs in these two modes (see Fig. 6) so that the observed shifting of 
phase might be expected to occur. 

Fig. 8 shows portions of the records for December 1, 1951 from two of the 
instruments, Each record has several parallel hourly traces. Thus, for 
instrument No, 3 the records shown from 5:52 p.m. to 6:02 p.m. overlap to 
some extent the parallel records for the corresponding 10-minute intervals 
beginning successively at 1:52, 2:52, etc. until 7:52 p.m. 

The record for instrument No. 3 at the southeast quarter-point of the main 
span shows a uriform displacement of 0.98 in. along one side which was 
caused by the damper striking and resting momentarily upon the bottom of the 
damping pot during the downward displacement of the bridge. Due to the high 
damping of the instruments and the low frequency of the bridge in relation to 
that of the instrument, there was no rebound effect and the other half of the 
record, made while the bridge was above its normal position, was not affected 
and is therefore reliable. In view of the essential symmetry of the records, 
the downward displacements may be determined from the preceding and suc- 
ceeding upward displacements. 

The record for instrument No. 3 from 5:57 p.m. to 6:01 p.m. (interpolated 
between 5:52 and 6:02) shows an average amplitude for the upward dis- 
placement of 1.48 inches. Applying the magnification factor computed for the 
frequency of 8.4 cpm and critical damping indicates an average single ampli- 
tude of 6.2 feet. This average was equalled or exceeded sixteen times in this 
four-minute period, the maximum peak corresponding to a single amplitude of 
7.1 feet. 

Single amplitudes of 6.4 ft. and 4.0 ft. were recorded simultaneously sever- 
al times at the southeast and southwest quarter-points, respectively. They 
were frequently at opposite phases so that the transverse slope was 
(6.4 + 4.0)/84* indicating a single amplitude in torsion of 7.5°. 

Fig. 8 also shows a portion of the record for instrument No. 2 at the center 
of the south side span. This shows the symmetric mode at 16.8 cpm during 
the same period, 5:52 to 6:02 p.m. 


b) Other Indications of Movement. Marks on the paint on the fingered road- 
way expansion joints between the main span and the towers showed that, at the 
curbs, the main span had reached positions 27 inches and 25 inches nearer to 
the north and south towers, respectively, than it occupied six days later on 
December 7, 1951 when measured in practical calm and at a temperature of 
about 40°, considerably lower than on the day of the storm. These measure- 
ments should be corrected for the effect of horizontal rotation due to lateral 
deflection which was not measured or estimated on December 1, 1951. Ina 
similar storm on February 11, 1941, R. G. Cone reported the maximum later- 
al deflection as 5.6 ft. (See Appendix I). The end rotation caused by this much 
deflection would cause a longitudinal movement of about 2 inches at the curb 
line. Subtracting this from the measured displacement would leave 25 inches 
and 23 inches respectively. 

It was evident from the damage to the wind shoes that the bearing blocks 
which transmit the lateral forces from the main span to bearings attached to 


* The instruments were located on the sidewalks 3 ft. inside the trusses and 
42 ft. from the centerline of the roadway. 
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the towers had moved longitudinally the full distance permitted by the design. 
The plans indicate clearance for movements from the normal position of 18 
inches toward the tower and 21 inches toward the channel. Presumably these 
clearances are at 70°, assumed the normal temperature for the design.(10) 

At 40°, on the day of the inspection, the contraction of the suspended structure 
alone would increase the clearance toward the tower to about 23 inches. This 
agrees reasonably well with the marks on the roadway expansion joint noted 
above. 

Fig. 9 shows dimensions in the vicinity of the hangers at mid-span. The 
hangers pressed against the edges of the steel platform with sufficient force 
to crimp the edge angles. Mere contact would account for movements of 27 
inches in each direction, and the maximum movement must have been several 
inches more than this. Since the span oscillation was largely torsional, the 
two cables frequently moved in opposite directions and their pulls on the sus- 
pended structure were opposed. It must be presumed that the maximum 
movement of either cable relative to the deck occurred when its movement 
was opposite to that of the suspended structure. Its value would be the maxi- 
mum end movement permitted by the slot plus the negligible change in half 
span length due to the curvature of the vibration wave plus the longitudinal 
displacement of the cable center. The first two factors probably account for 
about 20 inches. The approximate maximum relative movement of the cable 


and deck at the center of span for various quarter-point amplitudes would then 
be: 


Quarter-point Amplitude Deck Movement Cable Movement Total 


5 ft. 20" 17" 37* 
6 ft. 20" 20.5" 40.5" 
7 ft. 20" 24" 44" 


Further indication of relative longitudinal motion of cable and deck ap- 
peared where light standards Nos. 65 and 75 rubbed paint off the east cable. 
These standards are 375 ft. south and north, respectively, of midspan at which 
locations the cable is at about the height of the standard as shown in Fig. 10. 
Angular displacement, explained below, brought them into contact with the 
cable with enough force to rub off the paint for a distance of 44 to 45 inches, 
extending toward mid-span from points directly opposite the normal positions 
of the standards. If it is assumed that the full width of the light standard made 
contact with the cable at the position of maximum displacement, 5 inches 
should be deducted from the length of the mark to determine the maximum 
relative movement. 

Fig. 11 illustrates how the mark on the cable near standard No. 65 could 
have been made. When the south half of the east cable moves down, the cable 
near mid-span moves south a distance which depends upon the quarter-point 
amplitudes as shown in the tabulation in the preceding paragraph. If the sus- 
pended structure also moved south the same distance contact would be made 
on the cable directly opposite the normal position of the standard at point ‘A” 
on Fig. 11 where one end of the mark on the paint was observed. When, on the 
other hand, the suspended structure moved in the opposite direction while the 
cable moved south, the contact point would be at the other end of the line at 

A relative rotation 5.9° was necessary in order to bring the light standard 
into contact with the cable as shown in Fig. 10. This rotation resulted from 
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the torsional oscillation and from the transverse slope of the hangers caused 
by the static pressure of the wind on the suspended structure as detailed 
below. 

The wave form of the first asymmetric torsional mode is essentially a sine 
curve whose equation is 6 = 99 sina x/2100, with x measured from midspan 
and 6. representing the amplitude at the quarter-point. At the light standard 
x = 375 ft. and @ = 0.53205. For 0, = 7.5°, a value often reached during the 
storm, 6 would be 4.0°. 

The transverse wind pressure against the suspended structure is resisted 
by the cabies and by the lateral truss. The division of the loading between 
them follows the general form of a probability curve, the truss sustaining 
100% of the load per foot near the towers and the cables sustaining more than 
100% (serving as an intermediate support for the truss) at mid-span. At 
some section near the center the cables just resist all of the wind load per 
foot. It is perhaps near the truth to assume that this is the case near these 
two light standards. The dynamic pressure of a 70 mph wind is 12.5 lb. per 
sq. ft. Assuming a drag coefficient of 2.8 for this wide structure gives a wind 
load of 2.8 x 12.5 = 35 lb. per sq. ft. With an exposed area of about 15 sq. ft. 
per lin. ft. the transverse wind load would be 525 lb. per lin. ft. The dead 
load of the suspended structure is 14,600 lb. per lin. ft. and the transverse 
slope of the hanger would be 525/14,600 = 0.036 corresponding to an angle of 
about 2°. 

Thus the static wind load and the torsional oscillation account for angular 
rotation of 2° and 4°, respectively, or a total of 6° which is about the minimum 
required to make contact between the cable and the lamp standard. The ad- 
ditional displacement needed at light standard No. 75 to stretch the anchor 
bolts 5/16 in. can be accounted for by one of a dozen movements of greater 
amplitude than used in the above calculations. 


Correlation of Bridge Movement with Wind 


A detailed study was made of the recorded vertical motion on the west 
sidewalk at the middle of the main span and the associated wind velocity and 
direction for the period December, 1944 to March 1946, and for a number of 
later individual storms. 

In Fig. 12 a graphic representation of the relation of amplitude to wind 
velocity and direction is shown on polar coordinate paper on which the points 
of the compass and the direction of the bridge axis had been previously indi- 
cated. An outline map of the vicinity of the bridge was also sketched in with 
the center of the bridge at the origin of coordinates. 

The wind velocity was scaled outwardly from the center and plotted on the 
radius representing the direction of approach of the wind. At this point the 
observed amplitude was recorded. A marked correlation of amplitude with 
wind velocity and direction became apparent so that connecting points of equal 
amplitude resulted in consistent amplitude “contours” as shown in Fig. 12. 

The contours show very forcefully the effect of wind direction. Moderate 
winds blowing at right angles to the bridge excite greater oscillation than do 
violent winds nearly parallel to the bridge. The contours cover little more 
than the quadrant between west and south; first, because the wind was usually 
from that quadrant, and secondly, because other winds had little effect on the 
bridge. During one storm (October 30, 1945) a south wind reached a gust ve- 
locity of 75 mph and ranged consistently around 50 mph; yet it produced a 
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double amplitude of only 12 inches.* A north-northeast and a south-southeast 
wind of 50 mph each produced only 2 inches of double amplitude. A 74 mph 
wind from north-northwest caused a double amplitude of 15 inches. 

It should be noted that if only the normal component of the wind velocity 
were effective in producing oscillation of the bridge the contours would all be 
parallel to the bridge axis. This is almost true of the actual contours, 

It was noted in plotting the points that occasionally a low amplitude from 
one storm would fall upon a higher amplitude from another storm of the same 
velocity and direction. Such variation arises from the degree of steadiness of 
the wind. In plotting this chart account was taken of the gustiness or steadi- 
ness of the wind and the velocity plotted was not necessarily the maximum 
peak velocity. In general, for fairly steady winds, the points were plotted at 
the gust velocities which were fairly frequent and uniform although much of 
the time the velocity might be 6 mph less than the peaks. For very gusty 
winds the point was usually plotted somewhat below the irregular and in- 
frequent peak gusts, It would appear that an unusually steady wind might give 
higher values than indicated by the contours. 

Fig. 13 shows to enlarged scale a graph of the envelope of the upper and 
lower limits of velocity taken from the anemometer record between 5:00 and 
6:40 p.m. on December 1, 1951. Also shown is the single amplitude of the 
movement of the southeast quarter-point for the duration of the asymmetric 
oscillation in the same period. The curve shows a running average of the 
amplitude neglecting the higher peaks. The characteristic frequency of this 
mode could not be identified on the record until the amplitude reached 10 
inches, This was due to the interference of the symmetric mode and random 
oscillation which preceded the asymmetric mode. 

It may be noted that the amplitude built up and diminished in response to 
the increase and decrease of the velocity but the response lagged some two 
minutes or more at several significant changes in velocity and showed only a 
moderate response to sharp gusts. This is to be expected in a fully developed 
mode of oscillation since the increment of energy contributed by a brief gust 
is not great in comparison with the total energy of vibration of the entire 
bridge. This is in contrast with the quick and brief response on previous 
records where the motion was somewhat random and not so much controlled 
by synchronized motion of the entire structure. 

Because of the time lag in the response of the bridge to changes in velocity 
and the reduced effectiveness of brief gusts the recorded amplitudes do not 
indicate accurately the relation between velocity and amplitude. However, a 
fairly consistent indication was found by first sketching from Fig. 13 the curve 
of average velocity and then plotting amplitudes against the velocity recorded 
2.5 and 5.0 minutes earlier (see Fig. 14). A smooth curve was then drawn 
through these points indicating the highest amplitude for a given velocity, 
following the reasoning that these more nearly indicated what could be built 
up at that velocity while the lower values showed the restricting effect of a 
too brief duration of the velocity for which the points were plotted. 

Fig. 15 shows response curves in which the amplitude of oscillation is 
plotted against the velocity of the wind which caused the oscillation. The co- 
ordinates shown are for the bridge itself. The curves from model tests have 
been adjusted to prototype coordinates by proper scale factors. 


*During this storm there was only the one instrument on the bridge, at mid- 


span. An asymmetric oscillation of greater amplitude could have occurred 
without detection. 
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Curve A was obtained by plotting amplitude against wind velocity as taken 
from Fig. 12 on a line extending westward from the origin at right angles to 
the axis of the bridge. It represents an average response curve of the bridge 
in the symmetric mode having a frequency of about 8 cpm. 

Curve B is the response curve for the first asymmetric mode of December 
1, 1951, taken directly from Fig. 14. The scale on the right indicates the 
torsional double amplitude involved in response curve B, assuming the verti- 
cal amplitudes of the southwest quarter-point at 0.6 of those plotted for the 


southeast (leeward) quarter-point. The actual ratio varied from about one- 
half to two-thirds. 


Correlation of Bridge Behavior with Section Model Tests 


Section model tests of the Golden Gate Bridge which were made in con- 
nection with the investigation of the Tacoma Narrows Bridge are described in 
Chapter VII of Part IV of the report on that bridge and in Chapter IV of Part 
v.(11) Curves, C, D and E of Fig. 15 are torsional response curves from 
section model tests plotted to the prototype scale for the asymmetric mode at 
3 different structural damping levels. Curve C is from Fig. 85 of Part IV, 
for B = 0° (horizontal wind) and 6 = 0.016 (logarithmic decrement of model in 
still air). Curve D is from Fig. 82(c) for 8 = 0° and 6 = 0.083. Curve E is 
from a test, the results of which have not been published. In this test 8 = 0° 
and 6 = 0.011. 

Curve F represents a reasonable interpolation of a section model response 
curve for a logarithmic decrement of about 0.031. If the structural damping 
of the Golden Gate Bridge in this mode were known a fairly accurate pre- 
diction of the maximum response curve for the prototype could be computed 
from model tests by the method described in Chapter IV of Part V of U. of W. 
Bulletin No. 116. Lacking data on damping, curve F represents the best avail- 
able basis for correlating the laboratory indications and the field observations. 


The difference between curves B and F can be attributed principally to two 
factors as follows: 


(1) Time required to build—The curves coincide at the lower amplitudes 
which can be built up rather quickly by the wind but the higher velocities 
did not continue long enough in the field to build the maximum ampli- 
tudes of which they were capable. A very rough idea of this factor may 
be gained from “Time-to-build” studies on the model setup from which 
curve C was taken.* These indicated that a wind just capable of build- 
ing a double amplitude of 18°, for example, required more than a minute 
to build to 16° from an amplitude of about 10°. The time scale of the 
model is the square root of its linear scale, in this case ¥'75 or 8.66. 
Applying this scale the time required for the same build-up on the 
prototype in a steady wind would be 10 minutes or more. It is evident 
from Fig. 13 that the fluctuations in wind velocity over a period of 10 
minutes must greatly reduce its effectiveness toward building up the 
maximum amplitudes attainable by a steady wind. 

(2) Dynamic damping—It is believed that the often opposing direction of the 
longitudinal motion of the cable and suspended structure had a strong 
dynamic damping effect on the dominant asymmetric oscillation, absorb- 
ing part of the energy which would otherwise have contributed to building 


*See Fig. 88 of Part IV, U. of Washington Bulletin No. 116. 
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up higher amplitudes. A significant amount of energy was also dissipated at 
the wind shoe bearing plates. In addition, the influences of this longitudinal 
motion as it changed phase vith respect to the asymmetric oscillation 
constantly altered the relative amplitude and phase of the vertical and torsion- 
al components of the coupled oscillation so that they could not maintain the 
optimum relation conducive to the most effective action of the wind. Further- 
more, the symmetric vertical oscillation of the side spans and part of the 
main span at an observed frequency double that of the dominant asymmetric 
oscillation, caused primarily by the action of the shorthangers near the center 
of the main span, was incongruous with the slower coupled motion which the 
wind tended to impose and thus must have increased the aerodynamic damping. 
It will be noted from Fig. 13 that this symmetric oscillation prevailed from 
5:40 to 6:30 including the period of high wind velocity when the motion might 
have been expected to reach greater amplitudes than did occur. 

In addition to the two sharply defined factors above, the effectiveness of the 
wind was further weakened by any non-uniformity of the wind velocity along 
the span, 

The thought naturally arises as to why the bridge oscillation in the more 
prevalent symmetric mode never approached the amplitude developed in the 
asymmetric mode on December 1, 1951, even though the wind velocity was 
often comparable to that storm. In 1945 the Washington Toll Bridge Authori- 
ty’s first section model test on this bridge indicated the possibility of much 
greater motion than had been recorded in the field. It was suspected that this 
might be a result of non-uniformity in the prevailing west wind which was de- 
flected sharply downward after passing over the high hill just west of the north 
side span. A simple laboratory test demonstrated the soundness of thistheory. 
The full model of the first design for the new Tacoma Narrows Bridge, with 
solid deck, was subjected to a wind condition which produced a vigorous oscil- 
lation. Then the vanes controlling the vertical inclination of the wind were 
adjusted over the greater portion of one side span to produce a sharp down- 


ward inclination. This had the effect of quickly reducing the oscillation to a 
small, somewhat random motion. 


Correlation with Bleich Flutter Analysis 


Dr. Friedrich Bleich applied the mathematical theory of flutter to sus- 
pension bridges(6) and extended it to take into account the modifying wind 
forces caused by the blunt edges, (12) He applied this to section model test 
data for the Golden Gate Bridge and predicted a critical velocity (at which 
flutter begins) of 9.2 fps at a frequency of 107 cpm for 6 (structural damping) 
of 0.003. Converted to prototype scale for the first asymmetric mode the 


critical velocity would be 38 mph, corresponding quite closely to curve C and 
E of Fig. 15. 


SUMMARY 


The observations on the Golden Gate Bridge have served to demonstrate 
the essential reliability of the formulas derived to show the wave forms and 
frequencies of the various modes of motion of a suspension bridge and also to 
give considerable confidence in the methods that have been developed for pre- 
dicting the behavior of suspension bridges in the wind by means of section 

model tests. This has proved true notwithstanding the presence of structural 
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factors not taken into account in the analysis and of non-uniformity and un- 
steadiness in the wind, factors whose influence can be judged qualitatively but 
not, as yet, on a quantitative basis. 

The provision for increasing the time scale of the recordings should in- 
crease their usefulness in verifying theoretical analyses of motion. Closer 
correlation with test indications might be expected if additional anemometers 
were provided and if the vertical angle of the wind were measured to give a 
more comprehensive conception of the wind structure associated with the 
measured bridge response. However, the bottom lateral system which was 
completed in 1954 is expected to restrict materially the torsional motion of the 


bridge and may reduce all oscillation to random movements of no signifi- 
cance.* 
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APPENDIX I 


Observations of Bridge Oscillation 


R. G. Cone 


Storm of February 9, 1938 (Report of January 7, 1941 at San Francisco) 


About 1:00 p.m. on the afternoon of February 9, 1938, a wind of unusually 
high velocity was blowing through the Golden Gate with the direction normal 
to the axis of the Bridge as nearly as could be determined. The force of the 
wind was so strong that it was impossible to stand erect on the sidewalk or on 
the roadway of the Bridge. 

The writer drove to the San Francisco tower in a closed car and was able 
to open the door on the leeward side and get out on the roadway. By crouching 
and standing in the lee of the west leg of the San Francisco tower, it was 
possible to cross the roadway. Sighting along one of the offsets in the tower 
to the Marin shore, it could be seen that the center of the Bridge was deflected 
between eight and ten feet from its normal position and was holding this de- 
flected position. It was also observed that the suspended structure of the 
Bridge was undulating vertically in a wavelike motion of considerable ampli- 
tude. The undulations were fairly rapid, in the neighborhood of 20 to 30 
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vibrations per minute. Because of their rapidity their amplitudes could only 
be estimated but it appeared that the stiffening truss was being distorted as 
much as two feet vertically in 300 feet of bridge. The wave motion appeared 
to be a running wave similar to that made by cracking a whip. The truss 
would be quiet for a second, and then in the distance one could see a running 
wave of several nodes approaching. The force of this wave was taken up in 
the movement of the expansion joint and the rocker arms at the top. 

While this movement was going on one of the electricians, Mr. F. L. 
Pinkham, came by driving a truck. The writer motioned to Mr. Pinkham to 
stop and get out of the truck. He did so and attempted to climb over the curb 
slightly north of the tower but the force of the wind blew him back over the 
curb and down onto the roadway. He crawled over to where the writer was 
standing in the protection of the tower, and was asked to observe the 
movements and actions of the Bridge, as a witness to substantiate what had 
been seen since the oscillations and deflections of the Bridge were so pro- 
nounced that they would seem unbelievable. Mr. Pinkham stood for some time 
observing the Bridge. 

The writer then decided to try and secure a record of some of the de- 
flections on film and returned to the office for his camera. Returning with the 
camera a few minutes later, he noticed that the roadway had stopped oscillat- 
ing but it was still deflected out of line. He then went to the top of the tower 
to get a record of the deflection in the camera but by that time the wind had 


diminished so that the deflection was not pronounced enough to show clearly 
in the picture. 


Storm of February 11, 1941 (Report of March 14, 1941) 
Direction of Wind: 


From southwest at an angle of approximately 45 degrees with the axis 
of the Bridge. 


Velocities at top of San Francisco Tower: 


1:00 p.m. 48 mph 
1:15 p.m. 60 mph 
2:50 p.m. 62 mph 
3:10 p.m. 60 mph 
3:30 p.m. 60 mph 


First Series Readings: 


12:45 to 1:15 p.m. interrupted by heavy rain obscuring sights. Transit 
was set up in shelter of east leg of San Francisco tower on a line 50 inches 


east of the east edge of the top chord. 12:45 p.m. foresight on same point, 
Marin tower. 


Deflections laterally: 


12:45 p.m. 1.75 ft. 
12:50 p.m. 2.30 ft. 
1:00 p.m. 2.90 ft. 
1:05 p.m. 4.90 ft. 
1:15 p.m. 4.00 ft. 
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During this period vertical oscillations were observed simultaneously 
at three points along the span by sighting on the stadia hairs in the transit. 
The period of these oscillations varied from 5-1/2 seconds to 6-1/2 
seconds for the full amplitude up and down. At the quarter point and at the 
third point (San Francisco) the magnitude of the oscillations was from 9 


inches to 20 inches as determined from using a 5-inch casting on the chord 
for a measure. 


Second Series Readings: 


2:50 p.m. to 3:30 p.m.—wind southwest. 


Latera: Oscillation 
Time Deflection Velocity Quarter point Center 
2:50 p.m. - 5.6 ft 62 mph 2.0 ft 1.6 ft 
3:10 p.m. - 4.6 ft 60 mph 1.8 ft 
3:15 p.m. - 2.9 ft 60 mph 1.8 ft 


The period of the oscillations were in the neighborhood of 8 seconds for 
the full up and down movement. 


APPENDIX II 


Analysis of Accelerometer Records of Secondary Symmetric 
Vertical Oscillation of December 1, 1951 


Fig. A-1 illustrates the relation between the trace of the accelerometer 
chart and the actual motion when the latter is a combination of two simple 
harmonic motions, the frequency of one being twice that of the other. 

Curve (c) showing the acceleration, Fig. A-1(b), to a reduced time scale, 
is similar to the charts for instruments 5, 8, 1 and 2 made during the asym- 
metric oscillation of the bridge on December 1, 1951 and illustrated by Chart 
No. 2 of Fig. 8. It was derived from a curve for the motion shown by the full 
line in Fig. A-1(a) which is composed of two simple harmonic motions one 
having a circular frequency, 2w, and amplitude, A, and the other a lower 
frequency, w, and a greater amplitude, B. It will be noted that the motion 
corresponding to acceleration curve (c) shows sharp amplitude peaks nearly 
equal to A + B at a frequency, w, alternating with brief periods when there is 
scarcely any movement. This is quite different from the amplitudes which 
would be obtained by applying Eqs. (1) and (2) or the curve in Fig. 3 at the 
frequency, 16.8 cpm, to the alternating high and low amplitudes of the instru- 
ment charts. Motions have been seen on small suspension bridges and on 
models which exhibited the delayed motion and quick whip illustrated by Fig. 
A-1. 

To analyze the charts by this method it is first necessary to determine the 
acceleration acting at the peak amplitudes. For this purpose Eq. (3) is suf- 
ficiently accurate. The actual displacement, a; of the sprung mass of the 
instrument with respect to its support is obtained by dividing the single ampli- 
tude of the recording, a,, by .92 which is the ratio of the pen movement to the 
mass movement by virtue of the mechanical linkage. The static deflection, 
i.e., the deflection of the sprung mass due to the acceleration of gravity is 
11.2" as computed from the known mass and the spring constant. Introducing 
these values in Eq. (3) the acceleration, @ is found to be 
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= Bo x 32.2 = 3.12 a, ft/sec® (a) 
in which ap is measured in inches. This acceleration may be substituted for 
2 

at (Fig. A-1(b) at the high peak and the low peak giving the two equations 
which can be solved for A and B, which determine the amplitude of the motion 
of the bridge, sketch (a), The equations depend upon the relative amplitudes, 
A and B, and must be solved by trial. In studying the motions of points 5, 8, 
1, 2, 7 and 10 at about 6:08.4 during a period of relatively smooth oscillation 
it was found that the equations were roughly 


large peak, QC = w?(4a + .77B) (b) 
small peak, = - 77B) (c) 


in which w is the circular frequency, 27N, which in this case is *A x2 = 68. 


The circular frequency, 2w, corresponding to the higher frequency, 16.8 cpm, 
is 1.76. 

Values of @ for the large and small peaks at 6:08.4 were computed for 
points 5, 8, 1, 2, 7 and 10 and these were substituted in Eqs. (b) and (c) to ob- 
tain the values of A and B, which were then used to determine the amplitudes 
of the vertical oscillations of the bridge at the points indicated. The results, 
shown in the tabulation below, are reasonably approximate indications of the 
magnitude of the symmetric oscillation at the time selected for comparison. 


At approximately that time the single amplitudes outlining the dominant 
asymmetric mode were 


Station 4, South 1/4-pt. main span west 54" 
Station 6, South 3/8-pt. main span west 43" 
Station 9, North 1/4-pt. main span west 57" 


The measurements made are not actually simultaneous—the time scale is 
too small to permit such measurements—and it is not known which of the high 
or low peaks of the tabulation may have synchronized with the main sine wave 
outlined quite well by stations 4, 6, and 9. It is seen, however, that the maxi- 
mum amplitudes at the center of the main span were about one-fifth of the 
maximum quarter-point amplitudes of the asymmetric mode, 


Station A B Approximate Amplitudes 
(inches (inches) 
Large Small 
5 center main-east 6.5 7.5 12 2 
8 "west 5 8 11 0.5 
1 side-east 6 12 16 
2 " west 5 10 13 
7 "N. east 3.20 7 9 
10 6 10 14 0 


| 
—| 
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APPENDIX III 


Frequency of Asymmetric Mode as Affected by 
Longitudinal Component of Hanger Pull 


The hangers are long enough to permit relative longitudinal movement of 
the cable and deck in either direction but they oppose the relative motion by 
the longitudinal component of the hanger pull when inclined from the vertical. 
If the first asymmetric mode of oscillation is pure torsion with the cables 
moving in opposite directions and their longitudinal motions therefore op- 
posed, the deck must either remain stationary with the hangers sloping in 
opposite directions to the two cables, or it must continuously change its 
motion with respect to the cables. If it swings with one of the cables, it will 
not oppose the motion of that cable but its opposition to the other will increase 
and this difference will result in the two sides oscillating at unequal frequen- 
cies so that their phase relation must change. Such action is indeed evident 
in the records for the storm of December 1, 1951, and considerable variation 
in frequency is found if only one or two cycles are measured. However, the 
frequency when averaged over 8 to 30 cycles is quite consistently found to be 
8.4 cpm which suggests that the average effect was that of the deck remaining 
stationary so that the hangers opposed the longitudinal motion of both cables. 

An approximate solution for the frequency of the first asymmetric mode 
was obtained by using Eq. (4.10) from Part Il, University of Washington Engi- 
neering Experiment Station Bulletin No. 116(8) after substituting for h the 
computed value of the sum of all the hanger horizontal components for half of 
the main span. To compute these horizontal components it was necessary to 
obtain a formula for the longitudinal displacement, { ,, at each hanger lo- 
cation. The formula was derived from Eq. (b) on page 15 of Part Il. It was 
solved for each hanger. Then the slope of each was computed as the dis- 
placement divided by the hanger length and the longitudinal component of the 
dead load pull in each hanger was computed. Summing these up for the half 
span, h was found in terms of 19, the quarter-point amplitude, and substituted 
in Eq. (4.10) which was then solved for 2 from which the frequency was com- 
puted from the relationships given in Part II. 

In Eq. (4.10) h is applied to the cable at midspan as a concentrated longi- 
tudinal load. This condition is approximated in the present analysis because 
the longitudinal component of the hanger load drops off rapidly as the hanger 
length increases and most of the horizontal force comes from a few hangers 
near midspan, 

The computed frequency for the torsional mode was 7.75 cpm which is still 
8% less than the observed frequency. Other elements of resistance not con- 
sidered in this approximate analysis must account for the difference. One 
such neglected element is the torsional resistance of the tower, not ordinarily 
considered a factor in the asymmetric modes but brought into limited play in 
this case by the action of the longitudinal component of the hanger pull. 
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SYNOPSIS 


Tests of pairs of large beams with conventionally spaced and with bundled 
longitudinal reinforcement are reported. The bundles of reinforcement used 
comprised groups of four No. 6, four No. 8 or three No. 9 touching bars. 
Pairs of beams were compared with respect to width of flexural cracks, 
steel stress distribution, deflection, and ultimate strength. No significant 
difference in behavior or ultimate strength was found for bundled as com- 
pared to spaced reinforcement. 

Tied columns were tested by concentric loading to compare spaced and 
bundled longitudinal reinforcement consisting of twelve No. 6 or twelve No. 8 
bars. Comparison with respect to ultimate strength indicated that bundling is 
a safe detailing procedure when adequate ties are provided, even for 6.6 per 
cent longitudinal reinforcement. Splicing of bundled reinforcement in columns 
was explored and found feasible. 


INTRODUCTION 


Bundled reinforcement in structural concrete refers to reinforcement 
placed in groups of touching bars. As compared to the minimum bar spacings 
commonly used in beams, for instance those given by the 1956 ACI Building 
Code Requirements for Reinforced Concrete, ACI 318-56 Section 505(a), 
bundling permits the necessary bars to be placed in much narrower sections. 
As a result, bundling permits construction of lighter, more graceful and more 
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economical beams of box, channel or T-section. In beams of normal width, 
the clear distance between bundles will be considerably greater than the dis- 
tance between individual evenly spaced bars. Particularly when heavy nega- 
tive moment reinforcement must be embedded in the top of beams, bundling 
greatly facilitates concrete placement and insertion of spud vibrators. In 
columns, bundled reinforcement permits a reduced concrete cross section, 
which may be an important advantage in the lower stories of tall buildings. 
Bundling also permits interior ties to be omitted, so that concrete placement 
is facilitated. Finally, bundled bars in beams and columns may be a Satis- 
factory alternate solution for large sizes of specially rolled reinforcing bars 
which are occasionally used in very large structures. 

Practical use of bundled reinforcement in beams has been pioneered by 
Mr. H. M. Hadley; the list of references at the end of this paper indicates 
that he has reported extensively regarding experience with bundling. Several 
structures with bundled reinforcement have been built,(2,4-7) and good service 
records have been reported. 

Laboratory tests that have been reported(1,3) concern principally bundles 
of four 1/2-in. square bars in beams, and there may be some question regard- 
ing the performance of bundles of larger bar sizes. No tests of bundled rein- 
forcement in columns have been reported. An experimental investigation was 
therefore carried out in the Research and Development Laboratories of the 
Portland Cement Association during 1955-57 to investigate the performance 


of large deformed bars placed in bundles as longitudinal beam and column 
reinforcement. 


Test Beams 


As compared to spaced bars, bundling may be questioned primarily with 
respect to the bond integrity of beams. The most serious conditions may then 
be expected for bars placed as negative reinforcement near the top of deep 
and short beams. Previous tests have clearly indicated that, due to adverse 
effects of settlement, the bond resistance of top bars is less than that of bot- 
tom bars; and a short beam span leads to high bond stress before the flexural 
ultimate strength is developed.(8) The test specimens for this investigation 
were therefore short, deep beams, with the tension steel at the top as cast. 

In the beam designations, the first number shows the number of bars and 
the second number their size; the letter S indicates spaced bars and B indi- 
cates bundled bars; H indicates high-strength steel. 

The five test beams with spaced reinforcement shown in Fig. 1, Beams 
8-6S, 8-6SH, 8-8S, 8-8SH and 6-9S were designed by first determining the 
minimum beam width for a chosen group of bars. By ACI 318-56, 505(a), 
this width is governed by a minimum protective cover of 1-1/2 in. and, for 
1-1/2 in. maximum size aggregate, a clear distance of 2 in. between parallel 
bars. The beam depth was chosen so that the ratio of reinforcement was 1.5 
per cent. Finally, the distance from the face of a centrally located column 
stub to the beam support was chosen as twice the effective beam depth. The 
beams and all bars were extended 6 in. beyond the supports. The gross con- 
crete dimensions of Beams 8-6S, 8-8S and 6-9S excluding the column stubs 
were 13x21x97 in., 14.5x33.5x146.5 in. and 11.5x38.8x163.8 in. 

The beams with bundled reinforcement, Beams 8-6B, 8-8B and 6-9B were 
identical to the corresponding beams with spaced bars except for the bar 
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arrangement. The effect of decreasing the beam width for bundled reinforce - 
ment was investigated through Beams 8-6BH and 8-8BH. For these two 
beams, and their companions with spaced reinforcement, high-strength rein- 
forcement was used to delay flexural failure and develop very high bond 
stress. 

The column stubs of all beams were reinforced with four bars of the same 
size as the longitudinal beam reinforcement, and these bars were extended 
through the beam. Vertical stirrup reinforcement was provided as shown in 
Fig. 1 to prevent diagonal tension and shear failures. The stirrups also 
served the function c. preventing horizontal splitting that might otherwise 
have been caused by high bond stress. Beams 8-6SH and 8-6BH had two No. 6 
bars placed as compression reinforcement to prevent flexural compression 


failure, Beams 8-8SH and 8-8BH had two No. 8 bars placed as compression 
reinforcement. 


Properties of Bundles 


The external perimeter for a bundle of four bars as shown for Beam 8-6B 
and 8-8B in Fig. 1 is 37D, that is, 25 per cent less than for the same bars 
spaced in the usual manner. On the other hand, a single very large bar with 
the same cross-section area as a bundle of four bars with diameter D would 
have a diameter of 2D and a perimeter of 27D. Accordingly, the bundle of 
four bars has an exposed perimeter 50 per cent greater than that of the single 
large bar. The bundle of three bars used for Beam 6-9B similarly has an ex- 
posed perimeter 16.7 per cent less than that for the same spaced bars, and 
44 per cent greater than that of a larger bar of the same area. 

These geometric properties indicate that bundling leads to only a moderate 
increase in bond stress as compared to spaced bars. Replacing a single very 
large bar by a bundle with the same area, leads to a reduced bond stress. It 
should be noted, however, that lug height by ASTM A-305-53T would be greater 
for a single very large bar than for a bundle of bars, and bond resistance for 
top bars is known to increase with increasing lug height.(8) 


Materials 


A laboratory blend of Type I cements was used. Sand and gravel aggre- 
gates were combined to gradations within the limits given by ASTM C-33-55T 
for 1-1/2 in. maximum size. The concretes were mixed in 6 cu ft batches, 
and had a slump of 2-1/2 to 3 in. The cement content was 4.0 to 4.2 bags per 
cu yd, and the w/c ratio by weight was 0.62 to 0.70. The compressive 
strengths reported with the beam test data in Table 1 are averages of twelve 
to sixteen 6x12 in. cylinders made, cured and tested with each beam. 

All reinforcing steel used conformed to ASTM A-305-53T for deforma- 
tions and had the yield points given in Table 1. Lug height was 0.048, 0.055, 
0.067 in. and lug spacing was 0.35, 0.48, 0.50 in., respectively, for No. 6, 8, 
and 9 bars. The steel was purchased in such lengths that each bar vielded 
two beam lengths and a test coupon. By using one of these lengths in each 
beam of a pair, variation in steel properties between pairs was minimized. 
The No. 3 stirrups used in beams with eight No. 6 longitudinal bars had a 
yield point of 52,100 psi, the No. 4 stirrups used in the remaining beams 


48,700 psi. All reinforcement for each beam was tied together without weld- 
ing. 
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Casting 


All beams were cast in plastic-coated plywood forms with the column stub 
down so that the longitudinal reinforcement was embedded near the top as it 
is near the columns in continuous beams. Concrete was placed by spud vi- 
brators in the test beams and their companion 6x12 in. cylinders. It was 
found, of course, that bundling of the bars eased the placing operation sub- 
stantially. After testing the beams it was found that mortar had penetrated 
into and filled the cavity between the bars of all bundles. The test beams and 
their companion cylinders were cured one week under moist burlap. The 
beams were then turned upside down and stored in the laboratory until they 
were tested at the ages given in Table 1. 


Test Method 


All beam specimens were tested as shown in Fig. 2, supported on a roller 
at one end and on a rocker at the other, with a 6-in. overhang at each end. 
The testing machine load was applied in increments to the column stub 
through a 2-in. steel plate. The total duration of each beam test was about 
two hours. 

Deflection dial gages were mounted directly below the column edges and 
mid-way between these edges and the supports. The widths of all cracks 
were measured by a graduated microscope at the level of the centroid of the 
longitudinal reinforcement. 

To minimize the bar surface area isolated from bond by water-proofing of 
electric strain gages, SR-4 Type A-12 gages were placed in the intermediate 
grade bars in milled slots 3/32 in. wide, 3/8 in. deep, and about 6 in. long. 
The high-strength bars could not be milled. A gage was cemented to the side 
of each slot, lead wires were attached, and the slot was filled to the bar sur- 
face with wax waterproofing. Tension tests indicated that gages so placed 
yielded measurements in close accord with mechanical strain measurements 
over the same reduced section. The stress at any bar load was 5 to 9 per 
cent higher in the slotted section than in the full bar section. Eight strain 
gages for each beam were located as shown in Fig. 1. Two gages were 
placed on each of four bars symmetrically about mid-span. The measured 
strain without correction was assumed to represent the average strain in all 
bars at the location of the gage. Therefore, measured strain is reported 
herein as stress obtained by multiplying the average strain in the two half 
spans by the modulus of elasticity for the full section of the various bar 
sizes as obtained in tension tests. 


Beam Test Results 


All beams with intermediate grade steel, Beams 8-6S through 6-9B in 
Table 1, failed by yielding of the longitudinal reinforcement followed by large 
deflections and final crushing of the concrete compression zone at the column 
face. As shown in Fig. 2, both flexural and diagonal cracks tended to extend 

upward toward the corner at the column stub so that it was hardly possible to 
differentiate between “flexural” and “diagonal” cracks. No indication of bond 
failure was found in any of these beams, and no visual difference in behavior 

was noted for beams with bundled as compared to spaced bars. 
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Three beams with high strength reinforcement failed in bond as indicated 
by large amounts of bar slip at the beam ends; for Beam 8-6SH the tension 
steel yielded following bar slip at the beam ends. 


Steel Stress and Deflection 


Measured steel stress and deflection are shown for the three beam pairs 
with intermediate grade steel at various load levels in Figs. 3, 4, and 5. Both 
steel stress and deflection are given as an average of two measurements 
symmetrical about midspan for each beam. 

Distribution of measured steel stress along a longitudinal reinforcing bar 
in a beam specimen may be expected to reflect bond distress. Preceding a 
final destruction of bond, an abnormal rise of steel stress should take place 
toward the beam ends. Figs. 3, 4, and 5 show that the distribution of steel 
stress was closely alike for the two members of each pair of test beams even 
at high steel stress. It may be noted, on the other hand, that the steel stress 
for all beams was practically uniform at high loads in the middle third of the 
span. This was certainly caused by a stress redistribution resulting from 
the deep-beam type of crack pattern shown in Fig. 2. It is also seen that the 
overhangs contributed to the bonding action because the steel stress of all 
beams is not zero over the supports at high loads. It is felt that this be- 
havior is related to local stress disturbances in the support region where 
heavy reaction forces entered the abnormally short beams. 

The deflection curves are also similar for all beam pairs. Accordingly, 
both steel stress and deflection measurements indicate that there was no 
significant difference in behavior between bundled and spaced reinforcement. 


Crack Width 


Crack patterns were closely similar within pairs for all tests. Bond dis- 
tress may be expected to open up a few wide cracks near the beam ends 
rather than to increase the width of all cracks. Crack widths are therefore 
given in Fig. 6, for the beams with intermediate grade steel, as the average 
width of the three widest cracks in the beam. Steel stress is given in the 
figure as values calculated from applied moment at the column face section, 
taking the internal moment arm as 7/8 times the effective depth. It is seen 
that there is no systematic difference between the crack widths for bundled 
and for spaced reinforcement. Furthermore, no crack opened suddenly be - 
fore yielding of the reinforcement was in progress. This indicates that even 
the high local bond stress which acts near cracks resulted in only the normal 
minor bond slip for bundled as well as for spaced reinforcement. 

For the four beams with high strength reinforcement, a similar lack of 
systematic difference was observed between crack widths for bundled and 
spaced reinforcement. However, for Beams 8-6BH, 8-8SH, and 8-8BH, a 
few cracks near the beam ends became very wide shortly before final bond 
failure took place as the ultimate load was reached. 


Flexural Strength, Beams with Intermediate Grade Steel 


It is seen in Table 1 that some of the beams with intermediate grade steel 
carried loads considerably above their yield loads, which yield loads are 
listed as detected by strain and crack width measurements. The calculated 
flexural ultimate loads, P.ajc-, were obtained by the equation for ultimate 
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internal moment M, = bd2t'. q(1-0.59q) given in ACI 318-56, A605(b), con- 
sidering the external moment at the face of the column. The ratio of 
measured to calculated ultimate load exceeds one for all beams. The aver- 
age ratio for the three beams with spaced bars is 1.13, and the average ratio 
for the beams with bundled bars is also 1.13. This indicates that there was 
no systematic difference in ultimate flexural strength developed by spaced 
and by bundled bars except that the beams with bundled bars were slightly 
stronger by virtue of the slight increase in effective beam depth. 

The average ratio of 1.13 also confirms previous findings that the equation 
for ultimate moment, which was developed essentially by tests of small beams, 
is indeed applicable also to the large beams of this investigation. It is be- 
lieved that the excess of measured ultimate loads over the calculated load re- 
sulted principally from strain hardening of the reinforcement. A biaxial state 
of stress at the column stub appeared to delay crushing of the compression 
zone so that large steel strains were developed locally at the major flexural 
cracks. 

Bond stresses are given in Table 1 as calculated at ultimate load by divid- 
ing the shearing force by the external perimeter of the bars times 74/8. These 
bond stresses, for the beams with intermediate grade steel, were sustained 
without any indication of bond failure. They are in no way to be regarded as 
ultimate bond stresses. To develop higher bond stresses with intermediate 
grade steel it would have been necessary to make special test beams with part 
of the tension zone removed, or to make the beams so short that they would 
act as walls rather than beams. Both of these cases were thought not to 
represent practical conditions under which bundled bars may be used. High- 
strength steel was therefore used to study ultimate load stress. 


Bond Strength, Beams with High-Strength Steel 


Table 1 shows that the beams with high-strength reinforcement failed at 
ultimate loads close to the flexural strengths calculated by ACI 318-56, 
A606(a). Beam 8-6SH failed in flexure after bond slip had been observed at 
the beam ends. The remaining three beams failed at loads below the com- 
puted flexural ultimate. Failure was in bond, as indicated by large amounts 
of bar slip observed by dial gages as a relative movement between bar ends 
and the concrete surface at the beam ends. 

Bar slip is plotted as a function of calculated bond stress in Fig. 7. Bond 
stresses at ultimate strength, calculated by dividing shearing force by ex- 
ternal bar perimeter times 7d/8, are also shown. It is seen that bond slip 
was in progress when Beam 8-6SH failed in flexure at a bond stress of 520 
psi. By comparison with the slip records for Beams 8-8SH and 8-8BH in 
Fig. 7, both of which failed in bond, it must be expected that Beam 8-6SH 
would have failed in bond at a stress only slightly greater than 520 psi if 
flexural failure had been prevented by a yet higher steel yield point. Hence, 
the ultimate bond stress for spaced No. 6 bars must be expected to exceed 
only slightly the value of 513 psi observed for bundled bars. For No. 8 bars, 
the ultimate bond stress for spaced bars was 337 psi, which is slightly less 
than the stress of 391 psi observed for bundled bars. It should be noted in 
Fig. 7, however, that bond stress for a given slip value was always lower 
for bundled than for spaced bars. 

It can be concluded that, when only external bar perimeter was used to 
calculate bond stress, there was no systematic difference in ultimate bond 
stress developed by spaced and bundled hars. 
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BEAMS AND COLUMNS 


Test Columns 


The beam tests indicated that bundling of tension reinforcement is a satis- 
factory detailing procedure. A series of 10 tied columns was then designed to 
study bundled compression reinforcement. Concentric loading was chosen. 

An outline of the test program is shown in Fig. 8. All columns were 12x12 
in. with a height of 6 ft. Two amounts of longitudinal reinforcement were 
used, namely, 6.58 and 3.67 per cent made up by 12 No. 8 and 12 No. 6 bars, 
respectively. The 1/4 in. tie diameter used is the minimum permitted by 
ACI 318-56, 1104 (c). The corresponding maximum tie spacing of 48 tie 
diameters is 12 in., which is also the maximum spacing as governed by the 
12-in. column size and by 16 times the diameter of the No. 6 bars. 

Six columns with 12 No. 8 bars were tested. Column 12-8S contained bars 
spaced in the regular manner and surrounded by a square tie. The interior 
bars were held firmly by two interior rectangular ties. All ties were spaced 
12 in. Column 12-8B-1 contained bars bundled at the corners. Interior ties 
were Omitted, and the exterior tie spacing was maintained at 12 in. For 
Column 12-8B-2, the exterior tie spacing was decreased to 6 in. A splice 
was provided at mid-height of Column 12-8B-3 as shown in Fig. 8. The 
spliced bars were cut by a saw, and each bar was touching its longitudinal 
extension. The bars of Column 12-8B-4 were cut at the splice with a hydrau- 
lic bar cutter so that the bar ends were wedge-shaped. A clear spacing of 
1/4 in. was provided between the two parts of each longitudinal bar. It should 
be noted that the splice lap is only 10 bar diameters as compared to the mini- 


mum amount of 20 diameters given by ACI 318-56, 1103(c). A similar group 
of six columns with 12 No. 6 bars was tested. 


Materials 


A laboratory blend of Type I cements was used. Sand and gravel aggre- 
gates were combined to gradations within the ASTM C-33-55T limits for 3/4- 
in. maximum size. The mix ratio, cement: sand: gravel, was 1:3.58:2.38 by 
weight, and the w/c ratio by weight was 0.64 to 0.68. 

Two concrete batches were used for each column. In previous column 
tests(9) it has been found that failure generally takes place near the top of 
vertically cast columns. To explore this phenomenon, the bottom batch of 
some columns was made with a slightly higher w/c-ratio than the top batch. 
Compressive strengths representing averages of three to four 6x12 in. 
cylinders for each batch made and cured with the corresponding columns are 
given in Table 2. 

All reinforcement was intermediate grade steel and was tied into cages 
without welding. The ties were 1/4-in. plain bars, the longitudinal reinforce- 
ment conformed to ASTM A-305-53T for deformations and had the yield 
points reported in Table 2. All reinforcing bars were cut by a saw to a length 
tolerance of 1/32-in., and 3/4-in. bearing plates were placed touching the 
bars at the top and bottom of the columns. The lower plate was placed in the 
form before casting, the upper plate was set in a thin layer of high-strength 
plaster after the concrete was cured. The heavy tie reinforcement shown in 


Fig. 8 prevented failure at the column ends by possible local non-uniform 
stress conditions. 
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Casting 


All columns were cast in a vertical position in plywood forms protected by 
an epoxy resin paint. Concrete was placed in columns and companion cylinders 
by spud vibrators, and it was noted that the absence of interior ties in the 
columns with bundled reinforcement eased the concrete placing operation 
substantially. By inspection after testing the columns, it was found that 
mortar had filled the cavity between the bars of all bundles. The test columns 
and their companion cylinders were cured four to five days under wet burlap. 
They were then stored in the Laboratory until they were tested at the ages 
given in Table 2. 


Test Method 


All columns were tested under concentric loading as shown in Fig. 9, with 
both ends fixed against rotation. Spherical bearings permitted rotation at both 
column ends until a load of 20 kips was applied, and the bearings were then 
blocked by steel wedges. Electric strain gages applied at mid-height of all 
four column faces were monitored by a continuous strain recorder. The 
spread between the four gage readings was less than 15 per cent even at ulti- 
mate strength, which indicates that a closely concentric loading was obtained 
for all columns. In addition to electric strain measurements, the total 
shortening over the entire column height was measured by a dial gage. A 


continuous loading speed of 160 kips per minute was maintained for all 
columns. 


Column Test Results 


The column test results were evaluated essentially in terms of the equation 
for ultimate strength of concentrically loaded tied columns established during 
the ACI Column Investigation in the 1930’s, Py = 0.85 f, (Ag = Agt) + Ast fy- 
This equation has been confirmed by several recent investigations(9) and 
is used in Section A608(b) of ACI 318-56. A comparison of measured and 
calculated ultimate loads is given in Table 2 together with concrete and steel 
properties. 

The test data were also studied in terms of the relationship between ap- 
plied load and total column shortening expressed as strain. The load- 
shortening curves for the columns with 12 No. 8 bars are given in Fig. 10. 


Type of Failure 


All columns failed by crushing of the concrete followed by buckling of the 
longitudinal reinforcement. Failure took place in the upper half of all 
columns except three. The typical nature of such a failure is shown in 
Fig. 9(a). The strength of the concrete placed in the lower half of some 
columns was reduced to explore the phenomenon of top failure, which has also 
been observed in numerous previous tests.(9) Though the cylinder strength 
of the bottom batch for Columns 12-8B-1, 12-8B-2, 12-8B-3, 12-6S, 12-6B-1 
and 12-6B-2, was 4 to 14 per cent less than that of the top batch, failure took 
place in the upper half of the columns. For Column 12-6B-3, cylinder 
strength of the bottom batch was 22 per cent below that of the top batch, and 
in this case failure took place in the lower half of the column as shown in 
Fig. 9(c). Even so, the measured ultimate load exceeded by 24 per cent the 
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value calculated on the basis of the cylinder strength for the bottom batch. 

These findings confirm the previously reported observation that the 
column strength of the concrete placed in the lower half of columns is in- 
creased, probably by the improved compaction afforded when the upper half 
is cast. Similarly, the column strength of the concrete placed in the upper 
half may decrease somewhat by water gain from below. To evaluate effects 
of bundling, therefore, measured ultimate loads were compared to calculated 
values based on the average cylinder strength for the top and bottom batch of 
each column. 

The two spliced columns 12-8B-4 and 12-6B-4, which had 1/4 in. clear 
between bars at the splice, failed in the splice region at mid-height as shown 
in Fig. 9(b). 


Effect of Bundling 


The ratios between measured and calculated ultimate load given in Table 2 
exceed one for all except two columns. This confirms previous findings that 
led to the ACI Column Investigation equation. 

The load ratios of Columns 12-6B-1 and 12-6B-2 were 1.03 and 1.04 as 
compared to a value of 1.04 for Column 12-6S that had conventionally spaced 
bars. For 3.67 per cent longitudinal reinforcement, therefore, no detrimental 
effect of bundling was found regardless of tie spacing. 

For the columns with 6.58 per cent reinforcement, which exceeds the maxi- 
mum value of 4 per cent given by ACI 318-56, 1104 (a), the 12-in. tie spacing 
of Column 12-8B-1 led to a load ratio of 0.94 as compared to 1.09 for Column 
12-8S with spaced bars. By reducing the tie spacing to 6 in. for Column 
12-8B-2, the load ratio increased to 1.00. Furthermore, Column 12-8B-3, 
which had a 6-in. tie Spacing and failed above the splice, had a load ratio of 
1.04. Even for 6.58 per cent reinforcement, therefore, no detrimental effect 
of bundling was found when the tie spacing for bundled bars was reduced to 6 
in., that is, 24 tie diameters, or one-half of the least dimension of the column. 

Bundling did not Significantly affect the relationship between applied load 
and column shortening. This is shown for the columns with No. 8 bars in 


Fig. 10. The results for the columns with No. 6 bars indicated a similar lack 
of effect of bundling. 


Splicing 


Bundled reinforcement placed in the corners of a column section may be 
spliced in the same manner as Single corner bars. The bars from below may 
be offset to a position inside the bars above the splice, and a proper amount 
of lap may then be provided. The bars may also be butted and welded. 

In these tests, a splice particularly suitable for bundled bars was explored. 
As shown in Fig. 8, the splices of the three bundled corner bars were stag- 
gered a distance of five bar diameters, and a fourth splice bar 35 bar diame- 
ters long was added at each corner. For Columns 12-8B-3 the bars were cut 
by a saw, and each bar was touching its longitudinal extension. The contact 
was not perfect; after testing, a mortar layer 1/32 to 1/16 in. thick was 
found between the bars. Both of these columns failed outside the splice, and 
the measured ultimate loads exceeded the calculated values by 4 and 16 per 
cent, respectively. 

To simulate less accurate manufacture of reinforcement, the bars of 
Columns 12-8B-4 and 12-6B-4 were cut at the splice by a hydraulic bar cutter 
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with 60-degree cutting edges. The bar ends were wedge-shaped by the cutting 
to an angle of about 90 degrees, and a clear distance of 1/4 in. was provided 
between the bars when the reinforcing cages were tied. Both columns failed 
at mid-height in the splice. However, in spite of the unfavorable conditions 
for a direct stress transfer between bars in the longitudinal direction, the 

No. 6 bar column developed an ultimate strength 5 per cent over the calcu- 
lated value. The ultimate strength of the No. 8-bar column was only 6 per 
cent below the calculated value. As shown for the No. 8-bar columns in 

Fig. 10, the splice did not significantly change the relationship between load 
and column shortening. 

It was planned in subsequent tests to strengthen the splice by longitudinal 
welds between the four bars at each splice. Even without welds, however, 
three out of four spliced columns developed an ultimate strength in excess of 
the calculated values, and it is obvious that the strength of splices with longi- 
tudinal welds would exceed that of the three bars outside the weld. Therefore, 
no columns with welded splices were made. 

It is believed that the short lap used in the splices did not suffice to trans- 
fer stress by bond. The mortar between meeting bars was probably sub- 
jected to a triaxial state of stress so that a compressive strength far in ex- 
cess of the cylinder strength could be developed. To assure that the 
longitudinal bars do not buckle in the splice, the reduced tie spacing used in 
the tests, 24 tie diameters or one-half the column dimension, may be neces- 
sary. If a splice of the type studied is used in eccentrically loaded columns 
so that tensile stress may be developed in the longitudinal bars, the mortar 
between bars cannot be expected to transfer stress and welds, or a longer 
lap, are obviously necessary. 


CONCLUDING REMARKS 


The test results reported confirm the previous findings of Mr. H. M. 
Hadley and others that the use of bundled reinforcement is a sound detailing 
procedure. 

It can be expected that bundling of tension reinforcement in beams will not 
lead to detrimental consequences as compared to spaced bars, if: (1) there 
are not more than four touching bars in each bundle, (2) bond stress computed 
on the basis of external bar perimeter is limited to the values now permitted 
for spaced bars, (3) each bar in a bundle is a deformed bar and is individually 
well anchored, and (4) stirrup reinforcement is provided in regions of high 
bond stress. 

Bundling of compression reinforcement in tied columns can also be used 
even for high ratios of longitudinal reinforcement, if the provisions of ACI 
318-56 regarding other details are strictly complied with. For large amounts 
of longitudinal bundled reinforcement, it is advisable to reduce the maximum 
tie spacing to about one half of that given by the ACI Building Code. 

Since bundling of reinforcement was found to be a safe detail in tests 
involving the extreme cases of bending alone and compression alone, 


bundling should also be a satisfactory detail for members subject to com- 
bined bending and axial load. 


-- 
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MATRIC FORMULATION OF SLOPE-DEFLECTION EQUATIONS 
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SUMMARY 


A simple and systematic way is outlined so that the slope-deflection 
equations may be put in matric form, With the use of the electronic computer, 
the method suggested here will expedite the exact analysis of many complex 

rigid frames under many loading conditions, including sidesway. 


NOTATION 


The important notations used in this paper are listed below. 


X = a column matrix showing values of all joint rotations and all sidesway 
displace ments 


$= a column matrix showing values of elastic rotations at ends of all 
members due to end moments 


M = a column matrix showing values of moments acting on ends of all 
members, exclusive of fixed-end moments 


P = a column matrix showing values of unbalanced moments at all joints 
and of all unbalanced sidesway forces 


A = the statics matrix, expressing the joint (moment) conditions and the 
sidesway (shear) conditions in terms of the end moments 


B = the geometry matrix, expressing the elastic rotations at ends of all 
members in terms of the joint rotations and sidesway displacements 


S = the stiffness matrix, expressing the end moments in terms of the elastic 
end rotations 


Note: Discussion open until March 1, 1959. To extend the closing date one month, a 
written request must be filed with the Executive Secretary, ASCE. Paper 1819 is 
part of the copyrighted Journal of the Structural Division, Proceedings of the 
American Society of Civil Engineers, Vol. 84, No. ST 6, October, 1958. 


1. Prof, of Architectural Eng., Univ. of Illinois, Urbana, Ill. 
1819-1 
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K = ASB = asal = the spring matrix 


INTRODUCTION 


The slope-deflection method of analyzing statically indeterminate beams 
and rigid frames has been well established. It involves the solution of a sys- 
tem of simultaneous equations, equal in number to the unknown joint rotations 
and sidesway displacements. The moment-distribution method, in general, 
consists in solving the simultaneous equations in the slope-deflection method 
by successive approximations. Equilibrium of joints and balance of sidesway 
may be carried in alternate cycles in the same distribution or by combining 
results from separate distribution tables, 

With the recent advance in electronic computation, a matric method is out- 
lined here so that complex rigid frames, such as continuous gable frames or 
tall building frames, may be analyzed quickly for many different loading con- 
ditions including sidesway. The fundamentals of matrix algebra, and its 
application to the analysis of continuous beams, have been treated by 
Benscoter.2 As electronic computing machines are becoming more widely 
available and used, they may be, in ordinary cases, asked only to produce the 
inverse of a square matrix, of which the number of rows or columns is equal 
to that of the unknown joint rotations and sidesway displacements. 


Derivation 


Consider a statically indeterminate beam or rigid frame. Let m be the 
total number of members, j the total number of unknown joint rotations, and s 
the total number of unknown sidesway displacements. Let P be a column 
matrix of (j + s) rows showing values of unbalanced moments at the joints and 
unbalanced sidesway forces, due to loads acting on the structure and settle- 
ments of supports. Let M be a column matrix of 2m rows showing values of 
moments acting on ends of all members, exclusive of fixed-end moments due 
to loads and settlements of supports, Define the statics matrix A of (j + s) 
rows and 2m columns as conditions expressing the balancing moments at the 
joints and balancing sidesway forces in terms of the end moments. Then, by 


definition, 


Let X be a column matrix of (j + s) rows showing values of all joint 
rotations and all sidesway displacements. Let? be a column matrix of 2m 
rows showing values of elastic rotations at the ends of all members, as 
caused by the end moments M, exclusive of fixed-end moments, Define the 
geometry matrix B of 2m rows and (j + s) columns as conditions expressing 
values of elastic end rotations ¢ in terms of values of joint rotations or 
sidesway displacements X. Then, by definition, 


2. Stanley U. Benscoter: Matric Analysis of Continuous Beams, Tran. ASCE, 
V. 112, 1947, P. 1109, 
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Consider a member jk, the j-end of which is connected to the ith joint. The 
end moment M., which acts clockwise on the member jk and counterclockwise 
on 1 ith joint, will balance an externally applied clockwise moment 

= MjAij = Mj at joint i. Thus Ajj = +1. Geometrically, a clockwise rotation 
a Xj of joint i will cause a clockwise rotation of oj = = X; iBji = = X; at end j of 
member jk. Thus Bjj = +1, and Ajj = Bj; 

Next consider the effect of the end shomont M; on sidesway. The free body 
for sidesway is usually a joint or a horizontal member and the force acting on 
this free body and resulting from Mj may be determined by considering 
member jk or a group of members like jkl, etc. as a free body. If, as an 
example, the positive direction of the externally applied force Pj on the ith 
free body for sidesway is horizontal to the right, the positive direction of the 
balancing force MjAij resulting from M; should be to the left on the ith free 
body for sidesway, but is again to the right on the free body for member jk or 
the group of members like jkl, etc. As the member jk or the group of 
members like jkl, etc. is given a rigid-body motion of Xj horizontal to the right 
at the point where the force MjAjj; acts, the counterclockwise rotation at end 
j of member jk is defined to be $j = X{B,, (A counterclockwise rotation of the 
axis of a member will add to the clockwise elastic rotation), During this 
rigid-body motion, the positive work done by the force Mj Aj; to the right in 
going through Xj to the right must be numerically equal to the negative work 
done by the clockwise moment Mj in going through the counterclockwise 
rotation $j = XiBjj- Or, 


5) = (My) 


from which 


By virtue of both considerations discussed above, it is seen that 


(3) 


Consider a member AB with variable cross section, as shown in Fig. la, 
subjected to transverse loads P; and P2, and with end B displaced to B' such 
that angle BAB' = R. The member shown in Fig. la is equivalent, both 
statically and geometrically, to the superposition of those shown in Fig. 
lbed. 

In Fig. 2, define S, 4 and Spp as the stiffness factors at A and B, and Cap 


and Cpa the carry-over factors from A to B and from B to A, respectively. 
Then, referring to Fig. 1d, 


Map Saad, + 


(4) 
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in which 


San = = (5) 
AB Spa * = 


Also, referring to Fig. Ic, 


* (SpptSpa) 


For prismatic members, 


= = )EI/L 


(7) 


The stiffness matrix S is a square matrix of 2m rows and 2m columns in 
which the end moments M are expressed in terms of the elastic end rotations 
@. Entries in the stiffness matrix are those as shown in Eqs. (4) and (5) for 
members with variable cross section and in Eq. (7) for prismatic members. 


Then, by definition, 
(s] 2mx2m {5} 


Substituting (2) and (3) into (8), 
= salx 


Substituting (9) into (1), 


P = 


Defining the spring matrix K = ASAT, 


P = KX 


x = Kp (12) 


Procedure 


The procedure for analyzing a statically indeterminate beam or rigid frame 
by the matric slope-deflection method may be summarized as follows. 

1. Draw the P-X diagram. Number the unknown joint rotations (j) and the 
unknown sidesway displacements (s) consecutively from 1 to (j +s). Draw 
arrows (clockwise for joint rotations) on a diagram of the unloaded structure 
and label them Pj-Xj, Po-Xo, etc. 


3. For full treatment of the statements in this paragraph, see C. K. Wang, 


Statically Indeterminate Structures, McGraw-Hill Book Co., 1953, Arts. 
44 and 54. 


(10) 
(11) 
and 
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2. Draw the M-@¢ diagram. Number the unknown end moments (exclusive 
of fixed-end moments due to loads or settlements) from 1 to 2m. Draw clock- 
wise arrows to act on ends of members and label them Mj- 41, Mg- 9, etc. 

3. Draw equilibrium diagrams of all joints and all sidesway conditions. 
Number these free-body diagrams from 1 to (j + s). 

4. Construct the A matrix of (j + s) rows and 2m columns, expressing P 
in terms of M by observing the equilibrium of the free bodies drawn in step 3. 

5. Construct the B matrix of 2m rows and (j + s) columns, expressing ¢ 
in terms of X by geometric considerations. 

6. Check and see that B = AT. 

7. Construct the S matrix of 2m rows and 2m columns, expressing M in 
terms of 


8. Compute the SAT matrix of 2m rows and (j + s) columns, expressing 
M in terms of X. 
9. Compute the K = ASAT matrix of (j + s) rows and (j + s) columns, ex- 
pressing P in terms of X. Note that the K matrix is symmetric. 
10. Find the inverse of the K matrix. Use the Crout? method on a desk 
calculator or the electronic computer. 
11. Compute all fixed-end moments and unbalanced sidesway forces, due to 


known loading systems and/or known settlements of supports. Construct the 
column matrix P of (j + S) rows. 


12. Find the X matrix from X = K-1P, 

13. Find the M matrix from M = SA? Xx, 

14. Combine the fixed-end moments with the M values to obtain final end 
moments. 

15. Draw the free-body diagram of each member and check the consistency 
of the two final end moments with the values of X pertinent to this member by 
the moment-area method as is ordinarily known to structural engineers. 


Examples 


Three representative examples are shown. The conventional slope- 
deflection and moment-distribution solutions of these examples are contained 
in Wang’s Statically Indeterminate Structures, McGraw-Hill Book Co., 1953. 
Refer to Examples 55, 57, 71, and 73 in this textbook for Example 1, to 
Examples 62, 65, 78, and 81 for Example 2, and to Examples 66, 67, 82, and 
83 for Example 3. 

Example 1. Analyze the continuous beam shown in Fig. 3a due to the loads 
shown in Case 1, Fig. 3b, and due to a 1/2 in. settlement at support B shown 
in Case 2, Fig. 3c. 

The complete solution is shown in Tables 1.1 through 1.9 inclusive and in 
Fig. 3d through h. Note that in this problem j = 3, s = 0, and m = 3. 

Example 2. Analyze the rigid frame shown in Fig. 4a due to the loads 
shown in Case 1, Fig. 4b, and due to a rotational slip of 0.002 rad clockwise 
and a vertical settlement of 0.45 in. at joint D shown in Case 2, Fig. 4c. 

The complete solution is shown in Tables 2.1 through 2.9 inclusive and in 
Fig. 4d through h. In this problem j = 2, s = 1, and m = 3. Note that the A 


4. P.D. Crout, A Short Method for Evaluating Determinants and Soving 


Systems of Linear Equations with Real or Complex Coefficients, Tran. 
AIEE, V. 60, 1941, P. 1235 
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A=L 
Variable I 


(a) 


Sap=CpaSpp 
SBB 


Variable I 
(b) 


FIG. 2 


matrix actually expresses the two joint conditions and the one shear condition, 
The B matrix is determined entirely by geometric considerations; for in- 
stance, when member BC shifts a horizontal distance Xg to the right, the @ 
values at ends 1 and 2 of member 1-2 will be decreased by X./15 and those 

at ends 5 and 6 of member 5-6 decreased by X3/10. It is indeed interesting 
and desirable to check and see that A = BI, 

Example 3. Analyze the gable frame shown in Fig. 5a, due to the loads 
shown in Case 1, Fig. 5b, and in Case 2, Fig. 5c. 

The complete solution is shown in Tables 3.1 through 3.9 inclusive and in 
Fig. 5d through j. In this problem j = 3, s = 2, and m = 4. Note that H and V 
in Fig. 5f are found by solving the two simultaneous equations }Mp = 0 
(BC as a free body) and 5 Mp = 0 (CD as a free body), The Williot construc- 
tion shown in Fig. 5g and h is made in order that the last two columns in the 
B matrix may be found entirely by geometric considerations. Note also that 
the unbalanced H forces on joints B and D are found from simultaneous free- 
body considerations of members BC and CD. In Case 1, H and V are found to 
be +24 and 0 (Fig. 5i), respectively; and they are +1.5 and +0.75 (Fig. 5j), 
respectively, in Case 2. 
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+2/3 +1/3 
1/3 +2/3 


Table 1.5 


K = ASAT (Pekx) 
Multiplier = El, 


+8/3 +5/6 
+5/6 +7/3 +1/3 
+1/3 +2/3 


Table 1.7 


| | case 2| 1 | Case 2° 


+120 +180.9 
2 
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Table 1.2 


sal 
Multiplier = EI, 


(x=x~1p) 
Multiplier = 1/EI, 


70.4262 -0.1639 40.0820 
0.1639 +0.52h6 -0.2623 
#0.0820 -0.2623 41.6310 


Table 1.8 


= 


+71.63 
-85 
+45 63 


Table 1.1 
[lan 
P +1 41 
gi +1 
PS +1 +1 
+1 
+1 
Table 1.3 Table 1.4 
S 
Multiplier - El, 
Ny +1/2 +1/2 
41/2 +1 +1 
5/2 35/3 «5/6 
Table 1.6 
P x 
AL 
P 
P 
X3 #252.0 
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Table 1.9 Table 2.1 


+1 
A -1/15 -1/15 aha -1/10 


Table 2.3 


2/15 
2/15 4/15 


1/5 2/5 


Table 2.) Table 2.5 


saT (m=SATx) K = ASA? (P=KX) 
Multiplier = EI Multiplier = EI, 


+1,/15 -2/75 
-2/75 -3/50 +7/450 


| +#35.82 | 4112.6 
My | +71.63 | #225.3 
-82.38 | -652.0 
41.62 | -252.0 
Table 2.2 
2 
+1 -1/15 | 
1 
+1 2/3 1/3 
41 -1/lo 13 23 
%6 -1/10 2/5 1/5 
6 
| 
| 42/5 -2/75. 
-2/75 
+2/3 -+1/3 
+1/3  +2/3 
+2/5  -3/50 
 -3/50 
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Table 2.6 Table 2.7 


Multiplier = 1/EI, 


#1221356 -0.33487 + 0.78877 
1.28967 + 
+0.78877 +82.611 


“1/20 -1/20 +1/20 -1/20 -1/20 
“1/20 -1/20 +1/20 -1/20 -1/20 


ASCE 1819-17 
: : 2 P,| +10.08 | +520.83 
X) -36 +454.16 
Xo P3 + 7.776] + 20 
Table 2.8 Table 2.9 
Multiplier = 1/£1, 
cave [came 2 | 
My} - 9.06 | - 42.21 
X,| + 30.82] + 95.76 5.01 | + 23.89 
Xp} - 15.58] + 499.32 M #15.09 | +h96.95 
*491.9 | 4061.5 - 0.25 +498 .13 
-35-75 | - 43.96 
6 -32.63 -143.83 
Table 3.1 
My My M3 M Mc; Mg My Mg 
+1 +1 
+1 
P3 +1 +1 
Pi, 
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Table 3.2 
5 


Table 3.3 
Multiplier = EI, 
ty Is By Bg 


0.20000 0.10000 
0.10000 0.20000 
0.35777 0.17888 
0.17888 0.35777 
0.35777 0.17888 
0.17888 0.35777 
0.20000 0.10000 
0.10000 0.20000 


Table 3.4 
sal (M=saTx) 
Multiplier = EI, 
x X3 xX Xs 


+0.10000 -0.015000 
40.20000 -0.015000 
#0.35777 +0.17888 #0.026832 -0.026832 
#0.17888 +0.35777 #0.026832 -0.026832 


#0635777 0.17888 -0.026832 +40.026832 
40.17888  +40.35777 -0.026832 +0.026832 
+0.20000 -0.015000 
+0.10000 -0.015000 


+1 -1/20 
+1 #1/20 -1/20 
+1 #1/20 -1/20 
I’: , 
+1 -1/20 +1/20 
+1 -1/20 7 
-1/20 
BE 
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+0.55777 +0.17888 +0.011832 -0.026832 
+0.17888 +0.71554 +0.17888 

+0.17888 #0.55777 0.026832 +0.011832 
4#0.011832 0.026832 +0.0068664  -0.005366h 
~0.026832 #0.011832 -0.0053664  +0.0068664 


Table 3.6 


(x=x71p) 
Multiplier = 1/£I, 


-~0.78533 + 0.42277 + 12.080 + 19.335 

+1.7902 - 0.78533 - 7.8535 - 7.6535 

+ 0.42277 -0.78533 + 2.7187 + 19.336 + 12.080 
#12 .080 -7.8535 419.336 +54,2.76 4438.07 
#19335 ~-7.8535 +12 .080 +#4,38.07 #542.76 


Table 3.8 Table 3.9 


racer 


+#2802.7 | #5162.7 


ASCE 
Table 3.5 
K = ASA! (Ps=kX) 
Multiplier = EI, 
Py | | 
Po 
4 
x 
X a 
Xi 
Table 3.7 ee 
racers 
+40 -15 M) > 68.64 -81.31 
0 -5 Xi | + 266.0] + 119.0 + 95.2 | -69.41 
Py -2h 710.5 X3 266.0 + 218.7 -102.82 416.78 
| | 2802.7] | 102282 | -21.78 
My 95.2h -33.70 
Mg = 68.64 -55 .57 
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SYNOPSIS 


A general method is presented for calculating the deflection of members 
with variable stiffness to any desired degree of accuracy. The method is 
based on fundamental principles and concepts of engineering mechanics which 
are known to structural engineers. A member of variable stiffness, EI, 
loaded in some manner, may be replaced by selected members of uniform 
stiffnesses, which will have the identical deflection curve of the original mem - 
ber. The loading on these members are changed, but the boundary conditions 
and length are kept the same. The equivalent systems of uniform stiffness 
are derived and proved mathematically. To facilitate the mathematical com- 
putations, an accurate approximation is suggested by the authors which yields 
results to any desired degree of accuracy. 

This method may be applied to any statically determinate or indeterminate 
problem which deals with members of variable stiffnesses. 


INTRODUC TION 


The problem of calculating the deflection of structures with members of 
variable stiffnesses and acted upon by static or dynamic loads is of great im- 
portance in the analysis and design of structures. In this paper, only static 
loads will be considered. 

Several methods have been suggested in the past for the direct or indirect 
solution of the problems of deflections. Some of the best known are the double 7 
integration, the moment—area, the conjugate beam, the three-moment theo- 
rem, the Hardy Cross method of moment distribution, the slope deflectior and 
the energy methods. These methods, when applied to members of variable 
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stiffness, El, become very complicated and difficult to solve. However, it is 
possible to convert a member of variable stiffness to a member of uniform 
stiffness utilizing the theory of the equivalent systems. When this is accom - 
plished, any of the above methods may be readily applied for the solution of 
the problem. In this manner the computational work and the complexity of 
the problem can be greatly simplified. 

The purpose of this paper is to give first the mathematical derivation of 
the equivalent systems and then simplify its application to the engineering 
problems. The conjugate-beam method will be used in the mathematical 
proof of the theory retaining its original assumptions. In addition, it is as- 
sumed that the line joining the centroids of adjacent cross-sections of mem- 
bers of variable stiffness is straight and that all loads are applied in the plane 
of bending about a principal axis. If this line is curved, but the ratio of the 
section depth to the radius of curvature is small, straight beam theory will 
hold with good accuracy. 

Numerical examples will be presented in the later pages to illustrate both 
the mathematical derivation of the equivalent systems and the technique 
which greatly simplifies the application of the method. 


Derivation of the Equivalent Systems 


Consider a simply-supported beam, MN, acted upon by distributed loads 
which yield the deflection curve shown in figure l-a. The moment diagram 
divided by the stiffness, EI, is shown in figure 1-b. The area of the a 
diagram between the points M and N is designated by the letter A and the one 
between the point M and any section n, at distance, u, from the left support is 
designated by A,. X represents the distance between the center of gravity of 


the “ diagram and the right support N. 


The “Second Area Moment Proposition” states that the total displacement 
by of point N from the tangent MQ will be equal to the static moment of the 


a diagram between the points M and N taken about N. Referring to figure 
1-a, note that: 


L 
Ax = (L-x) 


and 


L 
My dx 
(t-*) 
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The actual slope, @, at any section, 1), is 


Mx dx 
J Er fs dx 
Oy-$= 4% - An 
Considering the = “ 


EI diagram as the loading on a beam of span L with sup- 
ports at M andN, a fictitious beam is derived which is referred to as the con- 


jugate beam. The left reaction of this beam is equal to AX or Ou and the 


slope, 6, at any section, n, is equal to the reaction minus the load A,, which 
is the shear in the conjugate beam at that section. 
The displacement y at n, figure 1-a, is equal to 6, - 59. But 


Mxdx _ 
4 = (Ou)(v) = u = 


U 
(u-x) 


Therefore, the displacement y at any section n may then be written as: 


For a beam with variable stiffness, El, the same expression for the de- 
flection curve may be used if E,lI,, the stiffness at any point, is expressed as 
a function of some reference stiffness value, E41. That is, 


Eyl = f (a+ bx") Ei]; (2) 


Substituting Eq. (2) into Eq. (1), the expression for the deflection, y, at 
any section, n, for the variable stiffness member may be written as: 
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dx (L-x) 


L 


(Mydx (u-x) 


If a beam of constant stiffness, E11, is selected so that the length and 
boundary conditions are identical to that of the original beam, but with load- 
ing which yields moment Mey, at any point x, the expression for the deflection 
at any section n, may then be written as: 


L 
Mexdx (L- x) 


U 


The deflection yj, of the variable stiffness beam will be identical to yg, the 
deflection of the equivalent beam of uniform stiffness Ej, if: 


Mx 
Me, f(a +bx") 


(4) 


At any point x. 


The moment diagram of the equivalent system may be determined by divid- 
ing the moment of the actual beam at each point by the value of the function, 
{ (a + bx"), of the stiffness at each corresponding point. 

The loading on the equivalent beam of constant stiffness can be found utiliz- 


ing the fundamental relations, <M =V and SY = W, Stated in ail texts of 


mechanics of materials. 
There is an infinite number of equivalent systems for any one variable 
stiffness beam, since the reference stiffness value selected was arbitrary. 


Application of the Theory of the Equivalent Systems 


When the equations for the moment and stiffness curves of a loaded mem - 
ber of variable cross-section are known, the exact equivalent systems may 
be obtained as prescribed in preceding pages. If the above curves cannot be 
expressed mathematically as a continuous function of x. approximations are 
necessary for the solution of the problem. The resulting accuracy in this 
case depends upon the nature of the approximation. 

It is known from the theory of mechanics of materials that the moment of 
a member at any point is a function of the load. If the moment and stiffness 


function diagrams are accurately plotted, then the t Te bay , which is the 


| 
— 
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moment diagram of the equivalent beam of uniform stiffness, E,l1, may be 
accurately plotted. The shape may be closely approximated by straight lines 
of length As. The number of As increments for an excellent accuracy depends 


M 
upon the shape of the f(a + bxi) curve. Usually three to five increments, 


judiciously selected, are sufficient. 

The equivalent system will then be a constant stiffness member of length 
and boundary conditions the same as the original system, acted upon by three 
or four concentrated loads determined from the equivalent shear force dia- 
gram in the usual manner. 

This method of determining an equivalent system has the advantage of be- 
ing applicable to any member of variable stiffness, whether or not the varia- 
tion of stiffness and moment can be expressed as a continuous function of x. 


M 
It is advisable to approximate the f(a + bxh) diagram quite closely in the re- 


gion corresponding to the least stiff portions of the original system of vari- 
able stiffness. 

When the equivalent system has been obtained, any of the known methods 
in mechanics may be used to determine the deflection at required points. If 
the problem is statically indeterminate, the redundants are first determined 
utilizing the theory of the equivalent systems. A new equivalent system may 
then be obtained as explained in preceding pages, which will be used to de- 
termine the deflection at required points. 

To demonstrate the application of the theory of the equivalent systems, a 
cantilever beam loaded as shown in figure 2-a will be considered. The stiff- 
ness function, f (a + bx"), and the moment diagram of the actual beam are 


M 
plotted in figures 2-b and 2-c, respectively. The f(a + bxn) diagram, which 


is the moment diagram of the equivalently loaded beam, is shown in figure 
3-a. The dashed line on this diagram represents the approximation of the 


i (a+ bxn) curve. The equivalent shear force diagram is presented in 


figure 3-b. For example, the shear at point B is equal to the difference in 
moment between the points C and B divided by the length of Axg. The equiva- 


lent system of uniform stiffness, E41; acted upon by four concentrated loads, 
is shown in figure 3-c. 


EXAMPLE 1: 


/. cantilever beam with rectangular cross-section of width B and depth d, 
given by the expression dy = (413% ys d, is loaded with a uniformly 


distributed load w lbs. per inch. The modulus of elasticity is constant and 
equal to Ey. The variation of the moment of inertia along the length of the 
beam is linear. Determine the exact equivalent system of uniform cross- 
section which will have a constant stiffness E41, and the length and boundary 


conditions of the original beam. Neglect the weight of the beam. For sketch, 
see figures 4-a and 4-b. 


8 
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Solution: 
The moment of inertia at any point x of the beam is: 


3x 
a 


from which the function 


F (e - 
L 


The moment of the beam at any point x is: 


The moment, Mex, at any point x of the equivalent system of uniform 
cross-section and constant stiffness E4lj is: 


6x 


The shear force equation of the equivalent system is: 


ax 


«<b 


The distributed load, w,_, at any point x of the equivalent system is: 


ex 


(au, (9) 


6x 


Equations 7, 8 and 9 are plotted in figures 4-c, 5-a and 5-b, respectively. 
The load, w, and the length, L, are taken as equal to one pound per inch and 
60 inches, respectively. 

The equivalent system of figure 5-b, of uniform cross-section and of con- 
stant stiffness E4l, will produce the identical deflection curve as the original 
beam of variable stiffness shown in figure 4-a. 


EXAMPLE 2: 


In example 1, determine an accurate equivalent system of uniform cross- 
section and constant stiffness E41; by approximating the shape of the 


f (a+ Curve. 


(7) 
ax x 
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Solution: 


The approximation of the Mg curve with three straight lines, 4s,, As2 and 
Ass is shown in figure 6-a. The equivalent shear force diagram is plotted by 
dividing the variation of moment by the corresponding increments of length 

Ax , Axg and Axg and it is shown in figure 6-b. The equivalent system of uni- 
form cross-section with constant stiffness, E,I, is shown in figure 6-c. 


CONC LUSIONS 


From the foregoing, we may conclude that the theory of the equivalent 
systems, derived from the fundamental principles and concepts of mechanics, 
may be applied to structures with variable stiffness members. This method 
provides a relatively simple procedure for determining the deflections of 
complicated structures without requiring prepared tables, graphs or charts. 
A slide rule, a firm understanding of the principles of mechanics and the en- 
gineer in industry may readily apply this theory to the most complicated 
problems. The method is very general and it may be extended for the solu- 
tion of staticaly indeterminate problems, such as continuous structures of 
variable stiffness, trusses, plates and shells without complicating the pro- 
cedure. It may also be used for the determination of the natural frequencies 
of vibration of the above structures. 


It should be emphasized that the approximation with straight lines of the 
Me diagram, greatly reduces the complexity of the solution and yields results 
to any desired degree at accuracy. 
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SUMMARY 


This paper pertains to certain problems dealing with the amplification of 
stress and displacement in guyed towers when changes in geometry are in- 
cluded in the analysis. A simplified method for determining the slack stress 


in the guy and horizontal displacements due to wind load has also been de- 
veloped and presented. 


SYNOPSIS 


A simplified method for the determination of the slack stress in the guys 
of a tower subjected to wind load has been presented. Equations for the 
determination of displacements at guy levels have been developed in terms 
of the wind load, initial guy tension, and slack stress. 

A simplified method which is useful in solving problems pertaining to the 
buckling or amplification of stress and displacement in guyed towers has been 
developed. Only fundamental concepts are used; these are extended to include 
the displacement of the tower axis due to external loads and internal distor- 
tions. The increase, or amplification, of the computed stress and displace- 
ment or bending moment at the center of a loaded section of the tower (when 
the displacements of the axis are included in the analysis) is developed and 
presented in graphic form. 

The amplification charts are useful in showing the entire range of struc- 
tural action for members subjected to axial compression and transverse 
bending. These charts are divided into (A) the ordinary range, where the 
usual methods of structural mechanics apply with little error; (B) the flexible 
range, where the ordinary methods of analysis must be extended to include the 
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effects of displacement; and (C) the dangerous range, where the amplification 
of stress is so large that there is a possibility of structural collapse. 

Charts of this type are valuable because they present graphically (a) the 
range of structural action for a given tower, (b) the magnitude of the amplifi- 
cation factor when displacements of the axis are included in the analysis, and 
(c) possible modifications that may be made to reduce the amplification. 


INTRODUCTION 


The guyed tower is a structural form which has become increasingly im- 
portant because of the need for high television and radio towers. A typical 
guyed tower is shown in Figure 1. This tower is an example of a typical tele- 
vision tower, and has the following properties: 


Tower: 
Height = 1000 feet 
Cross-section = 10' Triangular 
Legs = 5'' Round, Tl-Steel 
Girts and Diagonals = A7-Steel 
Total Weight = 440,000 pounds 


Bridge Strands 
For Single Guys For Double Guys 
Figure ld Figure le 


Initial Tension 
Size Size per Guy 


G3 1-3/8"' Dia. 2-1" Dia. 22.8 kips 
G2 1-9/16"' Dia. 2 - 1-1/8"! Dia. 28.7 kips 
Gl 1-5/8" Dia. 2 - 1-3/16'' Dia. 31.7 kips 


Considerable time is expended on the preliminary design and analysis of 
tall steel towers because of the lack of a simplified method for performing 
these operations. It is not permissible to use the ordinary methods of struc- 
tural mechanics for flexible steel towers, because the displacements (which 
are normally neglected) may increase or amplify the stresses which, in turn 
continue to increase the displacements. It should be realized that any refine - 
ments which are introduced into an analysis are included for reasons of safety 
rather than economy. 

In order to determine the stresses in a guyed tower, it is first necessary 
to assume a preliminary design. In general, the effect of variations can only 
be determined by analyzing several designs. Because of the great amount of 
time and labor required in this process, elaborate analysis is seldom practi- 
cal. A device useful in determining the influence of displacements on the 
preliminary design is termed an “amplification chart,” whereby the struc - 
tural action (as well as the effects of varying the tower proportions) can be 
visualized easily. This amplification chart can be modified for conditions 
other than those originally specified. Thus, one chart can serve for the 
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various tower types, proportions, materials, loading, and for various com- 
puted stresses. The amplification chart is developed from simple structural 
concepts. The results obtained by use of the chart are extended, by succes- 
sive correction, to include the displacement of the tower axis. These results 
agree closely with the results arrived at by the formal methods of structural 
mechanics. 


Horizontal Displacements Due to Wind 


The horizontal displacements of a tower depend upon the wind direction, 
slack stress, and initial tension in the guys. 

For the wind direction shown in Figure 2 (Wind I), the horizontal displace - 
ments of the tower axis at a guy level may be computed approximately by 
Equation 1. 


¢ + Pg - Py) L 


(1) 


Sina A, Ey 


where 4H, = Horizontal Displacement for Wind I 


Pin Wind Force in a Guy 


For single guys P = ——— (2) 


For double guys (3) 


Ps = Slack force in a guy 
P; = Initial force in a guy 


Wind I 
Pg 


Triangular Tower Section 
Single Guys 


P 


Pg 
Ps 
Double Guys 


Figure 2. Sketch Showing Wind I and the Resulting 
Forces in a Triangular Tower for Single 
and Double Guys 


= 
| 
L 
PittPs 
Pg 
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H = Total horizontal wind reaction at guy junction 
A, = Area of an individual guy 


Thus, before it is possible to determine horizontal displacements, it is 
necessary to determine the slack force Ps. The slack force depends upon 
the wind reaction, H, and the initial force, Pj, in the guy. The relation be- 
tween the various forces may be obtained from statics. 

Second order differences resulting from changes in guy geometry are 
negligible and for convenience are not included in the following development. 


From =0 
Then (Pry + 3 Pg - 3 Pi) Cosas 
But cos a # 0 


Therefore (Py + 3 Po - 3 =O 


and for a selected factor of safety, F. S. = 2.5, the design load, P,, may be 
computed as follows: 


Bw Pit Ps = BeSe = o.4 BeSe (4) 
FeSe 


where B.S. is the breaking strength of the cable. 
The slack force would be 


Pg, 165 Py - Oo2 = Py 1/3 Py (5) 


Similarly, for the wind direction (Wind II) shown in Figure 3, the horizon- 
tal displacement may be computed by Equation 6. 


P L+P 
Pg 
P 
Pg Wind IT LL+ Ps 
Single Guys 
PL 
Pg 
Pg 


Double Guys 


Figure 3. Sketch Showing Wind II and the Resulting 
Forces in a Triangular Tower For Single 
and Double Guys. 
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sn. 2 (6) 


The corresponding slack force for Wind II would be 


#3 0.8 = Py - 2/3 Pll, (7) 


Thus the slack stress for Wind I is not the same as for Wind II even 


though the wind load, P} 1» and the initial tensions, Pj, are equal in both 
cases. 


For example, consider an arbitrarily chosen value for the Initial Tension 
Pi 0.2 B.S. 
Then the corresponding value of the slack force for Wind I would be 


Ps = 0.1 B.S. and for Wind II, Pgy = 0 


The allowable value of the wind load would be 


= 0.3 B.S. 


The corresponding magnitude of the horizontal displacement, AH, for both 


Wind I and Wind II, as computed by Equations 1 and 6, would be equal for this 
case as follows: 


Au = Oco2 BeSe lg e PLL Lg 


Sina Ag Eg 3 sina Ag (8) 


Upon substituting Equation 2, in 8, the horizontal displacement for single 
guys would be 


2 4H Lg 
3 Sin* a Ag Eg (9) 
j 
And by substituting Equation 3 in Equation 8, the horizontal displacement 
for double guys would be 


1 H 
we (10) 
9 “9 

The displacements computed by Equations 9 and 10 are the same when the 
total cross sectional area for double guys is twice the area for single guys. 

However, this is a limiting condition for Wind Direction II as the slack 
force would be equal to zero. It is thus recommended that the guys be de- 
signed so that the slack force is slightly greater than zero for both wind con- 
ditions. 

The magnitude of the various guy forces as a function of the breaking 
strength (B.S.) of the guy for conditions, other than the example above, may 
be conveniently obtained from Figure 4. 

The chart of Figure 4 is entered with a design value of Pi in terms of the 


breaking strength and then vertically up to the corresponding value of P, and 
Py for either Wind I or Wind II. 


| 
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Pwt =PrL+ Pg 0.4 B.S. 


P Wind I & 
LL Wind II 


Figure 4. Chart Showing the Relationship between the Live- 
load Force, the Initial Tension and the Slack Stress in the 
Guys of a Guyed Tower. Secondary Effects not Included, 


0.4 
0.3 
| 
B.8. 
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PLL 
B.S. 
0.2 
0.1 
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Pi Pert (Wind IT) 
Pp, - 
Pu 
0 
0.1 ae 0.3 0.4 
B.S. 
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For a factor of safety of 2.5, Figure 4 shows that the larger working 
strength P,, results from Wind I and the smaller Slack Stress (Psp) results 


from Wind Il. The slack stress is zero for Wind II when the initial tension 
Pj = 0.2 B.S. The initial force (Pi) should probably be slightly greater than 
0.2 of the breaking strength. 

Knowing the value of Py,},, Ps, and Pj from Figure 4, the magnitude of 
the horizontal displacement corresponding to Wind I or Wind II may be deter- 


mined for towers with single guys from Equation 9 and with double guys from 
Equation 10. 


Horizontal Displacements Due to 
Guy Lengthening or Shortening 


Guyed towers are sensitive to changes in guy tension. Extreme caution 
must be used to insure simultaneous lengthening or shortening of all guys 
that meet at a particular tower level. 

Small changes in tension of a single guy will change the geometry of the 
tower axis and may cause buckling because the small horizontal displace - 
ments are amplified by the large axial forces. 

Consider the typical guyed tower shown in Figure 1. 

From fundamental relations between statics and geometry and the way in 
which they are related through the properties of the materials, the horizontal 
displacement may be determined from Equation 11, which is derived for guy 
lengthening or guy shortening, 4,, at a given guy level: corrections due to 
changes in sag are of a secondary nature and negligible for properly designed 
towers. 

For convenience, the cross-sectional area of the tower is assumed 
constant. Therefore: 


2 3. Ag Eg 
AH = Ae — (1 + 3 Cos~a he Ey (11) 


However, the above equation, based on the original geometry, is only ap- 
proximate as changes in tower geometry must be included in the analysis. 
A convenient technique useful to include changes in geometry is to multiply 
computed values by the amplification factor. Thus the computed transverse 
displacement 4H may be corrected for changes in geometry by multiplying 
by the amplification factor ® . 

The total horizontal displacement, 4H;, may then be computed as follows: 


2 
= de 


Ag Eo 
(1 + 3 Cos%a ) (12) 
Ay Ey 


Sina 


The value of @ for displacement may easily be obtained from Figure 5. 


Amplification Chart 


The amplification chart (Figure 5) shows in graphic form the increase or 
amplification, of displacements due to combined axial compression plus the 
transverse bending when the displacements of the axis are included in the 

analysis. This chart is based on a computed total stress, f = f. +fp, of 


ASCE GUYED TOWERS 


1821-9 
20 kips per sq. inch at the center of a tower section between points of inflec- 
tion which supports a uniformly distributed wind load throughout the tower 
height. This selected wind loading is used even though it does not produce 
the maximum bending unit stress, but it is easy to analyze and very nearly 
correct. For a computed stress, f = fy +f}, (in kips per square inch) at the 
center of the tower section between points on inflection which is not equal to 
20 kips per square inch, (that is, for the higher working stresses associated 
with the higher strength steels), the value of 1/r should be multiplied by 
Vf 7 20 

Transverse loads on a guyed tower cause the elastic curve to move away 
from the tower axis. The displacements resulting from transverse bending 
are the moment arms for the dead load and the live load axial forces. Al- 
though the displacements are small, the additional bending moments may be 
large because the live-load plus the dead-load axial forces are large. In 
flexible members the bending stresses and displacements are increased, or 
amplified, by the additional bending moment resulting from the product of the 
large axial forces and the transverse bending displacements. The total hori- 
zontal displacement, oH, , with changes in geometry included is: 


An, = 4H § (13) 


and the amplified total combined stress is 

14 
where 


1 
= 
1. 78 (- (15) 
neE 


where 4H is the horizontal displacement neglecting changes in geometry and 
f, is the axial stress and f, represents the bending stress before amplifica- 
tion and resulting from the dead load plus the live load. The amplification 
factor is represented by ® and may be obtained in Figure 5. In the equation 
for the amplification factor, n is a shape factor that depends upon the loading 
and properties of the member, and L is the slenderness ratio. In Figures 5 


r 
and 6 the value of n is 1/10, and E is equal to 30,000,000 per square inch. 
The amplification chart, Figure 5, shows the amplification factor for bend- 
ing moment or displacement, @ , plotted against the slenderness ratio L/r 
and the ratio of axial stress and bending stress fg/fp or fp/fa. By use of 
both fa/fp and fp/fa, extremes can be shown conveniently on one chart. In 
addition to demonstrating the magnitude of the amplification of bending 
moment or displacements when displacements are included in the analysis, 
this chart also shows arbitrary ranges of structural action. 
For convenience, the amplification factor for total stress $' is obtained 


from Figure 6. The amplification of total stress is obtained from the follow- 
ing equations: 


f, + fb J 


= (16) 
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AMPLIFICATION FACTOR ® 
for 
BENDING MOMENT OR DISPLACEMENT 


— 


B-CMtical Range 


AwOrdinary Range 


066 0.8 1,0 0.8 0.6 0.4 0,2 

fa=fp » fp=0 
b 


Figure 5. Amplification Chart Showing the Amplification 
Factor, , for Bending Moment or Displacement as a Function 
of the Slenderness Ratio as ordinate and the retio of the 

Axial and Bending Stress as Abscissus, 


Use of the Amplification Chart 


The amplification chart can be used (a) to determine the magnitude of the 
amplification factor for bending moment or displacement and stress when 

displacements of the axis are included in the analysis, (b) to show the influ- 
ence of modifications of design on the amplification of displacement, and (c) to 
describe, in a general manner, the structural behavior of guyed towers for 
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250 


AMPLIFICATION FACTOR 
for 
STRESS 


ite 


~~~--1,40- 
ByFlexible Ran¢ 


1.06 
A-Ordinary Rane ¢ 
0.4 O58 5.6 “ 
fa fp 
fy fe 


Figure 6. Amplification Chart for Total Stress 
as a Function of the Slenderness Ratio as Ordinate 
and the Ratio of Axial and Bending Stress as Adscissus 


the entire range of structural action. A similar chart has been developed 
for arches. (1) 

The ratios needed to determine the magnitude of the amplification factor 
at the center of any section between points of inflection of a guyed tower may 
be computed as follows: 


1. The slenderness ratio L/r is computed and then modified for the com- 
puted stress which depends upon the type of steel. Figures 5 and 6 are based 
on a steel with a computed stress of 20 kips per square inch. Higher 


tH 
\\ 
\ \o-Dangerou’ fs fa+fp $ 
+} 
150} \ iy, 
| 2,00 
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computed stresses may be allowed for the higher strength steels. If the com- 
puted stress, f = fy +fp, is other than 20 kips per square inch, the modified 
slenderness ratio (L/r) moq, needed to enter the chart is determined by 
multiplying the slenderness ratio for the stress other than 20 kips per square 
inch by the square root of the quotient of the other stress divided by 20 kips 
per square inch. The lengths between points of inflection can be determined 
approximately by sketching the elastic curve for the loads on the tower. It is 
important to note that the (L/r) values pertain to the bending stress at the 
center of a loaded tower section between points of inflection. These values 


vary with the loading, guy tensioning, guy lengthening, type of steel, and the 
properties of the cross section. 


2. The ratio of bending stress to axial stress, fp/fa, is computed. The 
slenderness ratio and the ratio of bending to axial stress are all that are 
needed to enter the amplification chart. Variations in guy location, depth of 
section, moment of inertia, and cross-sectional area are included when the 
ratio of axial and bending stress is determined. 

The amplification chart is entered from the abscissa with the ratios of 
either fa/fp or fp/fa (whichever is the smaller) and thence vertically up to 
the value of (L/r) mod. The intersection represents the range of structural 
action and the magnitude of the amplification factor for bending moment or 
displacement, Figure 5, and for total stress, Figure 6, at the center of the 
section when displacements of the axis are included in the analysis. 


Example 1. Amplification of stress due to Wind Load. Consider a steel 
tower of proportions similar to those shown in Figure 1. From the propor- 


tions of the tower and the loadings, the ratios needed to enter the chart are 
computed: 


L = 238 feet 
r = I/A = 0.405 (10) = 4.05 ft. 
L/r = 238/4.05 = 58.8 


Therefore f = sro0e ‘ 1,868,000 (8.78) _ 11,100 + 10,800 = 21,900 psi. 


Because the charts are based on a computed stress of 20 kips per square 
inch, it is necessary to modify the slenderness ratio for any other computed 
stress. The computed stress (from Equation 18) is 21.9 kips per square 


inch; (guy eccentricity not included) therefore, the modified slenderness ratio 
would be 


(1204) = 6165 


The ratio of the bending stress to the axial stress is then computed: 


fp 10.5 
fa llel 
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The chart shown in Figure 5 is entered with the ratio of f,)/f, equal to 
0.97 - thence, vertically up to the value of 1/r equal to 61.5 The intersection 
shows the range of structural action and the approximate value of the amplifi- 
cation of the bending moment or displacement when changes in the dimensions 
of the arch axis are included. The amplification of the displacement midway 
between points of inflection interpolated from the chart, is approximately 1.16 
and in the flexible range of structural action. Thus the actual stress is ob- 
tained from Equation 14 and is: 


= + fy = 11,100 + 10,500 (161¢) = 23,600 psi 


Similarly the total stress could also be obtained from Figure 6 where the 
amplification factor for total stress @', interpolation from the chart, is 
1.07. Thus the total stress as computed by Equation 17 is 


f, =f q * = 21,900 €1-07) = 23,500 psi. 


Example 2. Amplification of stress due to guy lengthening or guy shorten- 
ing. 

To illustrate the dangers of guy lengthening or shortening, suppose that 
after the tower has been erected it is discovered that the tension in the lower 
set of guys needs to be increased or decreased by 15% of the buy breaking 
strength. 

For the typical tower shown in Figure 1, with a set of single guys 1-5/8 
inch in diameter, at the first guy level, 15% of the breaking strength of 
300,000 lbs. corresponds to a change in tension of 45,000 lbs. This change 
corresponds to the following guy lengthening or shortening (changes due to 
sag are of a secondary nature and negligible here): 


45,000 x 42h x 12 


3073 ine 


AE 2.05 x 30,000,000 


The corresponding horizontal displacement of the tower at the first guy 
level is computed from Equation 11 as follows: 


Eg 
4H Ae qi + J 3 73 x 2 qa + Oc( 
e O6012) = 10.0 ine 
Sina Ay Et 0-707 
t 


If only one guy is lengthened or shortened instead of all three guys which 
meet at the same level simultaneously, the axis of the tower will change in 
accordance with the geometry shown in Figure 7. Secondary effects are 
small and not included here. 

The length between points of inflection has increased greatly from the 
original length of 238 feet (see Figure 1) to 430 feet (Figure 7) corresponding 
to a slenderness ratio of L = 120. 

The critical buckling stress for a steel tower with the slenderness ratio 

of 120 may be computed from the following formula: 


{, = ————- = 20,600 psie (19) 
(L/r)* 
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Tower Axis. 


H 10.6" 


Displaced Tower Axis b. Bending Noment Disagrag 


Figure 7. Sketch Showing The Displaced Tower Axis 
due to Guy Lengthening of G, 


It is noted that the critical stress, fx, is an ultimate stress with no factor 
of safety included. 


The portion of the bending moment due to a horizontal displacement may 
be obtained from the following formula: 


Mm = c+ SET 4 «1,030,000 1b, ft. 
2 
L 


Where c isa factor that depends upon the end fixity at the supports. For 
the problem under consideration c = 0.6. 

Additional bending moments, Mg, result from the axial forces times the 
relative horizontal tower displacements. The axial force above the section 
under consideration is composed of the tower weight, Pg = 264,000 pounds, 
and initial guy tension, Pj = 218,000 pounds. 


Thus 


= (264,000 + 218,000)(10.6/12) = 425,000 lb. ft. 


The computed bending moment is the sum of the above moments due to 


tower displacement plus the moment resulting from the axial force times 
displacement and is: 


M = Mp + M, (21) 


M = 1,030,000 + 425,000 = 1,455,000 lb. ft. 


The corresponding computed stress with changes in geometry due to loads 
not included, is obtained as follows: 
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+ fy 


¢ = 480,000 155,000 
59 99 


f = 8140 + 8450 = 14590 psi. 


Thus the stress computed in the usual manner is not the total stress in the 
tower as changes in geometry must be included. 

The actual total stress, amplification included, may be obtained by use of 
Equation 14 as follows: 


From Figure 5, for L/r = 120 and f{,/fp = 0296, $ = 1.7 


Thus +t 


and fy, = Blo + 8450 (1.7) = 22,510 psi 


Similarly, the total stress compution by use of Figure 6 would be 
{,=f o = 16,590 (1.36) = 22,560 psi. 


Thus the tower of this example will probably collapse since the computed 
total stress of 22,500 psi exceeds the critical stress of 20,600 psi. 


The corresponding amplified horizontal displacement as computed by 
Equation 13 is 


= AH ° 106 x le? 18.0 ine 
Possible remedies would include simultaneous lengthening or shortening 


of all the guys that meet at a given level and by increasing the number of guy 
levels from three to four or five. 


Influence of Modifications on the Amplification of Stress 


The amplification of stress in guyed towers subjected to combined axial 
compression and transverse bending is influenced by the slenderness ratio 
1/r, the type of steel or computed stress, and the ratio f,/fp- The amplifica- 
tion of stress can be reduced by decreasing the slenderness ratio. Fora 
given length between inflection points, the amplification of stress decreases 
for increasing values of the radius of gyration. 

A higher allowable computed stress is permitted for certain types of steel 
and for wind load. The relative amplification of stress is influenced by varia- 
tions in the allowable or computed stress (Equation 15). An increase in the 
computed stress results in a disproportionate increase in amplification 
particularly for larger values of the slenderness ratio. 


Use of the Amplification Chart to Describe Structural Action 


The amplification chart may be useful in describing the range of the struc - 
tural action of guyed towers. The chart has been divided into three arbitrary 


ranges: (a) the ordinary range, (b) the flexible range, and (c) the dangerous 
range. 
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The ordinary range shows that the analysis, by customary methods of 
structural mechanics, gives sufficiently accurate results. In the ordinary 
range, changes in the dimensions of the tower cross section and changes in 
the displacement of the tower axis have little influence on the amplification 
of stress; therefore, (to that degree) superposition of effects is applicable. 
In the ordinary range, for a given slenderness ratio of less than 70 and for 
values of fa/fp from 0.5 to 2.0, the amplification of stress remains sensible 
constant, irrespective of the loading and proportions of the tower. For de- 
creasing values of f,/fp and fp/fa less than 0.5, the amplification of stress 
decreases for constant values of 1/r. Figure 6 has been drawn for a total 
stress before amplification, of 20 kips per square inch. If the computed 
stress is greater than 20 kips per square inch, it is necessary to consider 
the corresponding disproportionate increase in the amplification. The ampli- 
fication of displacements, Figure 5, for constant L increases with increasing 

values of fg/fp. In the ordinary range, elastic buckling is highly improbable 
except for relatively small values of f/f, - at which values the amplification 
of displacements is relatively large. 

In the flexible range, the ordinary methods of structural mechanics must 
be modified and extended to include the effects of displacements. Superposi- 
tion of effects is not applicable, and it is not permissible to separate live load 
and dead load. The amplification of stress in this range has been arbitrarily 
limited from 10% to 50%. The area of constant amplification of stress for a 
given 1/r - ratio is less pronounced in the flexible range. The amplification 
of stress increases for increasing values of f,/fp up to a certain value and 
then decreases. The amplification of displacements continues to increase 
for increasing values of f,/f}. Variations in loading and in properties of the 
tower have more influence on the amplification of stress in the flexible range 
than in the ordinary range. Amplification of stress is more pronounced for 
small changes in the slenderness ratio. Attention is given to variations in 
the computed stress because increases in the computed stress result in dis- 
proportionate increases in the amplification. Elastic buckling may be possi- 
ble at the relatively larger values of 1/r and for values of f,/fp greater than 
unity. 

In the dangerous range, the amplification of stress is greater than 50%. 

In this range the amplification of stress and displacement is so large that 
there is danger of elastic or plastic buckling. The line separating elastic 
and plastic buckling depends, to some extent, on the type of steel used. The 
higher allowable stresses associated with silicon, nickel, and T-1 steel, re- 
sult in disproportionate increases in amplification. 

For an assumed tower design, the amplification chart shows immediately 


if danger exists, and what modifications are indicated to avoid this danger 
without elaborate computations. 


Amplification of Stress and Displacement 


The bending moment obtained from the wind live load only is not the actual 
bending moment in the structure. For convenience, the moments resulting 
from guy eccentricity have not been included in this discussion. The bending 
moment resulting from the transverse loads is increased (or amplified) by 
additional bending caused by the axial forces. Usually a guyed tower is 


_| 
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designed so that the live load deflections are small. The bending moment due 
to the axial force is computed by multiplying the axial force by the distance 
from the elastic curve to the original tower axis. 

If the displacements of the tower are neglected, the same bending moment 
will be determined for the wind load only as for the wind load plus axial 
forces. Thus, by the usual theory, it is permissible to separate the wind load 
from the guy tensioning and to add the partial effects to obtain the total effect. 

In a guyed tower the displacements of the tower axis, although small, may 
form a large part of the dead-load plus the wind-load bending moment, as 
shown in Figure 8. If the displacements of the tower axis are included, 
larger bending moments will be determined for the wind-load plus the dead- 
load than for the wind-load only. The transverse displacements are moment 
arms for the combined wind-load plus the dead-load axial forces. When the 
displacements of the axis are included, it is not permissible to separate the 
live load from the dead load. Not only are the bending moments and the dis- 
placements increased or amplified by the live-load plus dead-load axial 
forces, but the sums of the partial effects of the live-load and the dead-load 
are not equal to the total effect. However, it is often convenient to permit an 
arbitrary allowable increase, to the results obtained by assuming that the 
structure does not deform, and still be safe. 

The stages of additional bending continue and may be represented by a 
mathematical series which may or may not converge. If the bending con- 
tinues to increase (non-convergent series), the structure will collapse. In 
Figure 8 there is shown a straight member subjected to both axial and 


transverse loads. The equation for the displacement at the center of the 
beam is: 


(22) 


Wind 
Load 


Original 
Tower Axis 
Wind Only 


Figure 8. Sketch Showing Axial Forces and the Transverse 
Wind Load on a Tower Sectior Between Points 
of Inflection, 
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in which M is the ordinate of the bending moment diagram at the center of 


the beam; and n = Po for a simply supported beam with uniformly distributed 


wind live load. For a more complete development of the amplification factor 
than presented here see Reference (1). 

The bending moment for the transverse loads is increased by additional 
bending caused by the axial forces. The shape factor for the transverse 
loads np is not the same as the shape factor for the axial loads n,. For 


ordinary guyed towers, the error resulting from setting ng = nb is small and 
in the order of 3%. 


The additional bending moment (Figure 8) is 


Mag = Pa 4b (23) 
which causes a further displacement 


Na Pa Ab L? 


4aa * 


and an increased bending moment 


2 
Ng Ap L 
Mig Pa (25) 


This process is continued, and the total bending moment can be represented 
by: 
2 
Pa L Pa 
Ap 
EI EI 


M = My + Pg dy + Py Ap (26) 


and letting successive values of n equal n,, 


Pg Ab 


le 


M Mp + 


EI 
Equation 28 can also be written as 


2 


If Ng Pq 
EI (27) 
ny Pa my Ny 
EI EI 
(29) 
n, Pa L 
EI 
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By setting n, equal to np, an approximate expression for the amplification of 
displacements and bending moments is obtained. 


(30) 


and 


Slight errors less than 3% are caused by neglecting successive changes in ng. 


The quantity is the amplification factor @ . 


E I 


The total stress for a straight member or tower section, subjected to both 
axial and transverse loads can be computed from 


(31) 


The effects of temperature change, plastic flow, and movement of the tower 
base are not included in Equation 14. 
This investigation and the amplification formula for a horizontal, straight 


member are applicable also to a curved member with supports on the same 
elevation. ( 


SUMMARY 


Simplified equations for the determination of displacements of guy levels 
in a tower have been developed in terms of the wind load, initial guy tension 
and slack stress. Amplification charts for stress and displacement which 
have been found useful in the predesign of guyed towers have also been pre- 
sented. These charts demonstrate graphically the increase or amplification 
of stress and displacement at the center of a section between points of inflec- 
tion of a guyed tower when displacements of the axis are included in the 
analysis. In the development of the chart only familiar concepts were used. 
The value of the chart is that it shows immediately when elaborate 
methods of analysis are necessary in the design of flexible towers. The 
chart also indicates what modifications may be made, sc that the ordinary 
methods of structural mechanics give adequate results. 
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THE DEAS ISLAND TUNNEL# 


Closure by Per Hall, Troels Brondum-Nielsen, 
and H.R. Kivisild 


PER HALL, | A.M. ASCE, TROELS BRONDUM-NIE LSEN,” and H. R. 
KIVISILD?.—The authors wish to thank Messrs. Rooney and Eremin for 
their comments on the paper. 

As regards cc mparative costs of alternative bridge and tunnel structures, 
the writers wish to mention that the reasons for a tunnel structure being 
competitive are partly those mentioned by Mr. Rooney, viz. severe navigation 
clearance restrictions and poor ground conditions for bridge foundations, and 
partly the following: First the extremely flat, low lying terrain in the area, 
the rigid maximum grade specifications and the clearance requirements 
would necessitate much longer approaches for a bridge than for the tunnel. 
Secondly the price of heavy structural steel in British Columbia is consider- 
ably higher than elsewhere. 

In regard to the ventilation system, the writers should like to comment that 
the system employed in the Deas Island Tunnel is designated as semi- 
transverse. Air is exhausted laterally into the ventilation duct, for the first 
half of the tube and introduced into the roadway for the second half of the tube, 
as is shown in Figure 1. 

Basic ventilation is provided by four fans with a capacity of 1,000,000 cfm, 
which create the main circulation of air in the tunnel. A secondary circula- 
tion is created by the movement of traffic through the tunnel. The fans oper- 
ate in a pre-set, time-controlled pattern with automatic alteration controlled 
by carbon monoxide detectors, visibility meters and fire detectors with a 
manual override for use during peak traffic periods and emergencies. 

No spare fans are installed in the tunnel. When one fan is out of operation 
for any reason a damper at the midpoint of the ventilation duct can be opened 
so that the other fan can serve the full length of the tube. 

As all the fans are reversible it is possible, in case of fire in the roadway 
tunnel, to combine suction and pressure in such a way that the air movement 
above the fire will be in the direction of the normal traffic flow. In this way 
smoke and heat will be moved away from any cars held up by the fire. 

In regard to the question of automatic control, both ventilation and lighting 
are completely automatic, with provision for manual override. Lighting with- 
in the tunnel is controlled by the programming clock, with an overriding con- 
trol by photo cells mounted at the tunnel entrances. Traffic signals are 


. Proc. Paper 1436, November, 1957, by Per Hall, Troels Brondum-Nielsen, 
and H. R. Kivisild 


. Exec. Vice-Pres., Foundation of Canada Eng. Corp., LTD., Montreal, 
Quebec, Canada. 


2. Christiani and Nielsen. 
. Foundation of Canada Eng Corp., LTD., Montreal, Quebec, Canada. 
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operated manually from the control room, where the operator can detect 
emergency conditions visually over a closed circuit T.V. system. The traf- 
fic signals can also be operated manually from stations at the tunnel en- 
trances. In cases of emergency, with excessive carbon monoxide content, 
red lights at the tunnel entrances are automatically activated by the carbon 
monoxide meters to stop traffic entering the tunnel. 

Further control of operation and traffic under all conditions is provided by 
means of telephone and loudspeaker communication systems. These features 
complement the T.V. system to provide the operators with complete visual 
and vocal contact with motorists in case of emergency. During normal day- 
time operation and at night, all installations will operate automatically, each 
action being registered on chart recorders, while during peak hours the con- 
trol will be manual. A small staff will be employed on the tunnel to take care 
of inspection, manual operation and maintenance. 

Mr. Eremin’s comments on structural rigidity and protective measures 
are interesting. In regard to the protective armour, the writers would like 
to add that the blanket consists of rock ranging in size from 1500 lb. on top 
of the tunnel to 500 lb. on the river bed. This rock is placed on reinforced 
concrete mattresses which act as a filter against the underlying sand and as 
a reinforcement for the rock blanket during periods of heavy scour action. 


The settlements observed during placing of the tunnel elements were as 
follows: 


Element No. North End South End 
1 9.0 m.m. 11.5 m.m. 
2 1.0 9.1 
3 4.8 15.5 
4 2.4 6.8 
5 Not Available 
6 0.0 7.2 


These figures include rebound of the temporary foundation blocks as well 
as the settlement of the elements themselves. The rebound is estimated to 
be about 1.0 m.m. 

The reason for the higher values for South End movement is that in all 
cases the North End was the one connected to the previously placed element, 
the South End being more or less free. An average figure could be obtained 
by interpolation between the two figures shown. 

Before sinking the first tunnel element the sinking procedure was 
changed so that Fig. 12.2(d) of the paper is not up-to-date. A sketch of the 
sinking rig finally used appears in Figure 2. 
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ANALYSIS OF HELICAL BEAMS UNDER SYMMETRICAL LOADING* 


ALAN M. C. HOLMES, ! J.M. ASCE.—The writer thanks Mr. Hall for his 
discussion of the paper. Certainly the term “Modified Polar Moment of 
Inertia” leaves something to be desired, even if only brevity. 

Numerous terms were considered and debated with fellow engineers, 
among them “Torsional Moment of Inertia” and “Torsion Factor.” It was 
decided that the latter was not specific enough. The term gives no clue as 
to its units and could be used to apply to “x” (in the original paper), or to 
the empirical constant used in the torque-shear stress relationship for 
rectangular cross-sections. 

The writer is not aware of any established practice in this respect, and 
there was no intention of dictating what the term should be. 

Thanks are due to Mr. Bergman also for discussing at such length the 
reasons why the solution presented is not “exact.” A service has been ren- 
dered to the engineer in reminding him of the folly of thinking his answer 
“exact.” However, it seems Mr. Bergman misinterpreted the statement, 
that moments calculated by the method presented in this paper can be 
mathe matically “exact.” 

The words “can be” were intended to mean that there are certain cases 
where (under ideal conditions) the moments could be exact. A wire of small 
diameter bent into a helix of fairly large radius and securely clamped at 
both ends could be such a case. Surely the words “can be” do not mean “are 
always.” Mr. Bergman repeatedly used the word “exact” in quotes to denote 
that it is a manner of speaking. Reference to the original text will show that 
the word “exact” is in quotes. 

“A precise MATHEMATICAL solution” was the writer’s method of differ- 
entiating between a method that takes no “mathematical shortcuts” and one 
which uses a rough mathematical approximation, so that even if the “input 
data” of the problem were reasonably accurate, the answer could never be. 
The approximation of a helix by a ring lying entirely in one plane is an in- 
exact mathematical procedure. The intention was to present a mathematical 
tool, applicable to a fairly wide range of situations, and with limitations to be 
determined by the user. It seemed hardly necessary to point out that the 
“output data” could be no more accurate than the “input data.” 

The value of K is perhaps questionable, but pointing out that the stress is 
higher at the inner edge of a curved beam than the value calculated for the 
stress in a straight beam throws no light on the subject. 


a. Proc. Paper 1437, November, 1957, by Alan M. C. Holmes. 
1. Structural Designer, Skidmore, Owings & Merrill, Chicago, II. 
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For a curved beam such as the one Mr. Bergman describes (width = 
centerline radius, and presumably lying in one plane), it can be shown 
(Figure 17) that even though the stress is higher at the inner edge by 52%, 
the value of K (“the measure of resistance to rotation per unit length due to 
unit moment”) remains essentially the same. For a helical beam where the 
curvature is more complex, this may still be an approximation, but it is con- 
ceivable that there is some value of K which it would be correct to use in the 
solution, and that this value could be obtained by mathematical or empirical 
methods. For the small deviations that may occur in the K of a helical beam 
of reasonable proportions, it is simpler to recognize that the answer is ap- 
proximate and to make an allowance for it, or to determine the variation in 
maximum moments for the range over which K is likely to vary. It can be 
said from experience that the maximum calculated moments vary less (from 
this cause) than we know concrete to vary in its properties. 

On the basis of this solution, and the model test which showed close agree- 
ment, the writer has designed eight circular stairs (see Figs. 18 and 19). 
These are now built and performing satisfactorily. The largest, which has 
an arc of 360°, an outer diameter of 24 ft, and slab width of 7 ft, was recently 
load-tested. Under a total load of 30 tons (uniformly distributed), the mid- 
point deflection was 0.092 in., which is in line with the theoretical value. 
Strain gages indicated that stresses were about those anticipated, although 


a precise comparison is neither possible nor justified because of extraneous 
factors. 
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Fig. 18. Free-standing spiral stair at the Air Force Academy Viewed 
from below. 
Arc: 360 degrees Outer diameter: 24 ft. 
Height: 18 ft.4 in. Inner diameter: 10 ft. 


Free-standing spiral stair ut the Air Force Academy Viewed 
from side. 

Arc: 290 degrees Outer diameter: 24 ft. 

Height: 14 ft. Inner diameter: 10 ft. 


& 
Fig. 19. 
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A METHOD OF DESIGN OF REINFORCED CONCRETE SECTIONS* 


Discussion by Juergen Plaehn 


JUERGEN PLAEHN!.—Mr. Moliotis presents a method for the design of 
any section of reinforced concrete. This method reduces the calculation of 
any section with a symmetrical axis to that of the rectangular section. 

The writer believes that the derived design formulas are difficult to sur- 
vey. The designing engineer runs the risk of seeing only the factors and 
formulas and not seeing the basic ideas on which his work is based. This 
results in misinterpretation and reading and writing errors with all of their 
consequences. Such difficulties can be avoided, if the design problem is 
solved step by step instead of using implicit design formulas. With the 
generally accepted assumptions for the design of reinforced concrete sec- 
tions in mind each problem with one symmetrical axis can be solved quickly. 
One table set up for rectangular sections can be used for several purposes 
without having the user confronted with inconceivable coefficients. Table 1 
of the “Reinforced Concrete Design Handbook” of the American Concrete 
Institute may be used as the only necessary table. A similar table which is 
more adaptable to slide rule computation is given below. 


The assumptions are: 1. The Bernoulli Hypothesis that plane sections remain 
plane. 


. The Hooke Law that stresses are proportional to 
strains. 


. The Straight Line Law for both composite 
materials, derived from assumption 1 and 2. 


. The concrete section is cracked in the tensile area. 
The steel has to resist all tension forces. 


Particularly assumption 2 (therefore also 3) is not in agreement with the be- 
havior of concrete in reality the modulus of elasticity for concrete is not con- 
stant and decreases considerably at the breaking strength. The assumed 


~ 


n= EF differs more or less from the real value. Therefore it has no sence 


to exaggerate the accuracy of the calculation. 


The presented design table is set up according to the following derivation. 
Similar tables for other steel stresses or other n-values can easily be 
computed. 


a. Proc. Paper 1509, January, 1958, by Panagiotis D. Moliotis. 
1. Asst. Design Engr., Richardson, Gordon and Associates, Pittsburgh, Pa. 
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y Ys 
.. 
h d z defined: 
x = (1) 
(2) 
fs + 


from (1) and stress strain diagram: 


nef +f 


from (2) and stress diagram: 


z= (d-x/3) = d(l-k/3) = jed (2") 
it is: Ges (4) 
M = Tez (5) 


from (4) with (2'): 


Mi = = bek(1-k/3)f 


2 M M 


from (5) with (2'): 


M fg(d-x/3) = f.+d(1-k/3) (3") 


from (5') with (6): 


- r (8) 


from (8) with (6): 
M t 
A, = = = p*bed (9) 


The examples show how to use the tables for various purposes. 
If a given unusual flexural member is to be designed, the neutral axis can 
be determined for the allowable stresses from the lineary stress diagram 
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nx f, = d 

nx f, +f, 
each section element the resultant C; of the compression stresses and its 
moment Mj, = C; x z; about the center of gravity of the tensile reinforcement 
can be computed. 2 Mj; is the resisting moment of the concrete stresses, 


and LC, = T© is the force in the tensile reinforcement. If the external 
moment M exceeds the resisting moment of the concrete stresses, the 


= . The table can be used for this purpose x = kxd. Then for 


residual moment AM= M-2M, must be resisted by additional compressive 
and tensile reinforcement pst = cSt is the distance between 


the centers of gravities of the compression and tensile reinforcement. Then 
the necessary amount of reinforcement can be determined because the 
stresses at the center of gravity of the tensile as well as of the compressive 
are known from the stress diagram. 

If bending and axial load are combined the same procedure is used when 

“ > z. The external moment will be taken about the center of gravity of 
the tensile reinforcement; therefore, the line of acting of the normal force 
coincides with the center of gravity of the tensile reinforcement. The stress 
diagram is the same as for bending alone. The amount of the resultant ten- 
sile force will be reduced, if N is a compressive force, and increased if N 
is atension force: T = T° + TSt 4 N. If the resisting moment of the con- 
crete stresses exceeds the external moment M<Mj, the cross section or 
the concrete fiber stress may be reduced. By interpolation between two com- 


puted values of A, or f, the accuracy of the calculation is usually 
sufficient. 


Resultant of the Stresses for any Section or Section Element: 


Acompression 


dif >> | 
& 
i i i 
f. 
(10) 
i Ay 


October, 1958 


Location of the resultant: 


* 


Element -i 


Special Cases: 


From equation (10) is seen that for any section or section element the 
resultant of the internal forces is equal to the area of the section multiplied 
by the stress at the center of gravity of the section Cj = Aj = fo,i,C.G.° The 


distance of the resultant from the neutral axis can quickly be determined by 
dividing the moment of inertia by the section modulus, both in respect to the 
neutral axis. See equation (11). In most cases it is possible to set together 
the elements of the area and of the stresses diagram out of rectangular and 
triangular elements, - the location of the resultant is then known, see (12). 
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n= 15 


0.00138 0.939 
0.00164 0.934 
0.00191 0.92 
0.00221 0.924 
0.00251 0.920 
0.00284 0,916 
0.00318 0.912 
0.00353 0.908 
0.00389 0.904 
0.00427 0.900 
0.00465 0.896 
0.90595 0.893 
9.00546 0.890 
0.00588 0,836 
0.00631 0.883 
9.00675 0.88C 
0.00720 0.877 
0.00755 0.874 
0.00812 0.871 
0.00859 


0.00907 0,866 
0.00955 0.863 
0.01004 ‘,860 
0.01054 0.858 
0.01105 

0.01156 

C 


0.01261 
0.01314 
0.01367 
0.01421 
0.01475 
0.01530 
0.01586 
0.01642 
0.01698 
0.01755 
0.01812 
9.01869 
0.01927 
0.01985 


uj 


0.02281 
0.02402 
0.02463 
0.02524 
0.02585 


0.808 
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0.868 


ol 


20,000 1b/in? 
300 21.540 0.639 
330 19.813 0.642 b Lin] A 
360 18.371 0.646 
390 17.150 | 0.649 oa fin | 
420 16.101 0.652 = 
450 | 15.191 | 0.65 
4380 14.393 0.658 
540 13.060 0.664 
570 12.498 0.666 der 
600 11,990 0,669 
630 11.530 0.672 
660 11.111 0.674 
7 690 | 10.727 | 0.677 ad a 
720 10.375 0.679 
750 | 10.050 | 0,682 p= As 
80 5.750 0.004 bd 
810 9.471 0.686 
840 9.212 0.689 
870 8.970 0.691 
3900 9.743 0.6923 
30 8.530 0.695 
960 | 92331 | 02697 
y90 8.143 0.699 
1050 7.798 0.703 
O80 2640 2 
1110 7.489 0.707 0.454 0.848 
1140 7.346 0.709 0.461 0.846 
1170 7.210 0.711 0.467 0.844 
1200 O80 0.712 0.474 0.842 
1230 6.957 0.714 0.480 0.840 
1260 6.840 0.716 0.486 0.838 
1290 6.727 0.718 0.492 0.826 
1350 6.516 0.721 0.503 0.832 
1380 6.417 Oeloo 0.509 0.830 
1410 6.321 0.724 0.514 0.829 
1440 6.230 0.725 0.519 0,827 
1470 6.142 0.727 0.524 0.825 
1500 6.058 0.728 0.529 0,82 
1560 5.899 0.731 0.02103 0.539 0.820 
1590 5.822 0.733 0.02162 0.544 0.819 
1620 5.749 0.734 0.02221 0.548 0.817 
1650 5,678 0.736 9.553 0.816 
1710 5.544 0.738 0.562 0.813 
1740 5.480 0.739 0.566 0.811 
1 1770 5.418 0.741 0.570 0.809 
1800 5.358 0.742 0.574 a 
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Example 1: 


given: M = 320 Kft 
17 

2 
fs ,e11." 20,000 1b/in 


" 


fc,all.” 1,200 1b/in 
25 


determine: As 3 


from table find for £220,000 1b/in® and f,= 1,200 1b/in®: 
r= 7.080 ; l/a=0.712 ; j = 0.842 


a=r\VE = 7.08 255 = 30.7 in 


20 


Mil 2 
As= ada z 0.712 = Ta49 in 


ze 0.842+30.7 = 25.9 in 


Often it may be more economical to use a f, smaller than fc,aj).; this 
particular if the widta is large. 


Example 2: 


given: M = 300 Kft 
b= 50 in 
15 in 
d= 10 in 
2 
1,200 1b/in 


2 
911.7 20,000 1b/in 


a 
deter mine: d ; A 


used: 420 lb/in? 1/3 
r=16.101 ; l/a=0.652 ; k = 0.239 3; Jj = 0.920 


= 16.101 = 39.4 in 3; x = 0.239+39.4 = 9.44<10 


Ag= 0.652 = 4.97 in® = 04920+39.4 = 36.3 in 
zo 


| 

bee 
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fe given: M = 140 Kft 
N = 100 K 
E = 14 in 
in 
ad = 31 in 
d*= 2 in 
2 
20,000 1b/in 
1200 1b/in 
n= 15 


determine: 4, + f 


Moment about the tensile reinforcement: 


M* = 140 + 100+#14/12 = 256.5 Kft 
a 31 

r= 1[ 256.5 
17 


f= 1020 1b/in® 


Me N 256. 100 2 
= 0.701 - = 0.7 in 
Example 4: 
given: 
= 500 Kft 
| dt=xd''= 2 in 
= 20,000 1b/in® 
s,all. 
st rete 2 
determine: A, AL 3 


with fo all. and fs,all. find from table the resisting moment of the concrete 
stresses: 


«= 7.080 3; = 0.712 3: k 0.474 


Example 3: 
M 
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to M, belongs the tensile reinforcement A 


_ 326 
O.712 = 7.5 in 


x = k*d = 0.474"31 = 14.7 in 


the residual moment 4M = (M - M,) must be resisted by additional tensile 
and compression reinforcement: 


AM = 500 - 326 = 174 Krt 


= = 2 = 2.42 ft 


total tensile reinforcement: 
A, = 7.5 + 3.60 = 11.10 in” 


compression reinforcement: 
from statical moment about the neutral axis find As: 


d -x 
x- 


AS = 3.60 « = 14-7) s 4.62 in? 
(14.7 - 2) 


because: orf, = Ag 


of, 
15.55 K/in* 


If moment and normal force are combined the same procedure is used. The 
moment about the tensile reinforcement will be determined: M' = M + N=E. 
The tensile reinforcement is found from: 


Ag = Ag; + + As,3; where Ag 3 = 


M 174 
A = - = Te = 3.60 in* 
Foal, | 
Ag 
2 
f° = 
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Example 5: 


given: cross section as shown above, 


£. = 1,900 lb/in‘ ; M = 1,800 Kft 
20,000 1b/in H N = K 
n = 15 ; E= 2 ft 
assume: =» 5 im) « 10 in 
M* = M + NeE = 1,800 + 50022 = 2,800 Kft 
64«15* 
x a « 33.9 in 
2000  £5* 


Zor 67 = 15 = 52 in = 4,32 ft 


AM = 2,800 - 2,358 = 442 Kft 3; Co, = To, = 442/4.32 =102 K 


St 
ST = 603 + 102 = 500 = 205 K i 
= 20.0°28,1/30.1 = 18.65 K/in‘ 
s,C.G. 


f' = 15.8 K/in* 
s,C.G. = 

A_ = 205/18.65 = 11 in® 


At = 102/15.8 = 6.45 in® 


0.172 
| 0.394 
4 
i 
f a 
r 
fe 
[in*) [in]] K/in ~ un Kft 
5850 2202500 Je515 30 45,0/ <4) 
0.295 148 | 23.90 | 52.00 640 
1329 *10=139 6.95 | 0.305 43 9.25 37235 133 
III] 2°3*5*10= 50] 12.23 | 0.394 20 | 12.23 | 40.33 66 
0.148 7} 12.65) 40.75 25 
IV] 2*10*10=200] 8.90] 0.172 34| 8.90] 37.00] 106 
0,22] 4 10.57 | 38.67 143 
Cc, = 603 DM, = 2358 
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Second example from Mr. Moliotis: 


”~ fry given: cross section as shown 
850,000 1b in 
900 1b/in® 


f 20,000 1b/in? 


n = 15 
determine: A, 


K n in 
My 
12«8*3=48 | 5.96 |0.476 22.9 | 6.56 | 23.26] 45.2 
3¥12=394 1. ee 0.132 5.2] 2.20 {18.90 8.2 


AM = 9304000 - 53.4 = 17.6 Kft 


pot. 17.6.12 9.2K 


= 1.40 + 0.46 = 1,86 in® 


At = 0.46 ( 28 - 11,3) = 1.22 in® 


= C 


Mr. Moliotis gives a revised value: AS = 1.22 in?; the above method can 
be used for the remaining examples given by Mr. Moliotis. 


RE FERENCE 


1. “Bemessung der Stahlbetonquerschnitte” Prof. Dr.-Ing. Dischinger 
Technische Universitaet Berlin-Charlottenburg, 1948. 


PANAGIOTIS D. MOLIOTIS,! M. ASCE.— Corrections: 
Page 2: Equation (2) should read: 


Page 4: Equation (9) should read: 


1. Dr. Civ. Eng., Asst. National Technical Univ., Athens, Greece. 


| 
4 
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Page 7: Equation (33) should read: 


& F | 
— 
(14+ fs_) 
\ Page 17: Lines 5, 6 and 7 should read: 2 
<5 (1- ) 
04034 
603 C 1270001 
) 
28 
1? is 
A.— 0542 103 X 28 + 1 27x 
4 


‘ 

" 
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ANALYSIS OF BRACED FRAMES®# 


Discussion by I.M. Nelidov 


I. M. NE LIDOv,! M. ASCE.— In this paper the author has presented a 
method of analysing frames with hinged connections of the braces, based on 
the principles of H. Cross and deflection stiffness. A sample problem was 
given of a portal assymmetrical frame with hinged braces, loaded with a uni- 
form load over the top member. 

The writer has done some work on similar frames with braces with rigid 
connections, where he used Cross’ method and potential deformations of 
the frame, for more visual understanding of the deformed shape of the frame. 
In this discussion, the writer will repeat the example given by the author, 
using his own approach. 

There are two steps in this analysis, one with lateral sway prevented, and 
another, with lateral sway. Referring to the figure showing the frame on 
p. 1560-8, the joints B and C and the joints at the lower end of knee braces, 
which we will call B' and C', are held from lateral displacement, while the 
uniform load of one k/ft. is applied to the top beam. 

In this case, the intermediate supports of the top beam between joints B 
and C are also immovable and the top beam acts as a continuous beam of 
three spans and four supports. The fixed-end moments and end shears are 


2 
respectively, for span 10.ft. M = 1. 0. = 8.33'k and V =5.k.; for span 


12 
20.ft., M =1. ~ = 33.33 'k and V =10.k. The balancing of the moments 


is shown below, with stiffness values for end spans .75 and for the center 
span .50, giving distribution factors respectively .60 and .40. 


FEM. 8.33 8.33 -33.33 933.33 -8.33 28.333 
48.33  -8.33 

12,80 8.53 --8.53-12,80 
0,00 0,00 


4026 026 
72.55 -1.71 


286 286 
34 3h - 252 


o17 
210 207 - 207 = 


= 204 04 
202 202 202 202 


27099 427499 -27499 


a. Proc. Paper 1560, March, 1958, by K.H. Gerstle. 


1. Structural Engr., San Francisco, Calif. 
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The shears resulting from these moments are respectively for left and right 
supports of the end span at left: V; = -2.80k and V2 = +2.80k, the sign plus 
indicating upward direction; for center span: Vj = V2 = 0; for right hand span: 
Vi =+ 2.80k and V2 = -2.80k. The reactions of the beam at consecutive 
supports starting from left are Rj = + 2.20k, = Ry, Rg = 17.80k = R3, positive 
for upward direction. The lateral horizontal forces induced by these reac- 
tions thru the braces and applied at points B' and C' are: -11.86k and 17.80k, 
at minus sign -, acting to the left, and with plus sign -to the right. Actually, 
when the joints B and C are held immovable, the joints B' and C' move later- 
ally due to the action of the uniform load transferred thru the braces; there- 
fore, we will deflect arbitrarily the joints B' and C' firstly outward, i.e. 
joint B' to the left and joint C' to the right an amount A; this will displace 
the intermediate supports of the top beam an amount .667 A at left knee- 
brace and A at right knee brace. If we assume an arbitrary fixed-end 
moments in the right-hand column DC', M = -10.'k., then at the upper part 
of the same column the moment will be with one end hinged M = +5.'k, at the 
lower end of column C'C. For the left -hand column, since its height is 1.5 
times that of column DC, the fixed-end moments in AB' will be M = + 4.44'k, 
and the moment at the lower part of B'B will be M = -2.22'k. The fixed-end 
moments in the left span of the top beam BC will be M = -6.67'k; in its cen- 
ter span M = -.83'k, taking into account that this span is twice the height 
of column DC'; in the right -hand span M = +10.'k. Signs of the moments 
are positive when clockwise. 

The rigidities of the top beam BC being previously computed, the balancing 
of the moments is as given below. 


-6.67 -.83 -.83 #10,00 
+6.67 23.33 5,00 
2.50 41.67 -1.67 - 2.50 


0,00 


#50 - 33 = .50 
016 216 
+.10 206 - 206 =- Pete) 


203 203 
202 ? = = 202 


t e22 -1.88 41.88 


The shears resulting from these moments are, for left span: Vy = +.02 
and V9 = -.02k.; for center span: Vi= -+.08k and Vo = -.08k; for right span; 
Vi = -.19k and Vg = +.19k. The respective reactions respectively from the 
left supports are: Ry = +.02k., Rg = +.06k, Rg = -.27k and Rq= +.19k. The 
horizontal forces exerted thru the knee-braces at their lower ends due to 
these reactions are, at the left joint B': -.04k and at C' -.27k. The moments 
in column AB due to the fixed-end moments resulting from the displace ment 
A are obtained thru the balancing shown below, with stiffness of AB' taken 
as 1.00 and of B'B as .75, giving the distribution factors respectively .57 for 
AB' and .43 for B'B; the column is laid horizontally for convenience. For 
column DC' the balancing is shown below at right. 


0,00 10.00 10,00 > 5.00 0,00 
*2 


“e 63 -1.26 96 1.43 2.85 


43.81 03.18 8,57 7.15 * 7.15 


#10.00 
-10,00 
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The shears for the left column are, for AB': V, = +.46 and Vo = -.46k, posi- 
tive acting to the right, considering vertical column; for part B'B: 

Vy = -.21k. and Vo = +.21k. The respective reactions in the column are: at 
A Rj = +.46k, at B' R2 = -.67k and at B Rg = +.21k., positive acting to the 
right. 

For the column DC, we will have: for DC', Vy = -1.57k. and V2 = 1.57k; 
for C'C, Vy = +.72k. and V2 = -.72, at C' and C respectively; the resulting 
reactions are: at D Ry = -1.57k, at C' Rpg = + 2.29k and at C Rg = -.72k. 

Since the lateral deflections at joints B' and C' are not equal, these 
points should be deflected both in one direction, so as finally to achieve the 
inequality. We arbitarily deflect these points to the left an amount A; this 
will induce an upward displacement of the upper end of the right-hand knee 
brace of the amount A; at the upper end of the left-hand knee-brace the dis- 
placement will be downward to the amount .667A. If, as before, the fixed-end 
moments in the right -hand column DC' are taken arbitrarily at +10.00'k, the 
moments in the left-hand column AB' will be +4.44'k.; the moments in the 
lower end of the upper half of the columns respectively will be, in C'C M = 
- 5.00'k. and in B'B M = - 2.22'k. The moments in the top beam in its left- 
hand span will be: M = -6.67'k, in its center part M = +4.17'k and in the 
right-hand span M = -10.00'k, all fixed-end moments. The balancing of 
these moments in the beam BC is as given below. 


~6.67 4.67 #44217 -10,00 =10,00 
+6.67 23.33 5.00 #10,00 


,02 « .dl .02 


=3.96 43.96 4.38 


The resultant shears are, in the left -hand span: Vy = +.40k and 

Va = --40k; in center span: Vj = -.42k and Vo = +.42k; in right -hand span: 

V1 = +.44 k and Vo = -.44k. The respective reactions are Rj = +.40k, 

Rg = -.82k, Rg = .86k and Rg = -.44k. The forces induced thru these reactions 
horizontally at points B' and C' at lower ends of the braces are: at B': +.55k 
and at C' + .86k. 

Since the lateral deflection of the column AB in this case is the same as in 
the previous case, the moments, shears and reactions in this column are the 
same as were obtained then. For the right -hand column the deflection was 
taken to the left or in opposite direction compared with the previous case; 
therefore, the signs for the moments, shears and reactions should be re- 
versed. In order to determine the unbalanced force of reactions at the 
bases of the columns due to the imposed restraints of the joints, we will 
write the equations from which the coefficients can be determined, multiply- 
ing by which the data obtained in this step, the unbalanced force could be 
found. These equations are: 


~067 C =.67 C, - for the left-hand column, and 


#17.80 ~.27 C) + .86 Co = 2629 2.29 C, - for the right -hand column. 


| 


1827-26 ST 6 October, 1958 


By solving these, it is found that Cy = 11.58 and Co = 3.75 

The unbalanced force is therefore: H = .46 (11.58 + 3.75) -1.57(11.58-3.75) 
= -5.25, acting to the left; which means that the frame is acted upon with the 
excess force acting to the right, and therefore, we must deflect it primarily 
to the left, for equilibrium. 

In achieving this, we must firstly to deflect the lower parts of columns to 
the left, so that the upper portion of the frame would remain undistorted; 
then, to deflect the upper parts of the columns, with the inevitable distortion 
of the beam; and then, deflect the columns outwardly, as was done in the 
second part of step one. 

By deflecting the lower part of column AB, namely AB’, to the left, the 
fixed-end moments induced there are M = +4.44 'k, if the moments in the 
right-hand column DC' for the same deflection are taken at M = + 10.00'k. 
In the upper part of both columns the fixed-end moments are zero. Below is 
shown the balancing of the moments for both columns. 


0,00 0.00 +10,00 #10.00 0,00 0.00 
=1.26 -1.91 - 2.85 5-70 4.30 


#7215 4.30 


The corresponding shears in the left hand column are: in AB' Vq = +.34k 
and V2 = -.34k; in B'B Vy = -13k and Vg = +.13k. The reactions are, at 
A Rj = +.34k, at B' Rg = -.47k, and at B Rg = +.13 k. 

The shears in the column DC are: in DC' Vq = +1.15k and Vg = -1.15k; 
in C'C Vq = -.43k and Voy = +.43k. The reactions are, at D Ry = + 1.15k, at 
C' Rg = 1.58k. and at C Rg =+ .43k. 

In deflecting the upper portion of the columns to the left an amount A, the 
top of the left brace will rise .667 A, while the top of the right brace will go 
down an amount A. This will produce the moments, shears and reactions 
similar to those obtained for the beam BC in the second part of step one but 
with opposite signs. The moments in the upper parts of the columns with 
fixed end condition will be correspondingly for the right and left columns 
M = +10.00'k and M = +4.44'k. The balancing of these moments is given 
below, for the left and right columns respectively. 


6,00 0.00 the 0,00 0,00 +10,00 
-2 e22 mhe 4 = 5 . 00 

263 1.26 296 2.85 = 2015 


=.63 =1,26 01.26 - 2.85 22-85 


The shears in the left-hand column are: in AB' Vq = + 13k, and V2 = -.13k; 
in B'B Vq = +.08k. and V2 = -.08k. The reactions are: at A Ry = -.13k., 

at B' Rg = +.21k. and at B Rg = -.08k. The shears in the right hand column 
are: in DC' Vy = -.43k. and Vo = +.43k.; in C'C Vqy = +.29k and V2 = -.29k; 
The reactions are: at D Ry = -.43k, at C' Ro = +.72k. and at C Rg = -.29k. 
In deflecting the columns outwardly, the same data are obtained as were 
listed in 2nd part of step one. The equations necessary for determination of 
coefficients Cg, C4 and Cs are given below. 


q 
? 
22013 41.92 
#10,00 
-10,00 
0,00 
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Cy Co C, C, ~.67 - for the lower end of left-hand 
column, and 


=186 C, Cy = “1058 


2e29 C. - for the lower end of right 
hand column and 


5 


5 3 2 


The solution of these equations renders the values of coefficients C as 
follows: 


Cy = 5.88, C4 =4.38 and Cs =.91. 


After substituting these coefficients in the last of previous equations, it was 
found that the equation is satisfied. 

In order to demonstrate how moments and reactions are obtained with the 
aid of the above coefficients, a sample computation will be given for joint A. 


Thus Mg = 3.81 (11.58 + 3.75) + 3.18x5.88 -.63x4.38 + 3.81x.91 = 58.5 
+ 18.70 - 2.76 + 3.46 = 77.9'k. 


Ry = .08x11.58 - .42x3.75 + 0.00 + .42x4.38 + .08x.91 = .93 - 1.58 + 0.00 
+ 1.84 + .07 + 20.00 = 21.26 k. 


Other moments and reactions can be obtained in a similar way, and they check 
the values given by the author in his figure on p. 1560-11 within 1%. 

In concluding, the writer wishes to say that while the method of 
analysing braced frames by following their deformations may by not as 
short as that offered by the author, in the estimation of the writer it gives a 
visual process of realizing the deformed shape of the frame. 


NUMERICAL SOLUTIONS FOR BEAMS ON ELASTIC FOUNDATIONS# 


Discussion by Victor R. Bergman and N, J. Ryker, Jr. 


VICTOR R. BERGMAN, ! A.M. ASCE.—Mr. Malter has chosen to deal ex- 
clusively with the mathematical manipulations involved in the solutions, and 
with this decision there can be no argument. Still, in view of Mr. Malter’s 
own statement that the phase covered in his work is only one aspect of a 
tripartite problem, it may be admissible to recommend for the benefit of 
designing engineers who must fact the problem in its entirety, a paper by 
Karl Terzaghi, published in 1955.2 

In his laudably concise presentation of the application of numerical methods 
to the solution of beams on the so-called “Winkler foundation,” Prof. Malter 
has furnished his readers with two references , 7) for further study of the 
finite-difference method. This technique, as applied to the problem considered, 
is of such beautiful simplicity that consultation of references is unnecessary. 
The only feature of the procedure that may prove unfamiliar to the structural 


analyst is the replacement of the second derivative in the familiar beam equa- 
tion, = , by its finite-difference approximation, De 2yi + ¥1 


A simple derivation of this approximate expression is offered in Fig. 1. 


N. J. RYKER, gr.3,— In this paper Professor Malter has outlined in ex- 
cellent detail two numerical methods for solution of beai.s on elastic founda- 
tion. The paper will undoubtedly increase the popularity of the powerful 
successive approximation methods for the solution of these problems. The 
earlier work of Professor Popovi1) dealt primarily with the graphical solu- 
tions and therefore may have had limited acceptance by those who abhor 
graphics (an attitude apparently growing in popularity). Professor Popov 

did mention that his calculations could be set up in tabular form. 

The second innovation by Professor Malter and one which can contribute 
greatly to rapid convergence for long flexible beams, is the initial assump- 
tion of an approximate deflected shape rather than the rigid shape. Contrary 
to the authors statement, however, the successive approximation method us- 
ing the Newmark type solution is not restricted to relatively short beams. 
It is the purpose of this discussion to present a method(2) that will assure 


. Proc. Paper 1562, March, 1958, by Henry Malter. 
. Vice-Pres., Godwin Constr. Co., New York, N. Y. 
. “Evaluation of Coefficients of Subgrade Reaction” (Geotechnique, Dec., 
1955). 


3. Preliminary Design Engr., Missile Division, North American Aviation, Inc., 
Downey, Calif. 


ASCE 1827-29 
a 
| 
2 


1827-30 ST 6 October, 1958 


2 
Reguired: Ap expression vor 
LY x0, in perms of or- 


x 


Derwarror: 


The change of slope X = Z 


rate of Change G7 ¢s 


fia 1 Derivation of Ditferene Lxpression fer 


| 
| 
| 

4 

| 

é = 


ASCE 


DISCUSSION 


1827-31 


convergent solutions for very long beams on elastic foundations. 

In general, the method proposed by Professor Malter will provide a con- 
vergent solution as long as AL = 3, where L is the length of the beam and 
A4 =k/4El. For beams of greater length a convergence factor may be used 
to obtain a convergent solution. The general procedure is as follows: 


1. Assume an original deflected shape and place the beam in equilibrium 
with the applied loads and moments by rigid beam deflection and rotation. 
Let the balanced beam shape be denoted by yo. 


2. Using the resulting foundation pressures and the applied loads, cal- 
culate the new deflected shape by any convenient method. Add rigid beam 
deflections and rotations so that there is no net load or moment on the beam 
(remove all applied loads). Let this deflected shape be denoted by 6}. 


3. Calculate the deflected shape due to the foundation pressures of 6, and 
adjust as in step 2. Let this shape be denoted by 69. 


4. Calculate the quantity r = b “5 at each panel point of the beam, 
1- 


and average for the entire beam. 


5. The final shape is then: 
* 2, + r 


6. If the quantity r ) is large or if the convergence factor r varies 
greatly along the beam repeat steps 1-5. 
The development of this procedure is quite simple. In successive approxi- 
mation methods the solution can be represented as: 


where Yo, the original assumed shape, balances all applied loads and 64, 69, 
etc., which are corrections to this shape, are in static equilibrium without 
applied loads. As long as éy and all successive 6,’s are similar in shape, 


then 

and 

y*y,+ 


The quantity r can be calculated from two successive cycles; to minimize 
calculations it is proposed that only the first and second cycles be used even 
though these shapes generally are not similar. 

7 * 


Y"¥o*+ 8, +r 


Equating these two quantities, the convergence factor in step 4 is obtained. 
The validity of this procedure is dependent upon the similarity of 6, and 
69; since in general they will not be similar, the second of the two equations 
above will more closely approximate the true shape. 
Mr. Karol(3) in his discussion of Professor Popov’s paper suggested us- 
ing: 
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where P, = total reaction under 6}, if the ends of the shape 6, are placed at 
the surface of the foundation. 
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Po =total reaction under 69, if the ends of the shape bo are placed at 
the surface of the foundation. 


This equation is based upon the end pressures of 6; and 69 because for a 
constant foundation modulus 


Py x20 + | 


where qj is the foundation pressure for 6,;. Because it is based solely upon 
the end values, it fails to give a convergent solution in some cases that the 
suggested method does. 

Mr. Plumb, (4) also in a concurrent development of a convergence factor 
similar to Mr. Karol’s continued his successive approximation process until 
the ratios of successive differences were a constant. However, in the case of 
beams of AL = 6 the procedure must be repeated for at least four cycles and 


requires the use of 6-10 significant figures because the answer is on the 
order of 10-4 64. 


The general procedure outlined above can be applied directly to the pro- 
cedures outlined by Professor Malter. The principal difference being that 
64 and é9 are not obtained directly. 


1. Follow Malter’s steps (a) through (m) and then using M, calculated 
in (m) repeat steps (i) through (1). 


2. Correct the original deflected shape yg and the two successive values : 
y to place them in equilibrium with the applied loads using Malter’s 
correction configurations. Let these successive values be yo, y;, and 


3. Calculate 6; = yy -Yo and = yg - yq- 


To illustrate this procedure consider the first problem presented by 
Professor Malter. The first cycle remains as shown. However, the actual 
calculated values of Mp must be used to initiate the second cycle from which 
the deflection, shears and moments shown in Table I were calculated. 

After the moment and shear at the end of the beam, rows 3 and 4, were 
calculated, equal and opposite moments and shears were applied to balance 
the deflected beam, rows 5 and 6. The rigid deflections of Malter’ S correc- 
tion configurations ior shear (Yp) and moment (y A) were added to y' » row 2, 
producing yg. Similar corrections were computed for ya and y' in the first 
cycle producing y, and yy respectively. From the differences, the conver- 7 
gence factor r was computed. The variation of r along the beam indicates 
that 6, and do are not similar, but the average value of r = .788 produces 
better agreement with Hetenyi’s solution than does yo. It is apparent from 
the values of y2 and yj, that, because of the idealization of the beam, the 
solution will never exactly coincide with Hetenyi’s solution. 

A better example of the rapid convergence is supplied by Malter’s second 
example, table II. In this case using one cycle of the general procedure 
(which is equivalent to two cycles of Malter’s procedure) and an average 


| 
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TABLE I - Example 1 


Point 
Second Cycle A B Cc D 
Mo, h°/12 Ib-in. 14.82 -17.69 0 
y' , 1072 inches +0.819 
V, h/l2 lbs =9 90 -4.17 29.9 
M, h2/12 lb-in. fe) -9 290 -14.07 -17.29 
y'cl=(- Yp -0.512 40.128 +0.769 41,261 
y‘ = ( Ee) 0.221 40,025 04172-0369 
Balanced Deflections 
+ 40.528 +0,.153 +0.597 +1.859 
+0.580 +0.123 +0.59h 41.982 
Yo + +0.193 40.629 41.393 
$1°%1-Yo 40.126 -0.070 -0.035 +0.589 
-0.952  +0.030  +0,003 <1.23 
+0.708  +0.700 +0.920 +0.826 
Results 
y 0.539 0.147 0.596 +1.885 
y (Malter corrected) 0.528 0.153 0.597 +1.859 
y (Hetenyi) 0.517 0.114 0.53 +1.39 


*A numerical error in Malter's analysis has been corrected producing these 


results. 
TABLE II - Example 2, Foundation Pressures 
(1b/ft) 
Point 

A B c D 
Popov 2e7 L.5 59 6.5 
Valter (3 cycles) 2.5 Lek 5.9 6.7 
Ryker (2 cycles) 2.5 Leb 5.9 6.62 


r 0797 2859 2666 
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convergence factor of 0.710 results were obtained which were equivalent to 
Malter’s three cycles. 

In Malter’s third example convergence was less rapid because of the ex- 
treme flexibility of the beam. An original shape was assumed which seemed 
consistent with the applied loading and this was balanced producing the 
original shape y,. After two cycles the convergence factor was computed 
and a resulting shape computed. Another cycle was computed because the 
wide variation in the factor, and the result produced a maximum deflection of 
.057 inches which is comparable to the value of 0.056 shown by Malter. It is 
interesting to note that Malter could only obtain the solution to the problem 
by averaging the results; this is equivalent to using a convergence factor of 
0.5 and in this case after two cycles the computed average convergence fac- 
tor was 0.55. 

In extending the method to very long beams,(2) a comparison was made be- 
tween this numerical procedure and the exact solution for centrally loaded 
beams of AL = 8. In this comparison the maximum error in deflection at the 
end of the beam after 2 cycles of the general procedure (equivalent to 4 cycles 
of actual deflection calculation) was less than 5 percent of the maximum de- 
flection. The maximum moment was identical to the moment calculated by 
Hetenyi’s method. Without the use of the convergent factor the successive 
approximation solution was rapidly divergent. The calculated value of the 
convergence factor was .083 indicating that an averaging technique would also 
be unsatisfactory. 

The purpose of this discussion was to show that the use of a convergence 
factor(2,3, or 4) enables one to obtain convergent solutions by successive 
approximations for beams lengths of all practical interest. That is beams 
with a free length beyond the applied load or restraint of not more than 
AL = 4. As shown by Hetenyi(5) there is no appreciable effect on the beam 
at greater distance. Longer beams can be handled by assuming free beam 
length of AL = 4 for analysis and assuming the rest of the beam is unstressed 


and undeflected. This will generally result in errors no greater than 5 per- 
cent of the maximum values. 
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DEFLECTION LIMITATIONS OF BRIDGES* 


Progress Report of the Committee on Deflection Limitations 
of Bridges of the Structural Division 


Discussion by A. Kabaila 


A KABAILA, !.—This report is a very timely one. The engineering pro- 
fession holds A.A.S.H.O. Bridge Specifications in high esteem not only in 
U.S.A., but also in other parts of the world. For instance, in 1953 the 
Conference of the State Road Authorities of Australia have adopted A.A.S.H.O. 
Bridge Specifications almost without modification. The Committee therefore 
should be congratulated in their effort to rationalize clause 3.6.10. of this 
Specification. 

One part of clause 3.6.10., which is open to some criticism, seems to have 
been treated by the Committee as a matter of minor importance only: 


“When bridges have cross bracing or diaphragms sufficient in depth 
and strength to insure lateral distribution of loads, the deflection may be 
computed for the standard loading, considering all beams or girders as 
acting together and having equal deflection.” 


Usually plate girder bridges have cross bracing of the same depth as the 
girders. If the bracing members are of little heavier size than the minimum 
specified, then almost all moderate width bridges, except those of rather short 
spans, do fall under the category of “bridges with sufficiently rigid and strong 
cross bracing.” When high strength steel is used and the girder depth is 
somewhat restricted, stiffness considerations are often the critical design 
factor. However, the deflection of all girders is not anywhere nearly the 
same, provided that the live load is placed towards one side of the bridge. 

A numerical example below demonstrates this aspect: 

Figure 1-1 symbolically shows a general arrangement of a bridge, Figure 
1-2 -a typical cross section. All girders are identical. The permissible 
live load deflection at point “C” is 41 x 12 / 300 = 1.64 ins. Live load deflec- 
tion computed on the assumption that all girders deflect equally is equal to 
1.15 ins. (see step 2. below). 

If the deformations compatible with the restraint conditions and statics 
are considered, then the computed live load deflection is 1.72 ins. (see step 4. 
below). In the latter case the following assumptions are made: 


1. Cross frame at “C” is rigid. In this example, this approximation has 
an effect on maximum deflection of approx. 3%. 


a. Proc. Paper 1633, May, 1958. 


1. Lecturer in Civ. Eng., Royal Melbourne Technical College, Victoria, 
Australia, 
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2. Ignore other cross frames. Because of the first assumption this is of 
no consequence on the deflections at point “C.” 


3. Torsional rigidity of the deck and the girders is neglected. The tor- 
sional effect is negligible in this case. 


Following are the steps of the computation: 


1. Consider an imaginary vertical restraint (a fictitious “jack”) at point *C? 
(Figure 2-1). The reaction of the “jack for one lane load, including the 
impact of 17.5%, is 95.7 kips. 


2. Figure 2-2 shows the computed stiffness of one girder: a load of 41.6 k. 
at point “C” produces a deflection of 1. ins. Since the bridge can accom- 
modate 2 lanes of loading only, it follows immediately that, if in accordance 
with 3.6.10 equal deflection for all girders is assumed, the deflection is 
equal to: 95.7 x 2 / 41.6 x 4 = 1.15 ins. 


3. When the “jacks” are removed, the effect is equal to applying a force of 
95.7 k. per lane at point “C.” Cross frame at “C” deflects downwards and 
simultaneously rotates into the position shown in figure 2-3. It follows 
from the geometrical considerations that the deflections of the girders 
from left to right are: z+ 3y; z+y; z-y; and z - 3y (ins.) 
respectively, where z and y are constants to be determined from the 
considerations of equilibrium. 


4. Reactions in kip units on the cross frame from the girders are equal to 
the corresponding deflections in inches multiplied by 41.6. Hence, the 
condition that the sum of vertical forces is equal to zero, gives: 
4zx 41.6 = 191.4k. Therefore z = 1.15 ins. 

Sum of moments about the centre line gives: 
2 (12.75 x 3y + 4.25y) x 41.6 = 191.4 x 3.5 ft.k.; Hence y = 0.19 ins. 
It follows that the maximum deflection of the outer girder is equal to: 


z+ 3y = 1.15 +3 x 0.19 = 1.72 ins. 


and that of the inner girder: 


z+y = 1.15 + 0.19 = 1.34 ins. 


The difference between this computed deflection, which is consistent with 
the conditions of equilibrium, and the assumption that all girders deflect 
equally, is approx. 50%. Little consolation can be found in the argument that 
the live load is not likely to be on one side of the bridge: one truck over- 
taking another certainly would approach the condition assumed in this compu- 
tation and in most cases the truck load is as severe as the lane load. 

Incidentally, the computed deflection of the internal girder, assuming that 
each girder acts independently and using the fractional lane loads given by 
the Specification, is equal to 1.96 ins. 

Where the designs are carried out by the consultants on a competitive 
basis, the main purpose of the specification is to provide the competitors 
with the minimum requirements of the responsible authority. Though the 
limiting values of maximum deflection are rather arbitrary, the variation of. 
50% between the deflections computed on a rational basis and those obtained 
by the procedure outlined by the Specification, is most undesirable. Rather 
than to encourage a competitor to use more accurate rational design methods, 
the Specification penalizes him for doing so. 
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The problem of load distribution is inseparable from the problem of de- 
flections. The live loads transmitted to the girders are, of course, propor- 
tional to the deflections. Thus we find that the portion of one lane load car- 
ried by the outside girder is equal to: 


1.72 x 41.6 / 95.7 = 0.748; or 1.50 of a wheel load. 
And on the inner girder: 
1.34 x 41.6 / 95.7 = 0.583; or 1.17 or a wheel load. 


The fraction of the wheel load on each stringer in accordance with 
A.A.S.H.O. Specification is: 


Outside girder: sixth edition of the Specification 1.00 
seventh edition of the Specification 1.39 
Inner girder: both editions of the Specification 1.70 


In the latter case the live load on the inner girder is overestimated by ap- 
proximately 45%. Undoubtedly, many bridge designing engineers are aware 
of the above conclusions and, in accordance with the general note of 
Division III*, they use fractional loads obtained in a similar manner to this 
numerical example. 

A more definite statement by the Specification of the acceptable approach 
to this problem, however, would be desirable. When the flexibility of the 
cross bracing or the diaphragms is comparable to that of the stringers 
(either because of greater width of more flexible cross members), the prob- 
lem of load distribution is more involved and the difference between the 
values given by the Specification and those computed by a more rigorous 
analysis may not be as great. However, satisfactory methods of analysis 
are readily available.(2,4) To simplify the designer’s task, the specification 
could include design charts and an outline of the recommended solution. In 
the search for practicable, yet more rational design methods, the prepara- 
tion of such charts can hardly be considered as an obstacle. Surely, there 
are many engineers who would consider it a privilege, if a responsible 
authority invited them to undertake such an investigation. 
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¢. Discussion of several papers, grouped by divisions. 
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